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1 


SCOPE 


1.1  Scope .  This  Technical  Report  provides  methods  and 
parameter  values  to  shift  positions  between  approximately  113 
geodetic  datums  and  provides  methods  for  converting  between 
geodetic  coordinates,  Cartesian  coordinates  and  map  projection 
coordinates.  Guidance  is  provided  on  selecting  methods 
appropriate  to  the  application,  and  on  developing  and  testing 
software  to  implement  these  conversions. 

1.2  Applicability.  The  methods  provided  in  this  Technical 
Report  are  used  for  a  wide  range  of  Army  mapping,  charting,  and 
positioning  applications.  Equations  are  furnished  for  the  map 
projections  commonly  used  within  the  Army.  References  are 
provided  for  other  map  projections  that  occasionally  may  be 
encountered.  Special  high  accuracy  applications  such  as 
engineering,  construction,  and  real  estate  boundary  surveys  are 
outside  the  scope  of  this  Technical  Report.  Appendix  B  provides 
datum  shift  parameters  for  approximately  113  datums  that  are 
currently  available  from  the  Defense  Mapping  Agency  (DMA) . 

However,  many  more  local  datums  exist.  For  guidance  on  datum 
shift  applications  or  for  datum  shift  parameters  not  covered  by 
this  Technical  Report,  contact  TEC  [address  given  in  the  preface] . 

1.3  Application  guidance.  Transformation  methods  and 
parameters  should  be  appropriate  to  the  application.  A  system 
analyst  should  do  a  thorough  evaluation  of  the  system  requirements 
before  developing  detailed  specifications.  This  evaluation  should 
consider  the  intended  use,  data  types,  formats,  and  required 
accuracy  of  the  output  mapping,  charting,  and  geodesy  (MC&G)  data, 
as  well  as  the  source  and  accuracy  of  the  input  data.  Also,  any 
hardware  constraints  must  be  considered.  The  methods  presented  in 
this  Technical  Report  can  then  be  examined  for  how  well  they  can 
meet  or  be  adapted  to  meet  the  requirements.  Where  more  than  one 
method  will  satisfy  a  basic  requirement,  give  priority  to  Army- 
preferred  methods  for  standardization  purposes.  Additional 
application  guidance  is  furnished  in  Section  8. 
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2. 


REFERENCED  DOCUMENTS 


2 . 1  Government  documents . 

2.1.1  Specifications,  standards, _ and  handbooks .  The 

following  specifications,  standards,  and  handbooks  form  a  part  of 
this  document  to  the  extent  specified  herein. 

SPECIFICATIONS.  None 

STANDARDS . 

MIL- STD- 24 01  Department  of  Defense  World  Geodetic 

System (WGS)  (1994) 

STANAG  2211  Geodetic  Datums,  ellipsoids,  grids  and 

(Fifth  Edition)  grid  references  (1991) 


HANDBOOKS .  None 

2.1.2  Other  Government  documents,  drawings,  and 
publications .  The  following  other  Government  documents,  drawings,  and 
publications  form  a  part  of  this  document  to  the  extent  specified 
herein.  Unless  otherwise  specified,  the  issues  are  those  cited  in  the 
solicitation.  The  date  in  parentheses  indicates  for  each  document  the 
edition  that  was  used  in  preparation  of  this  Technical  Report . 

Document  ID/Agency  Title 

DMA  Department  of  Defense  Glossary  of 

Mapping,  Charting,  and  Geodetic  Terms, 
Fourth  Edition.  (1981) 


DMA  TM  8358.1  Datums,  Ellipsoids,  Grids,  and 

Grid  Reference  Systems,  (1990) 


DMA  TM  8358.2  The  Universal  Grids:  Universal 

Transverse  Mercator  (UTM)  and 
Universal  Polar  Stereographic  (UPS) . 
(1989) 


2 


DMA  Instruction  8000.1  Geodetic  and  Geophysical  Sign 

with  change  1,  4  Aug  92  Conventions  and  Fundamental  Constants. 

(1988) 


Geological  Survey  Map  Projections  -  A  Working  Manual. 

Professional  Paper  1395  (Snyder,  1987) 

NOAA  Manual  NOS  NGS  5  State  Plane  Coordinate  System  of  1983. 

(Stem,  1989) 

NSWC/DL  TR-3624  Map  Projection  Equations.  (Pearson, 

1977) 

U.S.  Department  of  Conformal  Projections  in  Geodesy  and 

Commerce,  Coast  and  Cartography.  (Thomas,  1979) 

Geodetic  Survey, 

Special  Publication 
No.  251 

Engineer  Technical  Engineering  and  Design  Conversion  to 

Letter  No.  1110-1-147  North  American  Datum  of  1983  (1990) 

Federal  Register  Notice  To  Adopt  Standard  Method  for 

v.  55,  no.  155,  Friday,  Mathematical  Horizontal  Datum 
August  10,  1990  Transformation  (NOAA,  1990) 

Docket  No.  900655-0165 

Copies  of;  DMA  Instruction  8000.1  and  GS  Professional  Paper  1395  are 
available  from  USGS,  Branch  of  Distribution,  Box  25286,  Denver  CO 
80225.  Copies  of  DoD  Glossary  of  Mapping,  Charting,  and  Geodetic 
Terms;  DMA  TM  8358.1;  and  DMA  TM  8358.2  are  available  from  the 
Director,  DMA  Combat  Support  Center,  ATTN;  CCOR,  6001  MacArthur 
Boulevard,  Bethesda  MD  20816-5001.  Copies  of  the  NOAA  Manual  NOS  NGS 
5  and  the  Coast  and  Geodetic  Survey  Special  Publication  No.  251  are 
available  from  National  Geodetic  Information  Branch,  N/CG174,  National 
Geodetic  Survey,  1315  East-West  Highway,  Silver  Spring  MD  20910. 

Copies  of  NSWC/DL  TR  3624  (Accession  Document  Number  =  AD  A037  381  on 
paper,  or  GIDEP  #E151-2353  on  cartridge)  are  available  on  from  Defense 
Technical  Information  Center,  Bldg.  5,  Cameron  Station,  Alexandria,  VA 
22304-6145.  Copies  of  ETL  1110-1-147  are  available  through  the 
Department  of  the  Army,  U.S.  Army  Corps  of  Engineers,  Washington,-  D.C. 
203140-1000.  Copies  of  STANAG  2211  and  MIL- STD  2401  are  available 
from  the  Defense  Printing  Service,  ATTN:  DODSSP,  700  Robbins  Ave,  Bid 
4D  Philadelphia  PA  19111. 
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2.2  Non-Government  publications. 

Heiskanen,  Weikko  A.,  and  Helmut  Moritz;  Physical  Geodesy; 
Institute  of  Physical  Geodesy,  Technical  University,  Graz, 
Austria;  1967. 

Krakiwsky,  Edward  J. ;  Conformal  Map  Projections  in  Geodesy; 
Department  of  Surveying  Engineering,  University  of  New  Brunswick, 
Fredrickton,  N.B.,  Canada;  September,  1973. 

Rapp,  Richard  H.;  Geometric  Geodesy  -  Part  I;  Department  of 
Geodetic  Science  and  Surveying,  The  Ohio  State  University, 
Columbus,  Ohio;  1984. 

Rapp,  Richard  H. ;  Geometric  Geodesy  -  Part  II;  Department  of 
Geodetic  Science  and  Surveying,  The  Ohio  State  University, 
Columbus,  Ohio;  1987. 

2.3  Order  of  precedence.  MIL-STD-2401  shall  take 
precedence  for  all  matters  concerning  WGS  84  and  datum  shift 
parameters.  Nothing  in  this  document,  however,  supersedes 
applicable  laws  and  regulations  unless  a  specific  exemption  has 
been  obtained. 
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3.  DEFINITIONS  AND  UNITS 


з. 1  Acronyms .  The  acronyms  used  in  this  Technical  Report 
are  defined  as  follows: 

a.  BIH  Bureau  International  de  1'Heure 

b.  CONUS  continental  United  States 

c.  CTP  Conventional  Terrestrial  Pole 

d.  DMA  Defense  Mapping  Agency 

e.  DoD  Department  of  Defense 

f.  DTED  Digital  Terrain  Elevation  Data 

g.  ECEF  earth-centered,  earth-fixed 

h.  GPS  Global  Positioning  System 

i.  LGS  local  geodetic  system 

j .  MGRS  Military  Grid  Reference  System 

k.  MC&G  mapping,  charting  and  geodesy 

l.  MRE  multiple  regression  equation 

m.  MSL  mean  sea  level 

n.  NAD  27  North  American  Datum  1927 

o.  NAD  83  North  American  Datum  1983 

p.  NGS  National  Geodetic  Survey 

q.  NGVD  29  National  Geodetic  Vertical  Datum  of  1929 

r.  NOAA  National  Oceanic  and  Atmospheric 

Administration 

s.  NOS  National  Ocean  Survey 

t.  NSWC/DL  Naval  Surface  Weapons  Center/Dahlgren 

Laboratory 

и.  TM  Transverse  Mercator 

v.  UPS  Universal  Polar  Stereographic 

w.  UTM  Universal  Transverse  Mercator 

x.  WGS  72  World  Geodetic  System  1972 

y.  WGS  84  World  Geodetic  System  1984 


3.2  Terms . 


3.2.1  Convergence  of  the  meridian  (y)  .  The  angle  between 
the  projection  of  the  meridian  at  a  point  on  a  map  and  the 
projection  North  axis  (typically  represented  by  the  y  axis  of  the 
projection) .  The  convergence  of  the  meridian  is  positive  in  the 
clockwise  direction. 

3.2.2  Coordinate.  Linear  or  angular  quantities  that 
designate  the  position  that  a  point  occupies  in  a  given  reference 
frame  or  system.  Also  used  as  a  general  term  to  designate  the 
particular  kind  of  reference  frame  or  system,  such  as  Cartesian 
coordinates  or  spherical  coordinates. 

3.2.3  Datum.  Any  numerical  or  geometrical  quantity  or  set 
of  such  quantities  specifying  the  reference  coordinate  system  used 
for  geodetic  control  in  the  calculation  of  coordinates  of  points 
on  the  earth.  Datums  may  be  either  global  or  local  in  extent.  A 
local  datum  defines  a  coordinate  system  that  is  used  only  over  a 
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region  of  limited  extent.  A  global  datum  specifies  the  center  of 
the  reference  ellipsoid  to  be  located  at  the  earth's  center  of 
mass  and  defines  a  coordinate  system  used  for  the  entire  earth. 

3. 2. 3.1  Horizontal  datum.  A  horizontal  datum  specifies  the 
coordinate  system  in  which  latitude  and  longitude  of  points  are 
located.  The  latitude  and  longitude  of  an  initial  point,  the 
azimuth  of  a  line  from  that  point,  and  the  semi -major  axis  and 
flattening  of  the  ellipsoid  that  approximates  the  surface  of  the 
earth  in  the  region  of  interest  define  a  horizontal  datum. 

3. 2. 3. 2  Vertical  datum.  A  vertical  datum  is  the  surface  to 
which  elevations  are  referred.  A  local  vertical  datum  is  a 
continuous  surface,  usually  mean  sea  level,  at  which  elevations 
are  assumed  to  be  zero  throughout  the  area  of  interest . 

3.2.4  Earth- fixed.  Stationary  with  respect  to  the  earth. 

3.2.5  Elevation.  Vertical  distance  measured  along  the  local 
plumb  line  from  a  vertical  datum,  usually  mean  sea  level  or  the 
geoid,  to  a  point  on  the  earth.  Often  used  synonymously  with 
orthometric  height. 

3.2.6  Ellipsoid.  The  surface  generated  by  an  ellipse 
rotating  about  one  of  its  axes.  Also  called  ellipsoid. of 
revolution. 

3.2.7  Equator .  The  line  of  zero  geodetic  latitude;  the 
great  circle  described  by  the  semi -major  axis  of  the  reference 
ellipsoid  as  it  is  rotated  about  the  semi -minor  axis. 

3.2.8  Eouipotential  surface .  A  surface  with  the  same 
potential,  usually  gravitational  potential,  at  every  point;  a 
level  surface . 

3.2.9  Geocentric  Cartesian  coordinates.  Cartesian 
coordinates  (X,  Y,  Z)  that  define  the  position  of  a  point  with 
respect  to  the  center  of  mass  of  the  earth. 

3.2.10  Geodetic  coordinates  (geodetic  position) .  The 
quantities  of  latitude,  longitude,  and  geodetic  height  (0,  X  ,h) 
that  define  the  position  of  a  point  on  the  surface  of  the  earth 
with  respect  to  the  reference  ellipsoid.  Also,  imprecisely, 
called  geographic  coordinates  (See  also  geographic  coordinates) . 

3.2.11  Geodetic  height  (ellipsoidal  height,  h) .  The  height 
above  the  reference  ellipsoid,  measured  along  the  ellipsoidal 
normal  through  the  point  in  question.  The  geodetic  height  is 
positive  if  the  point  is  outside  the  ellipsoid. 
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3.2.12  Geodetic  latitude  ((j))  .  The  angle  between  the  plane 
of  the  Equator  and  the  normal  to  the  ellipsoid  through  the 
computation  point.  Geodetic  latitude  is  positive  north  of  the 
equator  and  negative  south  of  the  Equator  (See  figure  4.2). 

3.2.13  Geodetic  longitude  (X.)  .  The  angle  between  the  plane 
of  a  meridian  and  the  plane  of  the  prime  meridian.  A  longitude 
can  be  measured  from  the  angle  formed  between  the  local  and  prime 
meridians  at  the  pole  of  rotation  of  the  reference  ellipsoid,  or 
by  the  arc  along  the  Equator  intercepted  by  these  meridians.  See 
figure  4.2. 

3.2.14  Geographic  coordinates.  The  quantities  of  latitude 
and  longitude  that  define  the  position  of  a  point  on  the  surface 
of  the  earth.  See  also  geodetic  coordinates. 

3.2.15  Geoid.  The  equipotential  surface  of  the  earth's 
gravity  field  approximated  by  undisturbed  mean  sea  level  of  the 
oceans.  The  direction  of  gravity  passing  through  a  given  point  on 
the  geoid  is  perpendicular  to  this  equipotential  surface. 

3.2.16  Geoid  separation  (N) .  The  distance  between  the  geoid 
and  the  mathematical  reference  ellipsoid  as  measured  along  the 
ellipsoidal  normal.  This  distance  is  positive  outside,  or 
negative  inside,  the  reference  ellipsoid.  Also  called  geoidal 
height;  undulation  of  the  geoid.  See  figure  4.1. 

3.2.17  Grid  reference  system.  A  plane -rectangular 
coordinate  system  usually  based  on,  and  mathematically  adjusted 
to,  a  map  projection  in  order  that  geodetic  positions  (latitudes 
and  longitudes)  may  be  readily  transformed  into  plane  coordinates 
and  the  computations  relating  to  them  may  be  made  by  the  ordinary 
methods  of  plane  surveying  (see  also  3.2.21). 

3.2.18  Map  projection.  A  function  relating  coordinates  of 
points  on  a  curved  surface  (usually  an  ellipsoid  or  sphere)  to 
coordinates  of  points  on  a  plane.  A  map  projection  may  be 
established  by  analytical  computation  or,  less  commonly,  may  be 
constructed  geometrically. 

3.2.19  Map  scale.  The  ratio  between  a  distance  on  a  map  and 
the  corresponding  actual  distance  on  the  earth's  surface. 

3.2.20  Meridian.  A  north-south  reference  line,  particularly 
a  great  circle  through  the  geographical  poles  of  the  earth,  from 
which  longitudes  and  azimuths  are  determined;  or  the  intersection 
of  a  plane  forming  a  great  circle  that  contains  both  geographic 
poles  of  the  earth,  and  the  ellipsoid. 

3.2.21  Military  Grid  Reference  System  (MGRS ) .  MGRS  is  a 
system  developed  by  DMA  for  usage  with  UTM  and  UPS  projections. 
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3.2.22  Orthometric  height.  See  elevation. 


3.2.23  Parallel .  A  line  on  the  earth,  or  a  representation 
thereof,  that  represents  the  same  latitude  at  every  point. 

3.2.24  Prime  (initial)  meridian.  A  meridian  from  which  the 
longitudes  of  all  other  meridians  are  reckoned.  This  meridian,  of 
longitude  0°,  was  traditionally  chosen  to  pass  through  the 
Greenwich  Observatory  in  Greenwich,  England.  For  new  refined 
coordinate  systems,  the  location  of  the  prime  meridian  is  defined 
by  the  International  earth  Rotation  Service,  Paris,  France. 

3.2.25  Reference  ellipsoid.  An  ellipsoid  whose  dimensions 
closely  approach  the  dimensions  of  the  geoid;  the  exact  dimensions 
are  determined  by  various  considerations  of  the  section  of  the 
earth's  surface  concerned.  Usually  a  bi-axial  ellipsoid  of 
revolution. 

3.2.26  Scale  factor  (projection) .  A  multiplier  for  reducing 
a  distance  in  a  map  projection  to  the  actual  distance  on  the 
chosen  reference  ellipsoid. 

3.3  Units .  In  this  Technical  Report  all  distances  are 
expressed  in  meters.  Angles  are  expressed  in  degrees,  radians  or 
the  combination  of  degrees,  minutes,  and  seconds. 

3.4  Sign  conventions.  The  sign  conventions  used  in  this 
Technical  Report  are  defined  as  follows: 

a.  Elevation.  Positive  outside  or  above  the  vertical  datum, 
negative  inside  or  below  the  datum. 

b.  Geodetic  height.  Positive  outside  the  reference 
ellipsoid,  negative  inside  the  reference  ellipsoid. 

c.  Geodetic  latitude.  Positive  in  the  northern  hemisphere, 
negative  in  the  southern  hemisphere. 

d.  Geodetic  longitude.  Zero  at  the  prime  meridian,  positive 
0°  to  180°,  eastward  from  the  prime  meridian,  and  negative 
0°  to  180°,  westward  from  the  prime  meridian. 

e.  Geoid  separation.  Positive  outside  the  reference 
ellipsoid,  negative  inside  the  reference  ellipsoid. 

f.  Cartesian  Z.  Positive  in  the  northern  hemisphere. 

3.5  Unit  conversion  factors.  The  standard  unit  of  distance 
measurement  is  the  meter.  However,  MC&G  products  in  many  parts  of 
the  world  were  produced  using  other  units  of  measurement.  The 
user  must  ascertain  the  units  used  in  non-standard  products  and 
appropriately  scale  them  for  use  in  the  equations  given  in  this 
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technical  report.  Some  commonly  used  conversion  factors  are 
listed  in  Table  3.1. 


TABLE  3 . 1  Conversion  Factors . 


1  meter 

1  International  foot 
1  U.S.  Survey  foot 

1  Int ’ 1  Nautical  mile 

1  Int ' 1  Statute  mile 


3.28083333333  U.S.  Survey  feet 
3.28083989501  International  feet 
0.3048  meter  (Exact) 

1200/3937  meter  (Exact) 
0.30480060960  meter 
1852  meters  (Exact) 

6076.10333333  U.S.  Survey  feet 
6076.11548556  International  feet 
1609.344  meters  (Exact) 

5280  International  feet  (Exact) 


1  degree 
1  degree 
1  minute 
6400  U.S.  mils 


=  60  minutes 
=  3600  seconds 
=  60  seconds 
=  360  degrees 


3.5.1  Decrees  and  radians .  A  circle  contains  360  degrees  or 
27T  radians.  If  an  angle  contains  d  degrees  or,  equivalently,  r 
radians,  then  d  and  r  are  related  by 


_  _7Cd_ 
r  “  180 


(3.1) 


3.5.2  Specifying  the  unit  of  angular  measure.  In  this 
Technical  Report  the  unit  of  angular  measure  will  be  specified  if 
it  matters.  A  general  rule  of  thumb  is  that  angles  must  be  in 
radians  in  an  equation  that  was  derived  using  calculus  and 
trigonometry.  Radians  are  the  usual  units  in  computer  evaluation 
of  trigonometric  functions. 
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4 .  FUNDAMENTAL  CONCEPTS 

4.1  Introduction.  An  understanding  of  geodesy  helps  one 
understand  how  to  represent  the  locations  of  points  on  or  about 
the  surface  of  the  earth.  Terms  printed  in  boldface  are  defined 
in  Section  3.2. 

a.  The  location  of  a  point  is  described  by  a  set  of 
coordinates.  Some  types  of  coordinates,  such  as 
Cartesian  (see  Section  4.3.1),  uniquely  identify 
three-dimensional  locations  of  points  in  space.  Their 
use  is  computationally  efficient  in  many  applications. 
However,  users  often  have  difficulty  visualizing  the 
relationship  of  Cartesian  coordinates  to  the  surface  of 
the  earth.  Geodetic  coordinates  (see  Section 
4.3.3)  make  visualization  easier  by  providing 
horizontal  coordinates  on  the  surface  of  a  figure 
approximating  the  shape  of  the  earth  and  a  vertical 
(height)  coordinate  above  or  below  the  reference 
surface. 

b.  Geodetic  coordinates  are  defined  by  a  geodetic 
reference  system.  A  geodetic  reference  system 
includes,  but  is  not  limited  to,  a  reference  surface 
(usually  an  ellipsoid)  and  a  set  of  parameters 
positioning  and  orienting  the  reference  surface  with 
respect  to  the  earth.  Over  the  years,  many  geodetic 
reference  systems  have  been  developed  for  various 
localities.  These  systems  are  often  called  local  datums 
and  are  not  rigorously  tied  to  each  other. 

c.  When  applications  require  that  two  or  more  points  be 
located  in  a  common  reference  system,  transformation 
between  two  coordinate  systems  may  be  needed.  For 
example:  both  weapon  and  target  must  be  located  in  the 
same  reference  system  if  the  indirect  fire  computations 
are  to  be  valid.  As  the  ranges  of  weapons  and  target 
acquisition  systems  increased  and  positional  data  were 
derived  from  multiple  sources,  commanders  often  were 
confronted  with  the  problem  of  using  positions  on  more 
than  one  reference  system.  With  the  advent  of 
satellite  geodesy  one  can  establish  global  geodetic 
reference  systems  and  take  measurements  that 

allow  local  datums  to  be  related  to  the  global 
system.  The  standard  global  reference  system  for 
DoD  is  the  World  Geodetic  System  1984  (WGS  84) .  The 
process  of  relating  one  reference  system  to  another  is 
called  a  datum  shift,  or  datum  transformation. 
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d.  Coordinates  generally  are  represented  alpha - 

numerically  (e.g.  latitude  38°  51'  22.45"N  and 
longitude  77°  26*  34.56“W) .  Position  data  for  an  area 
often  are  represented  graphically,  as  on  a  map.  The 
planar  representation  on  a  map  provides  a  distorted 
view  of  features  on  the  curved  earth's  surface. 
Geodesists  have  developed  various  mathematical  methods 
to  minimize  certain  types  of  map  distortions,  depending 
on  the  intended  use  of  the  map.  The  process  of 
relating  coordinates  of  points  on  a  curved  surface 
(usually  an  ellipsoid  or  sphere)  to  coordinates  of 
points  on  a  plane  is  called  map  projection.  The 
preferred  map  projections  for  DoD  are  the  Universal 
Transverse  Mercator  (UTM) ,  used  over  most  of  the 
landmass,  and  Universal  Polar  Stereographic  (UPS) ,  used 
in  the  polar  regions. 

4.2  Reference  surfaces.  Reference  surfaces  are  used  as  the 
bases  for  geodetic  coordinate  systems  and  elevations. 

4.2.1  Reference  ellipsoid.  A  reference  surface  called  the 
reference  ellipsoid  (mathematically  described  in  Section  4.3.2) 
is  used  for  defining  geodetic  coordinates  (Section  4.3.3)  because 
it  is  regularly  shaped,  is  mathematically  tractable,  and  can 
approximate  the  surface  of  the  earth.  Since  the  earth  is 
flattened  at  the  poles,  an  oblate  ellipsoid  fits  the  shape  of  the 
earth  better  than  does  a  sphere.  An  ellipsoid  of  revolution  is 
the  reference  surface  most  commonly  used  in  geodesy  and  is  used 
throughout  this  Technical  Report.  A  geodetic  coordinate  has  a 
simply  defined  relationship  with  the  reference  ellipsoid  and  any 
other  geodetic  coordinate  because  of  the  regularity  and  relative 
simplicity  of  the  reference  ellipsoid's  shape.  The  parameters  for 
the  commonly  used  reference  ellipsoids  are  listed  in  Appendix  A. 

A  reference  ellipsoid  is  often  chosen  to  fit  the  geoid. 

4.2.2  Geoid.  The  geoid  is  a  closed  surface  of  constant 
gravity,  potential  approximated  by  mean  sea  level  (MSL)  and  the 
theoretical  extension  of  MSL  through  land  areas.  For  purposes  of 
this  technical  report,  the  geoid  is  considered  a  reference  surface 
from  which  elevations  are  measured.  A  reference  ellipsoid  may  fit 
the  geoid;  for  example,  the  maximum  difference  between  the  WGS  84 
geoid  model  and  the  WGS  84  ellipsoid  found  in  Appendix  D,  Table 
D.l  is  102  meters. 


4.2.3  Relationships  among  topography,  the  geoid.  and  the 
reference  ellipsoid.  The  relationships  among  topography,  the 
geoid,  and  the  reference  ellipsoid  are  illustrated  in  figure  4.1. 
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Figure  4.1  The  relationship  between  the  reference  ellipsoid,  the  geoid,  and  the  physical 
surface  of  the  earth.  Note  relationship  of  elevation  (H) ,  geodetic  height  (h) 
and  geoid  separation  (N) . 


4.3  Earth-fixed  coordinate  systems.  There  are  many 
coordinate  systems  used  to  represent  positions  on  or  near  the 
earth.  The  coordinate  systems  used  in  this  Technical  Report  are 
all  earth-fixed.  A  coordinate  system  is  earth-fixed  if  the  axes 
are  stationary  with  respect  to  the  rotating  earth.  Two  primary, 
earth-fixed  coordinate  systems,  Cartesian  and  geodetic,  are 
discussed  in  Sections  4.3.1  and  4.3.3,  respectively.  Conversions 
from  one  to  the  other  are  presented  in  Section  4.3.4. 

4.3.1  Cartesian  coordinates  (X,  Y,  Z)  .  The  Cartesian 
coordinate  system  used  in  this  Technical  Report  is  a  right-hand, 
rectangular,  three-dimensional  (X,  Y,  Z)  coordinate  system  with 
an  origin  at  (0,0,0).  When  the  origin  is  also  at  the  mass  center 
of  the  earth,  the  coordinate  system  is  said  to  be  geocentric. 


4.3.2  Cartesian  cooydinate_sy.s.t£in/re£grencg_g.liipg.g.id 
relationship.  The  relationship  of  the  Cartesian  coordinate  system 
to  the  reference  ellipsoid  is  shown  in  figure  4.2.  The  Z-axis, 
the  axis  of  rotation  (semi-minor  axis)  of  the  ellipsoid,  is  nearly 
parallel  to  the  axis  of  rotation  of  the  earth.  The  Z-coordinate 
is  positive  toward  the  North  pole.  The  X-Y  plane  lies  in  the 
equatorial  plane,  the  plane  swept  out  by  the  semi -major  axis  as 
the  ellipse  is  rotated.  The  X-axis  lies  along  the  intersection  of 
the  plane  containing  the  prime  (initial)  meridian  and  the 
equatorial  plane.  The  X-coordinate  is  positive  toward  the 
intersection  of  the  prime  meridian  and  equator.  The  reference 
ellipsoid  satisfies  the  equation 


2^  •  X! 

a2  +  a2  +  b2  “  1 


(4.1) 


Here,  a  and  b,.  respectively,  denote  the  semi-major  and  semi-minor 
axes  of  the  ellipsoid.  This  surface  can  be  generated  by  rotating 
an  ellipse  defined  by 


X2  zi_ 

a2  +  b2  " 


(4.2) 


about  the  Z-axis.  Ideally,  for  global  applications,  the  Z-axis 
coincides  with  the  mean  axis  of  rotation  of  the  earth  and  the 
origin  is  the  center  of  mass  of  the  earth.  Often,  in  practice,  an 
ellipsoidal  surface  is  fitted  to  a  local  arfea  so  that  the  origin 
and  axes  vary  slightly  from  ideal  coincidence. 
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Figure  4.2  The  geometric  relationship  between  Cartesian 
and  geodetic  coordinates. 
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4. 3. 2.1  First  eccentricity  and  flattening.  Two  important 
constants  are  the  (first)  eccentricity  8  and  flattening  f  that  are 
related  to  a  and  b  by 


Frequently  a  and  f  are  given,  rather  than  a  and  b,  as  the 
defining  parameters  of  a  reference  ellipsoid 

4. 3. 2. 2  Ellipsoid  parameters.  Appendix  A  contains  a  list  of 
reference  ellipsoids  and  their  parameters.  These  can  be  divided 
into  two  categories:  those,  such  as  the  Bessel  Ellipsoid,  that  are 
used  to  approximate  the  shape  of  the  earth  for  a  local  area  and 
those,  such  as  the  WGS  84  Ellipsoid,  that  are  used  to  yield  a 
global  approximation  to  the  shape  of  the  earth. 

4.3.3  Geodetic  coordinates  (<|),  A,  h)  .  Geodetic  coordinates 
consist  of  geodetic  latitude  ((|)) ,  geodetic  longitude  (A) ,  and 
geodetic  height  (h) .  Their  relationship  to  the  reference 
ellipsoid  is  shown  in  figure  4.2.  For  an  ellipsoid  that  satisfies 
equation  4.1,  the  normal  SP  intersects  the  ellipsoid  at  Q.  The 
angle  between  the  normal  SP  and  the  equatorial  (X-Y)  plane  is 
called  the  geodetic  latitude  ((}))  of  point  P.  The  meridian  plane 
containing  point  P  is  defined  as  the  half-plane  containing  the  Z- 
axis  and  point  P.  The  angle  between  the  prime  meridian  (X-Z) 
plane  and  the  meridian  containing  point  P  is  the  geodetic 
longitude  (A)  of  point  P.  Geodetic  longitude  is  not  defined  when 
P  lies  on  the  Z-axis.  The  distance  from  Q  to  P  is  called  the 
geodetic  height  (h)  . 

4. 3. 3.1  Latitude  and  longitude  limits. 

Latitude  limits: 

-90°  to  +90° 

-90°  at  the  South  Pole 
0°  at  the  Equator 
+90°  at  the  North  Pole 

Longitude  limits: 

0°  at  the  Prime  meridian 
180°W  (-180°)  to  180°E  (+180°) 

4.3.4  Coordinate  conversion.  Two  basic  earth- fixed 
coordinate  systems,  geocentric  Cartesian  coordinates  and  geodetic 
coordinates,  have  been  introduced  in  Sections  4.3.1  and  4.3.3. 

The  following  section  explains  how  coordinates  in  one  system  can 
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be  transformed  to  the  other.  In  Section  6,  transformations 
between  coordinates  on  two  different  reference  ellipsoids  are 
discussed.  Map  projections  are  treated  in  Section  7  and  the 
general  conversion  process  is  summarized  in  Section  8.1. 

4. 3. 4.1  Geodetic  to  Cartesian  coordinate  conversion.  If 
geodetic  coordinates  ((J),  X,  h)  are  known,  then  the  Cartesian 
coordinates  are  given  by 

X  =  (Rn  +  h)  cos  <J>  cos  X  (4.4) 

Y  =  (Rn  +  h)  cos  <|)  sin  X 

Z  =  (fi  Rn  +  sin 

where  RN,  the  radius  of  curvature  in  the  prime  vertical,  is  given 
by 


Rn  - 


a2 


Va2"  cos2  <j)  +  b2  sin2  <|) 


(4.5) 


_  _ a _ 

Vl  -  £2  sin2  (J) 

See  Figure  4.2.  Equation  sets  4.4  and  4.5  are  exact.  Calculating 
h  from  elevation  (H)  is  discussed  on  section  6.5. 


4. 3. 4. 2  Cartesian  to  geodetic  coordinate  conversion. 


4. 3. 4. 2.1  Finding  X.  In  the  reverse  case,  when  (X,  Y,  Z) 
are  known  but  (<j)  , X,  h)  are  unknown. 


for 

X  *  0,all  Y, 

X  =  0, 

X  =  0, 

X  =  0, 


the  longitude  (X)  is  given  bv 

Y 

X  =  arctan  —  (4.6) 

Y  >  0,  X  =  90°  E  (+90°) 

Y  <  0,  X  =  90°  W  (-90°  or  270°  E) 

Y  =  0,  X  is  undefined 


Equation  set  4.6  is  exact.  Take  care  to  ensure  that  the 
evaluation  of  the  arctan  function  yields  a  value  of  X  lying  in  the 
desired  quadrant. 


4. 3. 4. 2. 2  Finding  (j).  When  (X,  Y)  is  not  (0,  0),  the  general 
equations  for  <j>  and  h  are  more  involved.  The  method  presented  for 
finding  <|>  is  an  iterative  procedure  of  Bowring.  Following  Rapp 
(1984,  pp.  123-124),  one  finds  an  initial  approximation  to  a 
variable  (3  given  by 
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tan  p0  = 


bVx2  +  Y2 


Once  P0  is  obtained,  substitute  it  for  P  in  the  equation 


(4.7) 


tan  <})  = 


Z  +  8' 2  b  sin3 


/x2  +  Y2  -  a  £2  cos3  P 


Here  £2  and  £'2  are  given  by 


_  J)2 

£2  =  -  =  2f  -  f2 

a2 


E'2  = 


a2  -  b2 


£2 

1  -  E2 


(4.8) 


(4.9) 


This  approximation  to  <J>  is  substituted  into 

tan  P  =  (1  -  f)  tan  (j) 


(4.10) 


to  give  an  updated  approximation  to  p.  This  procedure,  of  using 
the  latest  approximation  of  p  (equation  4.10)  to  produce  an 
updated  approximation  of  (|)  (equation  4.8),  and  then  using  this  new 
value  of  (|)  to  update  P,  can  be  continued  until  the  updated  value 
of  <|)  is  sufficiently  close  to  the  previous  value.  For  many 
purposes,  one  iteration  is  sufficient  when  using  the  Bowring 
method,  so  that  (j)  can  be  found  directly  from  equation  4.8,  using 
the  original  value  of  P0.  Rapp  (1984,  p.  124)  states  that  for 
terrestrial  applications,  one  such  iteration  yields  values 
accurate  to  within  0.1  millimeter,  in  the  polar  case  (X  =  0, 

Y  =  0,  Z),  the  latitude  (<|))  is  defined  to  be  +90°  (90°N)  for 
positive  Z  values  and  -90°  (90°S)  for  negative  Z  values. 

4. 3. 4. 2. 3  Calculating  h.  Once  <J>  has  been  calculated,  in 
non-polar  areas  h  can  be  found  from 


Vx2  +  Y2  _ 

h  = - r~  -  Rn 

cos  <(> 


(4.11) 


In  polar  regions  it  is  preferable  to  use 


Rn  +  £2  Rn 


(4.12) 


where  RN  is  defined  in  equation  4.5 
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Note:  For  many  applications  h  is  not  required  and  need  not  be 
calculated.  An  exception  to  this  statement  is  the  computation  of 
height  or  elevation  for  the  Global  Positioning  System  (GPS) . 

4.4  Representation  of  Geodetic  Coordinates  Latitude  is 
given  before  longitude.  Latitude  and  longitude  usually  are 
numerically  represented  as:  decimal  degrees  (e.g.  44.44°  N); 
integer  degrees  and  decimal  minutes  (e.g. 44°  44.44'S);  or  integer 
degrees,  integer  minutes  and  decimal  seconds  (e.g.  44°  44'  44.44” 

E) .  For  clarity,  numeric  fields  representing  degrees,  minutes  and 
seconds  should  be  ended  by  the  appropriate  symbol,  respectively 

"  °  ’»  "  '  ",  and  Numeric  fields  with  a  value  of  zero  should  be 

indicated  by  a  "OO” ,  and  not  left  blank.  The  hemisphere  of  both 
latitude  and  longitude  must  be  shown,  preferably  by  placing  an  "N" 
or  “S"  at  the  end  of  numeric  latitude  and  an  "E"  or  "W"  after 
numeric  longitude. 

4.5  Height  Relationships. 

4.5.1  Elevations.  Throughout  this  Technical  Report,  the 
term  elevation  (H)  is  used  to  denote  the  distance  of  a  point 
above  the  geoid  or  vertical  reference  surface  as  measured  along 
the  plumb  line.  A  plumb  line  follows  the  direction  of  gravity  and 
is  perpendicular  to  all  equipotential  surfaces  of  the  Earth's 
gravity  field  that  intersect  it.  Points  lying  outside  the  geoid 
are  defined  as  having  positive  elevation.  For  the  purposes  of 
this  Techncial  Report,  elevation  and  orthoxnetric  height  are 
considered  equivalent.  In  practice,  the  reference  surface  for 
measuring  elevations  may  not  exactly  coincide  with  the  geoid. 

4.5.2  Geoid  separation.  Geoid  separation  (N)  is  the 
distance  from  a  reference  ellipsoid  to  the  geoid,  measured  along 
the  ellipsoidal  normal.  This  is  equivalent  to  the  geodetic  height 
(see  Section  4.3.3)  of  a  point  on  the  geoid.  Geoid  separation  is 
positive  when  the  geoid  lies  outside  the  ellipsoid. 

4.5.3  The  relationships  among  H.  h,  and  N.  Geodetic  height 
(h) ,  geoid  separation  (N) ,  and  elevation  (H)  are  depicted  in 
figure  4.1  and  are  related  by  equation  4.13. 

h  «  H  +  N  (4.13) 

4. 5. 3.1  Notational  use  of  H  and  h.  This  Techncial  Report 
uses  the  classical  notation  of  H  for  elevation  and  h  for  geodetic 
height . 
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5.  GEODETIC  SYSTEMS  AND  DATUMS 


5.1  Introduction.  A  geodetic  system  serves  as  the 
framework  for  determining  coordinates  of  points  with  respect  to 
the  earth.  Modem  global  geodetic  systems,  such  as  the  World 
Geodetic  System  1984  (WGS  84),  (see  Section  5.4),  have  been 
established  using  techniques  of  satellite  geodesy.  Typically, 
they  are  defined  by  a  geocentric  Cartesian  coordinate  system  with 
the  Z-axis  along  the  mean  rotation  axis  of  the  earth  and  the  X- 
axis  through  a  longitude  reference. 

a.  DoD  World  Geodetic  Systems  have  an  associated  reference 
ellipsoid  and  a  geoid  model.  The  geoid  model  is  sometimes 
used  as  a  reference  surface  for  elevations.  World 
Geodetic  Systems  provide  a  consistent  framework  for 
determining  geodetic  positions. 

b.  Local  geodetic  systems  are  established  by  a  variety  of 
techniques.  Typically,  they  consist  of  a  geodetic  datum 
(see  Section  5.2)  for  horizontal  control  and  a  vertical 
datum  (see  Section  5.3)  for  elevations.  Any  given  local 
geodetic  system  covers  only  a  small  fraction  of  the 
earth's  surface.  Local  geodetic  datums  and  vertical 
datums  generally  are  independent  of  each  other. 

5.2  Geodetic  (horizontal)  datums.  A  geodetic  datum  has  a 
earth-fixed  reference  ellipsoid  that  may  have  been  fit,  in  some 
manner,  to  the  surface  of  the  earth  in  the  area  of  interest. 
Geodetic  datums  have  parameters  that  define  the  size  and  shape  of 
the  ellipsoid. 

5.2.1  Background .  Historically,  local  geodetic  datums  have 
an  origin,  on  the  surface  of  the  earth,  that  relates  the  geodetic 
coordinate  system  to  the  ellipsoid.  For  these  local  datums,  the 
ellipsoid  semi-minor  axis  generally  does  not  coincide  with  the 
earth's  mean  rotation  axis.  However,  the  North  American  Datum 
1983  (NAD  83)  was  developed  using  the  same  satellite  geodesy 
techniques  used  for  global  geodetic  systems.  This  modern  datum 
has  an  origin  at  the  center  of  the  reference  ellipsoid,  and  the 
semi-minor  axis  coincides  with  the  mean  rotation  axis  of  the 
earth,  without  reference  to  the  word  vertical,  the  word  datum,  as 
used  in  this  Technical  Report,  means  geodetic  (horizontal)  datum. 

5.2.2  Geodetic  datums  and  WGS  84.  Over  the  years,  hundreds 
of  geodetic  datums  have  been  developed  for  various  locations. 
Currently,  113  local  geodetic  datums  have  been  tied  to  WGS  84. 
Coordinate  transformations  from  local  datums  to  WGS  84  are 
discussed  in  Section  6. 
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•-  5.3  Vertical  datums  and  elevations.  By  definition,  a 

vertical  datum  is  a  surface  of  zero  elevation.  Elevations  are 
measured  (positive  upward)  from  the  vertical  datum.  Ideally,  a 
vertical  datum  would  closely  approximate  the  geoid. 

5.3.1  Background.  As  a  practical  matter,  it  is  inpossible 
to  access  the  geoid  surface  directly  for  use  as  a  vertical  datum. 
Historically,  tide  gage  measurements  were  averaged  over  many  years 
to  establish  the  local  mean  sea  level  (MSL)  references  for 
vertical  datums.  The  National  Geodetic  Vertical  Datum  of  1929 
(NGVD  29),  which  is  based  on  tide  gage  measurements  and  precise 
leveling  surveys,  is  estimated  to  be  within  a  few  meters  of  the 
geoid.  There  is  greater  uncertainty  in  the  relationships  between 
the  geoid  and  other  local  vertical  datums. 

5.3.2  Vertical  datums  and  mean  sea  level  (MSL) .  The  exact 
relationships  between  different  vertical  datums  are  unknown. 
Therefore,  when  shifting  datums,  all  elevations  are  considered 
referenced  to  MSL  regardless  of  the  vertical  datum  used  to 
establish  them.  For  Army  applications,  two  exceptions  to  this 
rule  are : 


a.  Conversion  of  elevations  derived  from  WGS  72  geodetic 
height  measurements  to  WGS  84  coordinates  (see  Section 
6.7.1)  . 

b.  Civil  works,  engineering  and  construction  applications 
are  outside  the  scope  of  this  Technical  Report. 

5.3.3  Relationship  between  local  vertical  and  horizontal 
datums .  Generally,  vertical  datums  and  horizontal  datums  are 
independently  defined.  For  example,  NGVD  29  is  the  major  vertical 
datum,  while  the  North  American  Datum  1927  (NAD  27)  and  NAD  83  are 
the  major  horizontal  datums  used  in  North  America. 

5.4  World  Geodetic  System  1984  (WGS  84) .  A  world  geodetic 
system  provides  the  basic  reference  frame,  geometric  figure  and 
gravimetric  model  for  the  earth,  and  provides  the  means  for 
relating  positions  on  various  local  geodetic  systems  to  an  earth- 
centered,  earth-fixed  (ECEF)  coordinate  system. 

a.  WGS  84  is  the  current  standard  DoD  geodetic  system.  WGS 
84  is  the  latest  in  a  series  of  DoD-developed  world 
geodetic  systems.  It  is  a  replacement  for  WGS  72,  the 
previous  DoD  standard  ECEF  world  geodetic  system. 

b.  The  origin  of  the  WGS  84  coordinate  system  is  the  center 
of  mass  of  the  earth.  The  WGS  84  Z-axis  is  parallel  to 
the  direction  of  the  Conventional  Terrestrial  Pole  (CTP) , 
as  defined  by  the  Bureau  International  de  l'Heure  (BIH) . 
The  WGS  84  prime  meridian  is  parallel  to  the  BIH  Zero 
meridian.  The  X-axis  is  the  intersection  of  the  plane  of 
the  WGS  84  reference  meridian  and  the  CTP  equatorial 
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plane.  The  epoch  is  1984.0.  The  Y-axis  completes  a 
right-handed,  earth- fixed  Cartesian  coordinate  system. 

c.  The  WGS  84  reference  ellipsoid  is  a  geocentric  ellipsoid 
of  revolution.  The  WGS  84  coordinate  system  origin  and 
axes  also  serve  as  the  geometric  center  and  the  X,  Y,  and 
Z  axes  of  the  WGS  84  ellipsoid.  WGS  84  ellipsoid 
parameters  are  given  in  Appendix  A. 

d.  The  WGS  84  geoid  model  can  be  used  to  convert  between 
elevations  and  geodetic  heights  (see  Section  5.5). 

e.  Datum  shift  techniques  and  parameters  needed  to  relate 
positions  on  many  local  datums  to  WGS  84  were  also 
developed  as  part  of  WGS  84  (see  Section  6  and  Appendices 
B  and  E) . 

f.  Only  those  portions  of  WGS  84  related  to  the  purpose  of 
this  Technical  Report  were  described  above.  A 
comprehensive  view  of  WGS  84  can  be  obtained  from  DMA  TRs 
8350.2  (1991),  8350. 2 -A  (1987),  and  8350. 2-B  (1987). 

5.5  WGS  84  referenced  elevations.  Sometimes  elevations 
must  be  generated  from  geodetic  heights,  without  reference  to  a 
local  vertical  datum.  For  example,  satellite  positioning  systems, 
like  the  Global  Positioning  System  (GPS),  establish  coordinates 
geometrically  and  cannot  directly  measure  elevations.  An 
elevation  can  be  obtained  from  a  WGS  84  geodetic  height  using: 

H  88  ^WGS  84  “  ^WGS  84  (5.1) 

where  NWGS  84  is  the  value  of  WGS  84  geoid  separation  for  the  WGS  84 
geodetic  position  of  the  point.  Appendix  D  contains  a  method  for 
interpolating  NWGS  84  from  a  table  of  grided  values.  Table  D.l  is  a 
list  of  Nwgs  84  on  a  10°  x  10°  grid.  It  is  useful  for  moderate 
accuracy  mapping  and  charting  applications.  TEC  will  furnish 
qualified  users  with  either  digital  30'  x  30'  or  1°  x  1°  NWGS  84 
tables  or  spherical  harmonic  expansion  tables  for  survey  and  high 
accuracy  applications.  Requests  should  be  sent  to  the  address  in 
the  preface. 
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6. 


DATUM  SHIFTS 


6.1  Introduction.  Several  methods  are  available  to  shift 
coordinates  from  one  geodetic  datum  to  another.  To  develop  these 
shifts,  coordinates  of  one  or  more  physical  locations  must  be 
known  on  both  datums.  Typically/  for  the  local  datum  to  WGS  84 
shifts  presented  in  this  Technical  Report,  the  WGS  84  coordinates 
of  local  datum  survey  control  points  were  determined  using  Doppler 
satellite  observations .  Two  classes  of  datum  shift  methods  are 
geometric  transformation  and  polynomial  fitting  techniques.  Only 
geometric  transformations  are  discussed  here. 

a.  The  generalized  geometric  transformation  model  assumes 
the  origins  of  the  two  coordinate  systems  are  offset  from 
each  other;  the  axes  are  not  parallel;  and  there  is  a 
scale  difference  between  the  two  datums .  This  seven- 
parameter  model  is  discussed  in  Section  6.2.  Data  from 
at  least  three  well-spaced  positions  are  needed  to  derive 
a  seven -parameter  geometric  transformation. 

b.  The  number  of  parameters  in  a  geometric  transformation 
model  can  be  reduced.  The  commonly  used  three- 
parameter  model  neglects  rotations  between  coordinate 
systems  axes  and  scale  differences  between  datums.  Only 
the  origins  are  offset.  The  three-parameter  model  is 
discussed  in  Section  6.3.  Data  from  one  or  more 
positions  are  required  to  derive  a  three-parameter 
geometric  transformation. 

c.  The  WGS  72  to  WGS  84  transformation,  given  in  Section 

6.4,  is  a  special  five-parameter  geometric 
transformation.  When  WGS  84  superseded  WGS  72  as  the 
DoD-preferred  reference  system,  this  transformation  was 
created  specifically  to  transform  WGS  72  coordinates  to 
WGS  84  coordinates. 

d.  Heights  used  for  datum  shifting  are  treated  in  section 

6.5.  Shifting  of  vertical  datums  is  discussed  in  Section 
6.7. 


6.1.1  Shifting  between  two  local  datums.  The  datum 
shift  methods  given  in  this  Technical  Report  are  for  shifting  from 
a  local  datum  to  WGS  84  and  from  WGS  84  to  a  local  datum.  When 
shifts  are  needed  between  two  local  datums,  shift  the  coordinates 
on  the  first  local  datum  to  WGS  84,  then  shift  the  WGS  84 
coordinates  to  the  second  local  datum. 

6. 1.1.1  Shifting  from  NAD  27  to  NAD  83 . 

The  software  package  NADCON  has  been  recommended  as  the  standard 
method  for  transformations  between  the  North  American  Datum  of 
1927  (NAD27 )  and  the  North  American  Datum  of  1983  (NAD  83)  by  The 
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Federal  Geodetic  Control  Committee  (Federal  Register.)  Corpscon, 
a  package  incorporating  NADCON,  is  mandated  for  the  Army  Corps  of 
Engineers  (ETL  1110-1-147.)  Because  NADCON  is  a  survey  tool  it  is 
beyond  the  scope  of  this  Technical  Report.  Questions  can  be 
addressed  to  the  Survey  Division  of  TEC  at  the  address  in  the 
preface,  (703)355-2798,  or  to  the  National  Geodetic  Survey  at 
(301)713-3178. 

6.2  Seven-parameter  geometric  transformation.  The  seven- 
parameter  transformation,  the  most  general  transformation  between 
local  and  global  Cartesian  coordinates,  is  considered  in 
DMA  TR  8350. 2 -A  (1987) .  The  underlying  assumption  for  this  seven- 
parameter  model  is  that  the  local  and  WGS  84  Cartesian  axes  differ 
by  seven  parameters:  three  rotational  parameters  (00,  S,  \| /) ,  a 

scale  change  AS,  and  three  origin  shift  parameters  (AX,  AY,  AZ)  . 
See  Figure  6.1. 


FIGURE  6.1.  Relationship  between  coordinate  axes  in  a 
seven -parameter  geometric  transformation. 

6.2.1  Transformation  to  WGS  84  Cartesian  coordinates.  For 
the  seven-parameter  model,  the  local  and  WGS  84  Cartesian 
coordinates  are  related  by 


^WGS 

84  =  ^LGS 

+ 

AX  + 

®Ylgs  Y^lgs 

+ 

ASX^ 

^WGS 

84  =  Ylgs 

+ 

AY  - 

®^LGS  +  £  ^LGS 

+ 

asy^. 

^WGS 

84  =  ^*LGS 

+ 

AZ  + 

V^LGS  ~  £  ^LGS 

4* 

ASZLGs 
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The  subscript  LGS  designates  Local  Geodetic  System  coordinates  and 
the  subscript  WGS  84  designates  World  Geodetic  System  1984 
coordinates . 

6.2.2  Transformation  to  LGS  Cartesian  coordinates.  The 
equations  for  converting  WGS  84  geocentric  Cartesian  coordinates 
to  local  Cartesian  coordinates  are 


XLgs  =  ^WGS  84  “  AX  "  ®^WGS  84  +  V^WGS  84  ”  ASXwqs  84  (6.2) 

^LGS  =  ^WGS  84  -  AY  +  COXWGS  84  -  ®  “^WGS  84  ~  ASYwqs  84 

ZLGS  =  ^WGS  84  AZ  V^WGS  84  +  ®  ^WGS  84  ~  ASZwqs  84 


6.2.3  Parameter  values.  The  parameters  AX,  AY,  AZ,  CO,  £, 

\|f,  and  AS,  for  the  European  1950  and  Ordinance  Survey  of  Great 
Britain  1936  datums  are  tabulated  in  Appendix  E.  Use  of  seven- 
parameter  datum  shifts  for  some  applications  in  Europe  and  Great 
Britain  is  proscribed  by  STANAG  2211.  Since  the  values  of  the 
rotational  parameters  CO,  £,  and  \|/  are  given  in  the  table  in  arc 
seconds,  these  values  must  be  converted  to  radians  before  using 
equation  sets  6 . 1  and  6.2. 

6.2.4  Accuracy  of  seven-parameter  transformation.  The  RMS 
differences  between  WGS  84  surveyed  and  transformed  coordinates 
using  the  seven-parameter  datum  transformations  are  given  in 
Appendix  E.  For  some  datums,  the  accuracies  for  the  seven- 
parameter  transformation  are  not  an  improvement  over  those  for  the 
three-parameter  transformation.  Equation  sets  6.1  and  6.2  provide 
good  approximations  to  a  rigorous  orthogonal  transformation  when 

CD,  £,  \|/,  and  As  are  small. 

6.3  Three-parameter (AX,  AY.  AZ)  geometric  transformation. 
Frequently,  it  is  not  practical  to  determine  the  parameters  for  a 
seven -parameter  geometric  datum  shift  or  implement  such  a  shift 
method.  Often,  only  an  origin  shift,  three-parameter  model  is 
used  as  shown  in  figure  6.2.  The  three-parameter  datum 
transformation  satisfies  the  requirements  of  most  mapping  and 
charting  applications.  The  three -parameter  model  can  be  applied 
directly  to  Cartesian  coordinates  (see  Sections  6.3.1  and  6.3.2) 
or  used  to  compute  shifts  in  geodetic  coordinates  via  the 
Molodensky  equations  (see  Section  6.3.3) 

a.  The  datum  shift  parameters  (AX,  AY,  AZ)  are  the 
coordinates  of  the  origin  of  the  local  reference 
ellipsoid  in  the  WGS  84  Cartesian  coordinate  system. 

b.  Datum  shift  parameters  for  113  local  geodetic  systems  are 
listed  in  Appendix  B,  Table  B.l.  The  area  covered  by  the 
datum  shift  parameters,  the  local  datum  ellipsoid,  and 
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the  estimates  of  the  errors  in  the  datum  shift  parameters 
also  are  given.  Derivation  of  AX,  Ay,  and  Az  datum  shift 
parameters  is  discussed  in  DMA  TR  8350. 2-A  (1987). 


FIGURE  6.2.  Relationship  between  coordinate  axes  in  a 
three -parameter  geometric  transformation. 

6.3.1  Transformation  to  WGS  84  Cartesian  coordinates.  In 
certain  cases,  AX,  AY,  and  AZ  datum  shift  parameters  can  be 
applied  directly  to  convert  local  geodetic  system  coordinates  to 
WGS  84  Cartesian  coordinates  as  follows: 


^WGS 

84 

=  ^lgs 

+ 

AX 

(6.3) 

^WGS 

84 

=  ^LGS 

+ 

AY 

^WGS 

84 

=  zLGS 

+ 

AZ 

Equation  set  6.3  is  a  subset  of  equation  set  6.1  with  (0,  \|/,  8,  and 
AS  set  to  zero. 

6. 3. 1.1  Three-steo  method:  Transformation  to  WGS  84  geodetic 
coordinates .  Local  geodetic  coordinates  can  be  shifted  to  WGS  84 
coordinates  in  three  steps: 

a.  Convert  local  geodetic  coordinates  to  local  geocentric 
Cartesian  coordinates  using  equation  set  4.4.  Ensure 
that  local  ellipsoid  parameters  are  used. 

b.  Shift  the  local  geocentric  Cartesian  coordinates  to 
WGS  84  geocentric  Cartesian  coordinates  using  equation 
set  6.3. 
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c.  Convert  WGS  84  geocentric  Cartesian  coordinates  to  WGS 
84  geodetic  coordinates  using  the  methods  discussed  in 
Section  4. 3. 4. 2.  Ensure  that  WGS  1984  ellipsoid 
parameters  are  used. 

6.3.2  Transformation  to  local  geodetic  coordinates .  To 
convert  WGS  84  Cartesian  coordinates  to  the  local  geodetic  system, 
apply  AX,  AY,  and  AZ  datum  shifts  as  follows 


LGS  “ 

^WGS 

84 

-  AX 

(6.4) 

LGS  “ 

"^WGS 

84 

-  AY 

'LGS  ~ 

^WGS 

84 

-  AZ 

Equation  set  6.4  is  a  subset  of  equation  set  6.2  with  (0,  \|f,  8,  and 
AS  set  to  zero. 

6. 3. 2.1  Three-steo  method:  Transformation  to  local  geodetic 
coordinates .  To  shift  WGS  84  geodetic  coordinates  to  local 
geodetic  coordinates 

a.  Convert  WGS  84  geodetic  coordinates  to  WGS  84  geocentric 
Cartesian  coordinates  using  equation  set  4.4.  Insure 
that  WGS  84  ellipsoid  parameters  are  used. 

b.  Shift  the  WGS  84  geocentric  Cartesian  coordinates  to 
local  geocentric  Cartesian  coordinates  using  equation 
set  6.4. 

c.  Convert  the  local  geocentric  Cartesian  coordinates  to 
local  geodetic  coordinates  using  a  method  of  Section 
4 . 3 . 4 . 2 .  Insure  that  local  ellipsoid  parameters  are 
used. 

6.3.3  Molodenskv  shifts.  The  standard  Molodensky  method  is 
an  approximation  to  the  three-step  transformation  methods  of 
Sections  6. 3. 1.1  and  6. 3. 2.1.  To  use  the  Molodensky  methods  to 
transform  local  geodetic  coordinates  to  WGS  84  geodetic 

coordinates,  calculate  A<|),  AA,,  and  Ah  shifts  using  the  standard 
Molodensky  formulas  Rapp,  1987. 


A(j)  =  [-  sin  <])  cos  A,  AX  -  sin  0  sin  A,  AY  +  cos  <])  AZ 

E2  sin  6  cos  d)  . 

+ - fr - x  Aa 

W 

+  sin(|)cos<]) (2N  +  E'2  M  sin2  (j))  (l-f)Af]/[M  +  h] 
[-sin  A,  AX  +  cos  A,  AY] 


(6.5) 


AA,  = 


[(N  +  h)  cos  <|>] 
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Ah  =  cos  0  cos  A,  AX  +  cos  <})  sin  A  Ay  +  sin  (|)  AZ 
-  W  Aa  +  a(1w~  £)  sin2  (j)  (Af ) 


in  which 


W2 
M  = 


=  1  -  £2  sin2  <]) 

a(l  -  £2) 

W3 


4>  =  <t>LGS 


=  2f  -  f2 


a2  -  b2  _  e2 
b2  1  -  e2 


The  units  of  A(|)  and  AA  are  radians.  Note  that  N  =  RN  of  equation 
set  4.5. 


Note:  Most  applications  do  not  require  shifting  geodetic  heights, 
so  Ah  does  not  need  to  be  computed  or  applied.  Geodetic  heights, 
which  are  used  to  compute  datum  shift  parameters,  are  a  special 
construct  that  may  not  be  well  related  to  the  local  geodetic 
system.  Do  not  use  Ah  to  shift  local  geodetic  heights.  If  this 
is  a  requirement,  seek  competent  geodetic  council. 

a.  Apply  the  A<|>,  AA,  and  Ah  shifts  as  follows: 


^WGS  84 
^WGS  84 
fr-WGS  84 


=  ^LGS 
=  ^LGS 

=  b-LQg 


+  A^lgs 
+  A^lgs 

+  AhLGS 


(6.6) 


where  A(])LGS,  AA^g,  AhLGS  are  computed  from  equation  set 
6.5  using  <))LGS,  ALGS,  hLGS ,  semi -major  axis  and  flattening 
for  the  local  geodetic  system.  AX,  AY,  and  AZ  are  found 
in  Appendix  B.  Aa  and  Af  are  calculated  as  follows; 
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Aa  -  aWGS  84  a^g 

Af  ~  ^WGS  84  ”  ^liGS 

Note:  See  section  6.5  for  use  of  hLGS 

b.  The  Molodensky  formulas  can  also  be  used  to  convert  from 
WGS  84  to  the  local  geodetic  system: 


0LGS  =  ^WGS  84  +  ^WGS  84  (6.7) 

^LGS  =  ^VJGS  84  +  ^^WGS  84 
^LGS  =  hwGS  84  +  AhWGS  84 

where  A(|)WGS  84,  AA^g  84,  and  AhWGS  84  are  computed  using  (J)WGS  84, 
^wgs  84'  hwGs  84'  WGS  84  semi-major  axis  and  flattening,  and 
the  signs  of  Ax,  Ay,  and  Az,  as  given  in  Appendix  B,  are 
reversed.  Aa  and  Af  are  calculated  as  follows 

Aa  =  aLGS  -  aWGS  84 

Af  =  f lGS  ~  ^WGS  84 

Note  that  ellipsoid  parameters  a  and  f  are  found  in  Table  A.l, 
Appendix  A. 

6.3.4  AX,  AY,  and  AZ  shift  values.  Appendix  B  (Table  B.l) 

contains  the  AX,  AY,  and  AZ  values  for  the  three-parameter  methods 
of  converting  local  geodetic  systems  to  WGS  84.  The  accuracies, 
in  terms  of  error  estimates  for  AX,  AY,  and  AZ,  are  given  in 
Appendix  B.  A  further  discussion  of  the  accuracy  of  these 
parameter  shifts  will  be  given  in  Section  6.3.5. 

6.3.5  Accuracy  of  AX,  AY,  and  AZ  shift  parameters.  Error 

estimates  for  Ax,  Ay,  and  Az  are  tabulated  in  Appendix  B.  These 
estimates  include  the  errors  associated  with  the  Doppler  stations 
as  well  as  the  residual  differences  between  transformed  local 
system  coordinates  and  the  reference  WGS  84  coordinates  used  to 
develop  the  datum  shifts.  (A  Doppler  station  refers  to  a  position 
where  the  WGS  84  coordinates  were  obtained  by  observing  and 
processing  of  TRANSIT  satellite  data.) 

a.  For  38  datums,  datum  shifts  were  developed  using  a  single 
Doppler  station.  For  these  datums,  no  measure  of 
internal  consistency  of  the  datum  shifts  is  calculable, 
and  DMA  has  assigned  an  accuracy  value  of  25  meters  in 
each  coordinate.  These  error  statistics  do  not  reflect 
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any  errors  in  the  coordinates  used  to  compute  the  datum 
shift . 

b.  Datum  shifts  for  eight  datums  were  calculated  from 
sources  other  than  Doppler  stations  for  WGS  84  geodetic 
coordinates.  No  error  statistics  are  provided  for  these 
datums . 

c.  Geodetic  coordinate  difference  error  estimates  can  be 
approximated  from  rectangular  coordinate  errors  using  the 
following  spherical  relationships: 

0A(j,  =  •>/  (CT^x  cosA,)2  +  (GAy  sin(|)  sinA,)2  +  (CAz  cos  (|))2 

=  “\/  (  Gax  sinA,)2  +  (GAy  cosA,)2  (6.8) 

=  "V  (0>ax  coslj)  cosA,)2  +  (GAy  cos(|)  sinA,)2  +  (GAz  sin(|)) 2 

Note  that  these  relationships  neglect  any  correlations 
that  may  exist  between  the  error  estimates  for  Ax,  Ay 
and  Az. 

6.3.6  Local  datum  distortion.  For  some  local  datums,  such 
as  the  Australian  Geodetic  datums,  one  set  of  Ax,  AY,  and  AZ  datum 
shift  parameters  provides  a  consistent  measure  of  the  datum  shifts 
throughout  the  datum.  For  other  datums.  Ax,  Ay,  and  Az  shift 
parameters  for  the  whole  local  datum,  called  mean-value  datum 
shift  parameters,  poorly  fit  one  or  more  regions  within  the  datum. 
Regional  variation  in  the  datum  shift  parameters  may  result  from 
using  different,  poorly  connected  surveys  to  establish  the 
geodetic  control  on  the  local  datum;  from  rotation  of  the  local 
datum;  and  from  scale  differences.  An  extreme  example  of  regional 
datum  shift  variation  is  the  AY  component  for  the  ARC  1950  datum. 
The  mean-value  AY  is  -90  meters.  However,  Ay  ranges  from  -5  to  - 
108  meters  for  different  countries  covered  by  the  ARC  1950  datum. 
Appendix  B  (Table  B.l)  contains  regional  Ax,  Ay,  and  AZ  datum 
shift  parameters  for  17  local  datums.  For  these  datums,  one  should 
base  the  selection  of  mean  or  regional  AX,  AY,  and  AZ  shift 
values  on  the  accuracy  required  for  the  application.  Only  mean 
value  three-parameter  shifts  are  available  for  the  remaining 
datums  listed  in  Table  B.l. 

6.4  WGS  72  to  WGS  84  transformation.  WGS  84  has  replaced 
WGS  72  as  the  accepted  geodetic  reference  system  for  most  DoD 
applications.  The  WGS  84  system,  developed  through  a  more 
extensive  set  of  satellite-derived  and  surface  data  than  was 
available  at  the  start  of  WGS  72,  is  an  improved  geometric  and 
gravitational  model  of  the  earth.  Whenever  possible,  WGS  72 
coordinates  should  be  converted  to  the  WGS  84  reference  system. 
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6.4.1  Selecting  a  conversion  method.  Before  selecting  a  WGS 
72  to  WGS  84  conversion  method,  find  the  source  of  the  WGS  72 
coordinates.  Select  the  conversion  method  using  the  guidance  in 
Table  6.1. 


TABLE  6.1.  Selecting  a  WGS  72  to  WGS  84  Transformation  Method. 


Source  of  WGS  72  Coordinates 

1.  Doppler  satellite  station. 

2.  Local  datum  to  WGS  72  transformation. 
Local  datum  coordinates  known. 


3.  Local  datum  to  WGS  72  transformation. 

Local  datum  coordinates  unknown. 

4.  Local  datum  to  WGS  72  transformation 
using  localized  datum  shift  parameters 
derived  for  nearby  WGS  72  Doppler 
station.  Local  datum  coordinates  unknown. 

5.  Unknown.  Direct  (6.4.2) 

Accuracy  unknown 

6.4.2  Direct  WGS  72  to  WGS  84  transformation.  The  direct 
WGS  72  to  WGS  84  transformation  reflects  fundamental  changes 
between  the  WGS  72  and  WGS  84  systems.  Changes  include  a  shift  in 
coordinate  system  origin,  a  shift  in  longitudinal  reference 
(initial  meridian) ,  a  change  in  system  scale,  and  changes  in 
ellipsoidal  parameters.  For  WGS  72  coordinates  not  originating  in 
the  WGS  72  system  (i.e.,  non-Doppler-derived  coordinates),  take 
care  when  employing  the  direct  WGS  72  to  WGS  84  transformation,  as 
any  inaccuracies  and  uncertainties  inherent  in  the  WGS  72 
coordinates  will  be  directly  transferred  into  the  derived  WGS  84 
coordinates.  These  inaccuracies  and  uncertainties  are  generally  a 
result  of  local  geodetic  system  to  WGS  72  coordinate 
transformations  using  mean  rather  than  localized  datum  shifts. 
Thus,  every  effort  should  be  made  to  find  the  source  of  all  WGS  72 
coordinates  before  doing  a  direct  WGS  72  to  WGS  84  transformation. 


Transformation  Method 

Direct  (6.4.2) 

Local  datum  to  WGS  84 
transformation 
(6.2,  6.3,  6.6) 

Two-Step  (6.4.3) 

Shift  method  known 

Direct  (6.4.2) 


6. 4. 2.1  Conversion  equations.  It  is  important  to  note  that 
the  relationship  between  the  WGS  72  coordinates  and  the  WGS  84 
coordinates  cannot  be  described  by  a  three -parameter  model  (see 
Section  6.3) . 
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a.  The  equations  for  the  direct  WGS  72  to  WGS  84  geodetic 
coordinate  transformation  are 


rWGS  84  -  t WGS  72 


+  A(j) 


(6.9) 


^WGS  84  ~  ^WGS  72  +  AA 
h-WGS  84  =  -^WGS  72  +  Ah 


where 


4.5  cos<|)WGS  72  Af  sin  2<J)WGS  72 

A(p  =  - - ~z —  +  - r -  (arc  seconds) 

awcs  72  v  W 


AA  =  0.554 


(arc  seconds) 


Ah  -  4.5  sin([)WGS  72  +  9-wgs  72  Af  sin2 


Kwgs  72 


-  Aa  +  Ar  (meters) 


Af  =  0.3121057  x  10-7 


awGs  72  =  6378135  (meters) 


Aa  =  2.0 


(meters) 


Ar  =  1.4 


(meters) 


180  •  3600 

b.  The  equations  for  the  direct  WGS  84  to  WGS  72  geodetic 
coordinate  transformation  are: 


^WGS  72  =  <j>WGS  84  +  A(|) 

A(JGS  72  =  AwgS  84  +  AA 

h-WGS  72  =  h.WGS  84  +  Ah 


(6.10) 


where 


A(j)  = 


-4.5  cosOmgs  84  Af  sin  2(])WGg  84 

aWGS  84  Q  Q 


(arc  seconds) 


AA  =  -0.554 


(arc  seconds) 


Ah  =  -4.5  sin(|)  -  a, 


WGS  84 


Af  sin2  § 


WGS  84 


+  Aa  -  Ar  (meters) 


Af  =  0.3121057  x  10-7 
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6378137  (meters) 


aWGS  84  _ 

Aa  =  2.0  (meters) 

Ar  =  1.4  (meters) 


180  •  3600 

Elevation  based  on  the  WGS  84  geoid  model  will  be  different  from 
elevations  based  on  the  WGS  72  geoid  model  due  to  differences  in 
the  geoid  models  (see  Section  6.7.1  for  conversion). 

6.. 4.2.2  Estimated  errors  (G^.O^.G^)  .  The  estimated  errors 
in  a  direct  WGS  72  to  WGS  84  transformation  are 

°A<t>  =  °AX  =  3  m  (6.11) 

=  3  m  to  4  m 

These  error  estimates  do  not  include  errors  in  the  original  WGS  72 
coordinates.  The  accuracy  estimate  for  G^  reflects  discrepancies 
between  WGS  72  and  WGS  84  geoid  models. 

6.4.3  Two-step  WGS  72  to  WGS  84  transformation.  If  WGS 
72  coordinates  are  known  to  have  been  derived  from  local  geodetic 
system  coordinates  using  WGS  72  Ax,  AY,  and  AZ  datum  shift 
parameters,  the  preferred  method  for  transforming  them  to  the  WGS 
84  system  is  to  do  a  two-step  conversion.  First,  do  a  WGS  72  to 
local  transformation  using  the  inverse  of  the  transformation 
originally  used  to  shift  local  datum  coordinates  to  WGS  72.  This 
step  will  eliminate  the  inaccuracies  introduced  using  local  to  WGS 
72  datum  shifts.  Next,  do  a  local  geodetic  system  to  WGS  84 
conversion  to  produce  coordinates  within  the  WGS  84  system  (see 
Sections  6.2,  6.3,  or  6.6).  Seppelin  (1974)  provides  mean  value 
local  datum  to  WGS  72  datum  shift  parameters. 

6.5  Approximating  geodetic  heights  for  datum 
Transformation .  When  doing  the  Molodensky  approximation  (Section 
6.3.3),  ellipsoidal  height  (h)  is  often  unavailable.  MSL  height 
(H)  can  be  substituted  for  h  without  introducing  errors  in 
horizontal  position  significant  to  mapping  and  charting 
applications.  The  inpact  of  this  substitution  is  so  small  that 
zero  can  be  used  in  place  of  h,  if  H  is  unavailable. 

NOTE:  The  generation  and  usage  of  local  datum  ellipsoidal 

heights  are  outside  the  scope  of  this  Technical  Report.  For  such 
applications,  consult  TEC  at  the  address  in  the  foreward. 
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6.6  Multiple  regression  equations  (MRE) .  The  use  of  MREs 
for  datum  shifting  is  not  recommended. 

6.7  Vertical  datum  shifts.  Because  exact  relationships 
between  vertical  datums  are  unknown,  shifting  elevations  between 
different  vertical  datums  should  not  be  attempted,  except  for  WGS 
72  elevations  conputed  from  measured  WGS  72  geodetic  heights  (see 
Section  6.7.1.)  Otherwise,  the  following  are  assumed  to  be  the 
same: 


a.  Elevations  on  different  local  vertical  datums. 

b.  Elevations  computed  from  measured  WGS  84  geodetic  heights 
and  WGS  84  geoid  separations . 

It  is  estimated  that  errors  resulting  from  this  assumption  will  be 
one  to  two  meters.  Should  the  accuracy  be  insufficient,  contact 
TEC  at  the  address  listed  in  the  foreword  for  assistance. 

6.7.1  Satellite-derived  WGS  72  elevations  based  on  measured 
WGS  72  geodetic  heights.  The  WGS  84  geoid  model  is  a  closer 
approximation  to  the  geoid  than  the  older  WGS  72  geoid  model. 
Elevations  obtained  from  WGS  72  geodetic  height  measurements  and 
geoid  separations  may  be  improved  by  converting  to  WGS  84  using 

H  =  Hwgs  72  +  Nwgs  72  +  All  -  Nwgs  34  ( 6 . 12 ) 

where  HWGS  72  is  the  elevation  computed  using  the  WGS  72  geoid, 

Nwgs  72  is  the  WGS  72  geoid  separation,  and  Ah  is  computed  using 
equation  set  (6.9). 

6.7.2  Survey-derived  elevations.  Elevations  obtained  by 
conventional  leveling  surveys,  for  points  with  WGS  72  geodetic 
horizontal  coordinates,  should  not  be  converted  to  WGS  84 
elevations.  Simplified  procedures  to  convert  WGS  72-derived 
elevations  to  WGS  84  are  not  recommended. 

6.8  Unlisted  datums.  This  Technical  Report  provides  datum 
shift  parameters  between  113  local  datums  and  WGS  84.  Regional 
shift  parameters  are  also  given  for  portions  of  many  of  these 
datums.  Contact  TEC  at  the  address  given  in  the  foreword  for  any 
datums  and  their  respective  shift  parameters  that  are  not  listed 
in  this  Technical  Report. 
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7. 


MAP  PROJECTIONS 


7.1  Introduction.  Map  projections  are  used  to  represent  the 
earth's  features  or  points  on  a  map  or  in  a  plane  coordinate 
system.  A  map  projection  is  a  systematic  representation  of  a  part 
of  the  earth's  surface  on  a  plane. 

7.1.1  The  mao  projection  process.  The  map  projection 
process  can  be  done  in  three  steps: 

(1)  The  positions  of  points  on  the  surface  of  the  earth  (<j),  X, 
h)  are  reduced  to  positions  on  the  reference  ellipsoid  ((|), 

X,  h  =  0) . 

(2)  Next,  a  map  projection  is  chosen  and  used  to  transform 
the  feature's  geodetic  coordinates  (<j),  X,  h  =  0)  into 
planar  map  projection  coordinates  (x,  y)  . 

(3)  Finally,  the  map  projection  coordinates  of  the  features 
may  be  plotted,  on  a  map,  or  stored  for  future  use. 

7.1.2  Properties  of  projections.  Once  the  earth's  surface 
is  projected  on  a  plane,  one  can  observe  such  properties  as 
distance,  angle,  direction,  shape,  and  size.  None  of  the 
projections  described  in  this  Technical  Report  faithfully 
preserves  all  these  characteristics.  Some  type  of  distortion  is 
always  introduced  by  a  map  projection.  However,  a  given  map 
projection  can  preserve  some  of  these  properties,  and  if  a 
sufficiently  small  part  of  the  earth's  surface  is  displayed,  then 
other  properties  can  be  represented  approximately.  Map 
projections  presented  in  this  Technical  Report  have  areas  of 
coverage  beyond  which  the  projection  should  not  be  used  due  to 
increased  distortion. 

7.1.3  Mapping  equations.  The  map  projections  discussed  in 
this  Technical  Report  have  mapping  equations  that  are  defined  on  a 
reference  ellipsoid  and  are  functions  of  geodetic  coordinates 

(<|),  X)  .  Geodetic  height,  h,  does  not  enter  into  the  projection 
equations  since  h  =  0  on  a  reference  ellipsoid.  The  equations 
that  express  projection  coordinates,  x  and  y,  in  terms  of  (0,  A,) 
are  sometimes  called  mapping  equations. 

7.1.4  Conformal  projection.  Each  projection  presented  in 
this  Technical  Report  has  the  important  property  of  being 
conformal.  A  conformal  projection  is  angle  preserving.  If  two 
lines  on  the  ellipsoid  meet  at  an  angle  of  0,  then  in  a  conformal 
projection  the  image  of  these  lines  on  the  map  meet  at  the  same 
angle ,  0 . 

7.1.5  Scale  factor.  Two  quantities  associated  with  map 
projections,  which  are  often  confused,  are  the  scale  factor  of  the 
projection  and  the  map  scale.  The  scale  factor  of  the  projection 
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is  used  in  precise  direction  and  distance  calculations  using  a  map 
projection's  coordinates.  The  scale  factor  of  a  projection  is  the 
ratio  of  arc  length  along  a  differentially  small  line  in  the  plane 
of  the  projection  to  arc  length  on  the  ellipsoid.  This  number 
depends  on  both  the  location  of  the  point  and  on  the  direction  of 
the  line  along  which  arc  length  is  being  measured.  However,  for 
conformal  projections,  the  scale  factor  of  the  projection  is 
independent  of  the  direction  of  the  line  and  depends  only  upon  the 
location  of  the  point.  Since  all  projections  treated  in  this 
Technical  Report  are  conformal,  the  scale  of  the  projection  will 
be  referred  to  as  the  point  scale  factor  (k) . 

7.1.6  Mao  scale.  The  map  scale  is  the  approximate  number 
for  converting  map  distances  to  terrestrial  distances  and  is 
usually  published  on  the  map.  Technically,  the  map  scale  is  the 
ratio  of  a  distance  on  the  map  to  the  corresponding  distance  on 
the  ellipsoid.  Although  the  map  scale  changes  over  the  area  of  a 
map  sheet,  it  is  customary  to  assign  a  single  value  to  denote  the 
map's  scale.  This  value  is  the  true  map  scale  along  at  least  one 
line  on  the  map.  Usually,  the  map  scale  is  given  by  a  ratio  l:n, 
so  that  the  larger  the  map  scale,  the  smaller  the  value  of  n.  For 
military  maps  of  scale  greater  than  1:100,000  (n  <  100,000),  the 
map  scale  changes  only  slightly  over  a  map  sheet,  and  the  map 
scale  approximates  the  ratio  of  distance  on  the  map  to  distance  on 
the  ellipsoid. 

7.1.7  Convergence  of  the  meridian.  The  convergence  of  the 
meridian  (y) ,  often  called  grid  convergence,  is  used  to  convert 
between  direction  in  the  map  projection  (grid  azimuth)  and 
geodetic  azimuth.  The  convergence  of  the  meridian  at  a  point  on  a 
map  is  the  angle  between  the  projection  of  the  meridian  at  that 
point  and  the  projection  North  axis  (typically  represented  by  the 
y  axis  of  the  projection) .  The  convergence  of  the  meridian  is 
positive  in  the  clockwise  direction.  Convergence  of  the  meridian 
should  not  be  confused  with  declination,  which  is  the  difference 
between  geodetic  and  magnetic  north  at  the  point. 

7.1.8  Information  sources.  Information  on  map  projections 
can  be  found  in  Krakiwsky  (1973);  NSWC/DL  TR-3624  (1977);  U.S. 
Department  of  Commerce,  Coast  and  Geodetic  Survey,  Special 
Publication  No.  251  (1979);  and  Geological  Survey  Professional 
Paper  1395  (1987). 

7.2  Mercator  oro-iection.  The  Mercator  projection  is  a. 
conformal  projection  for  which  the  point  scale  factor  is  one  along 
the  equator.  The  equator  lies  on  the  line  y  =  0.  This  projection 
is  not  defined  at  the  poles. 

7.2.1  Meridians  and  parallels.  Meridians  and  parallels 
provide  the  framework  for  the  Mercator  projection.  Meridians  are 
projected  as  parallel  straight  lines  that  satisfy  the  equation  x  = 
constant.  Evenly  spaced  meridians  on  the  ellipsoid  project  to 
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evenly  spaced  straight  lines  on  the  projection.  Parallels  are 
projected  as  parallel  straight  lines  perpendicular  to  meridians 
and  satisfy  the  equation  y  =  constant.  Evenly  spaced  parallels  on 
the  ellipsoid  project  to  unevenly  spaced  parallels  on  the 
projection.  The  spacing  between  projected  parallels  increases 
with  distance  from  the  equator.  See  figure  7.1. 


FIGURE  7.1.  Meridians  (longitude)  and  parallels  (latitude) 

in  the  Mercator  projection. 

7.2.2  Mercator  mapping  equations. 

7. 2. 2.1  Finding  (x.v) .  The  following  discussion  of  the 
Mercator  mapping  equations  can  be  found  in  Geological  Survey 
Professional  Paper  1395  (1987) .  First,  choose  a  central  meridian 
(A,0)  that  represents  the  zero  point  on  the  x  axis.  Given  a  point 
((j),  X) ,  in  radians,  on  the  ellipsoid,  the  corresponding  point 
(x,  y)  in  the  plane  is 

x  =  a(X,  -  X0)  (7.1) 


fn  0Y  l  -  e  sin  0  ^ 

E 

2 

=  a  In 

tan  .  +  - — : — r 

v4  2^  i  +  £  sin 

_ 

or 


1  +  sin  <J)Yl  -  £  sin  0  ^ 

i 

W 

_JL  -  sin  <j)J^l  +  E  sin  <])  j 
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where  a  =  semi -major  axis  of  the  ellipsoid 
£  =  eccentricity 

XQ  =  longitude  of  the  central  meridian 

<j),  X  and  Xq  are  expressed  in  radians. 

7. 2. 2. 2  Point  scale  factor  and  convergence  of  the  meridian. 
The  equations  for  the  point  scale  factor  (k)  and  the  convergence 
of  the  meridian  (y)  are 


N  cos  <|) 

y  =  0 


(7.3) 


where 


N  = 


'1  -  £2  sin2  (j) 


7. 2. 2. 3  Finding  ((|).A,)  .  Next,  consider  the  inverse  problem 
of  finding  (()>,  X)  when  (x,  y)  is  known.  The  longitude  can  be 
found  from 


x  =  -  ♦  jl, 


(7.4) 


where  both  X  and  XQ  are  expressed  in  radians.  The  value  of  (])  can 
be  found  iteratively.  Let  an  intermediate  parameter  t  be  defined 


t  =  e 


(7.5) 


After  an  initial  value  of 


4>o  = 


2  arctan  t 


(7.6) 


iterate  the  following  equation  for  (|): 


,  „  fj(l  -  £  sin  6)1  | 

^n+1  =  J  -  2  arctan  t  — - - — ; — r— 

Tn+i  2  M  (1  +  £  sxn  <pn)_ 


(7.7) 
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Substitute  the  initial  value  (<J)0)  into  equation  7.7  to  find  the 
next  candidate  Similarly,  substitute  <{>]_  into  equation  7.7  to 

obtain  the  updated  candidate,  Continue  this  process  until  the 

difference  between  successive  values  of  <j)n  is  sufficiently  small. 
Then  set  ♦  =  4>n.  All  angles  are  measured  in  radians. 

7.2.3  Accuracy .  The  Mercator  equations  for  x,  y,  k,  y,  and  X 
are  exact.  The  iterative  equation  for  <J)  can  be  updated  until  any 
desired  accuracy  is  obtained. 

7.2.4  Area  of  coverage.  In  the  Mercator  projection,  as  the 
latitude  (<|))  approaches  the  poles,  the  y  coordinate  approaches 
infinity.  Area  and  length  distortion  increases  with  distance  from 
the  equator.  For  example,  the  point  scale  factor  is  approximately 
2  at  60°  latitude  and  5.7  at  80°  latitude. 

7.3  Transverse  Mercator  (TM)  projection.  The  Transverse 
Mercator  projection  is  a  conformal  projection  for  which  the  point 
scale  factor  equals  one  along  the  central  meridian.  The  line 
y  =  0  is  the  projection  of  the  equator,  and  the  line  x  =  0  is  the 
projection  of  the  central  meridian.  The  Universal  Transverse 
Mercator  (UTM)  projection  (see  Section  7.4)  is  a  modification  of 
the  Transverse  Mercator  projection. 

7.3.1  Meridians  and  parallels.  Both  the  central  meridian 
and  the  equator  are  represented  as  straight  lines.  No  other 
meridian  or  parallel  is  projected  onto  a  straight  line.  Since  the 
point  scale  factor  is  one  along  the  central  meridian,  this 
projection  is  most  useful  near  the  central  meridian.  Scale 
distortion  increases  away  from  this  meridian.  (See  Figure  7.2.) 
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FIGURE  7.2.  Meridians  and  parallels  in  the  Transverse 
Mercator  Projection.  A.0  is  the  central 
meridian . 

7.3.2  Transverse  Mercator  mapping  equations. 

7. 3. 2.1  Finding  (x.v) .  The  presentation  here  follows  U.S. 
Department  of  Commerce,  Coast  and  Geodetic  Survey,  Special 
Publication  No.  251  (1979).  A  point  with  geodetic  coordinates  (<j), 
X)  in  radians  has  Transverse  Mercator  coordinates  (x,y)  that  are 
given  by 

x  =  N  A  cos(|)  +  —  —  g °S — ^  (1  -  t2  +  T|2)  (7.8! 

+  N  (5  -  18t2  +  t4  +  14T|2  -  58t2TJ2) 


N  A2  .  .  . 

y  =  S^  +  — 2 —  sin  r  cos  t 

+  sin  0  cos3  <j)  ( 5  -  t2  +  9T|2  +  4T|4) 

N  A6 

+  ~^r, r-  sin  <|)  cos5  Cb  (61  -  58t2  +  t4  +  270T|2  -  330t2Tj2) 


where 


N  = 

t  : 

11  = 
£•2 

£2  = 


Vl  -  £2  sin2  <j) 
:  tan  <j) 

E’cos  (J) 

£2 

(1  -  e2) 
a2  -  b2 
a2 


S<j,  =  a[A0<}>  -  A2  sin  2(|)  +  A4  sin  4<|) 

-  A6  sin  6<j)  +  A8  sin  8<|)] . 

— — — eg  -IZ5 

25o  16384 

15pe  _  455 

128 


A0  = 


A,  = 


-  ^£2  _  -2-£4  _ 

4*  64* 


£8 


f  (£2  +  * 


4096 


€8) 


A,  = 


A«  = 


*/P6  4 ^£8 \ 

3072  *  32*  ; 


315 


£8 


~8  “  131072 

A  =  X,  -  A.0  in  radians 
A,0  =  central  meridian  in  radians 


(7.9) 


At  the  poles,  where  <J>  =  +  or  -  1Z/2,  the  Transverse  Mercator 
coordinates  are  x  =  0  and  y  =  S^,  a  quadrant  arc  on  the  ellipsoid. 

7. 3. 2. 2  Finding  «)).  X)  . 


7. 3. 2. 2.1  Finding  the  footooint  latitude  (<j)j)  .  To  find  ((j), 
X),  given  x  and  y,  first  compute  the  footpoint  latitude  using 

the  following  iterative  method  (See  Krakiwsky  1973).  The  equation 
to  be  iterated  to  yield  the  footpoint  latitude  ((j^)  is 

0'  =  4>  -  SV~— y  (7.10) 


where 


S'<(,  =  a[A0  -  2A2cos2(|)  +  4A4cos40  -  6A6cos6<|)  +  8A8cos8(|)]  (7.11) 

and  Aq,  A2,  A4,  Ag,  and  Ag  are  defined  in  equation  set  7.9. 

Let  the  first  candidate  for  <)>!  be  <|)  =  y/a.  When  this  value  of 
is  substituted  into  equation  7.10,  the  updated  candidate  for  is 
(jv.  Similarly,  the  substitution  of  (j)'  for  (|)  in  equation  7.10 
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gives  the  next  updated  value  of  <J)'.  When  successive  values  of  (j)' 
are  sufficiently  close,  then  the  footpoint  latitude  is  found  by 
setting  4>i  =  <|> *  ♦  A  non-iterative  method  for  finding  (j^  can  be 
found  in  NOAA  Report  NOS  NGS  5  (1989). 

7. 3. 2. 2. 2  Finding  (<[),A)  from  (j^.  Once  is  determined,  ((|), 
A,)  can  be  found  from  the  approximations 


X\4  B6  /X  NS!  - 

Nj  +  720  [nJJ 

(7.12) 

X°  *  SeC  ♦!  [  ^  ^  (j 

<  \3  B5  /x  \5  ' 

h)  + 120  w _ 

where 

B3  =  1  +  2ti  +  T\l  (7.13) 

B4  =  5  +  3ti  +  T\l  -  4T)i 
B5  =  5  +  28ti  +  24t}  +  6T|i  + 


B6  =  61  +  90ti  +  4611 1  +  45t}  -  252tiHj 
-  3J\l  -  66tJllJ  -  90tjT|i  +  225tJ  TjJ 


If  =  ±  then  <})  =  ±  90°  (a  pole)  and  A  is  undefined. 
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7. 3. 2. 3  Point  scale  factor  as  a  function  of  x  and  fo.  As  a 
function  of  x  and  (j^,  the  point  scale  factor  is  given  by 


*(fj 


where 


C2  =  1  +  %  (7 

C,  =  1  +  6T(i  +  9T|*  -  24t l  Tl[ 
and  tx  ,  TJ-l ,  and  Nx  are  given  by  equation  set  7.13. 

7.3. 2 .4  Point  scale  factor  as  a  function  of  ((j)A)  .  As  a 
function  of  (()),  X)  ,  the  point  scale  factor  is  given  by 


k  =  1  + 
in  which 


A2  cos2 


D2  + 


L4  COS4 


d4  + 


k.6  cos6 


(7.14) 


(7.15) 


(7.16) 


D2  =  1  +  Tl2 

D4  =  5  -  4t2  +  14T|2  +  13T|4  -  28t2T|2  -  48t2Tl4 
D6  =  61  -  148t2  +  16t4 

and  A,  t,  and  T]  are  given  in  equation  set  7.9. 


(7.17) 


7. 3. 2. 5  Convergence  of  the  meridian  in  terms  of  x  an 
The  convergence  of  the  meridian,  in  terms  of  x  and  (j^,  is 


t^x 

=  Ni  . 


i  &  (x. v  + 

3  In,  I  is  In,  ) 


/X  Y  E6  / x_Y  1 
W  "  315  [nJ  J 


(7.18) 


where 


E2  =  1  +  1 1  -T\l-  2T\i 

E4  =  2  +  5ti  +  21 +  3tJ  +  tit|i  +  9TlJ  -  7^} 

E6  =  17  +  77tJ  +  105tJ  +  45t* 
and  tlf  Tl!,  and  Nj  are  given  by  equation  set  7-13. 


(7.19) 


7. 3. 2. 6  Convergence  of  the  meridian  in  terms 
convergence  of  the  meridian,  in  terms  of  ( (j) .  AA  is 


i .  The 
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y  =  A 


sm 


1  +  F, 


A2 


cos^ 
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+  Fa 


A4 


cos- 
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15 


+  Ff 


Ae 


cos' 


1 


315 


(7.20) 


in  which 


F2  =  1  +  3TJ2  +  2T|4  (7.21) 

F4  =  2  -  t2  +  15TJ2  +  35T14  -  15t2TJ2  -  50t2114 
F6  =  17  -  26t2  +  2t4 


and  A,  t,  and  T)  are  given  in  equation  set  7.9. 

7.3.3  Accuracy .  The  Transverse  Mercator  equations  for  x, 
y,  <j),  X,  S^,  <|)1#  k,  and  y  are  approximations.  Within  4°  of  the 
central  meridian,  the  equations  for  x,  y,  <|),  and  X  have  an  error 
of  less  than  1  centimeter.  (See  Geological  Survey  Professional 
Paper  1395,  1987.) 

7.3.4  Area  of  coverage.  The  equations  given  here  for  the 
Transverse  Mercator  projection  can  be  used  within  4°  of  the  central 
meridian.  If  these  equations  are  used  with  UTM  coordinates,  then 
consult  Section  7.4.1  for  guidance  on  the  area  of  coverage  for  UTM 
coordinates . 


7.4  Universal  Transverse  Mercator  (UTM)  project ion.  The 
UTM  projection  is  a  family  of  projections  that  differ  from  the 
Transverse  Mercator  projection  in  these  ways  (see  Figure  7.3): 

a.  The  longitudes  of  the  central  meridians  (in  degrees)  are 

3  +  6n,  n  =  0, . . .59. 

b.  The  point  scale  factor  along  the  central  meridian  is 
0.9996. 

c.  The  y  value,  called  the  northing,  has  an  origin  of  0 
meters  in  the  northern  hemisphere  at  the  equator.  The  y 
value,  called  the  southing,  has  an  origin  of  10,000,000 
meters  at  the  equator  and  decreases  toward  the  pole  in 
the  southern  hemisphere. 

d.  The  x  value,  called  the  easting,  has  an  origin  of  500,000 
meters  at  the  central  meridian. 

7.4.1  UTM  zones .  For  the  UTM  grid  system,  the  ellipsoid  is 
divided  into  60  longitudinal  zones  of  6  degrees  each.  Zone  .number 
one  lies  between  180°  E  and  186°  E.  The  zone  numbers  increase 
consecutively  in  the  eastward  direction.  It  is  intended  that  a 
UTM  projection  should  be  used  only  in  one  of  the  6 -degree  zones, 
plus  the  overlap  area.  In  each  zone,  the  UTM  projection  extends 
for  a  40  km  overlap  into  the  two  adjacent  zones.  Moreover,  UTM 
coordinates  are  used  only  between  80°  S  and  84°  N  plus  the  overlap 
region.  The  UTM  projection  extends  to  84°  30 1  N  and  80°  30 1  S  to 
overlap  with  the  Universal  Polar  Stereographic  (UPS)  projections. 
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FIGURE  7.3.  Meridians  and  parallels  (dashed)  and  a 
Universal  Transverse  Mercator  Grid. 


7 . 4 . 1 . 1  Finding  the  UTM  zone  and  central  meridian .  To 
find  UTM  coordinates  from  geodetic  coordinates,  first  compute  the 
UTM  zone.  A  standard  zone  can  be  computed  directly  from  geodetic 
longitude.  If  the  geodetic  longitude  (A)  is  expressed  in  positive 
radians,  then  the  zone  (z)  is 


z 


z 


=  Greatest  Integer  < 
=  Greatest  Integer  < 


^31  + 

fl 80A 

i  6K 


for  0  <  A  <  71 


for  K  <  A  <  2n 


(7.22) 


Once  the  zone  (z)  is  known,  the  central  meridian  is  given  by 
equation  set  7.23.  If  a  longitude  lies  on  a  UTM  zone  boundary, 
then  the  user  must  choose  the  desired  UTM  zone.  If  a  map  is  being 
used,  then  the  UTM  zone  appears  on  the  map.  The  central  meridian 
of  zone  (z)  is 

A0  =  (6z  -  183)^  for  (z  >  31)  (7.23) 

A0  =  (6z  +  177)^q  for  (z  2  30) 

with  A0  expressed  in  radians. 
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7. 4. 1.2  Non-standard  width  UTM  zones.  Non-standard 
width  UTM  zones  are  used  in  North  Atlantic  regions  and  Norway  as 
delineated  in  Table  7.1  (see  DMA  TM  8358.1  (1990) .  These  non¬ 
standard  zones  may  also  be  applied  in  the  30'  UPS  overlap  region. 

TABLE  7.1.  Non-standard  zone  limits. 


Latitude 

Longitude 

Zone 

Lower 

West 

1  East 

Central  Meridian 

31 

s 

O 

in 

64°  N 

0°  E 

3°  E 

3°  E 

32 

56 

64 

3 

12 

9. 

31 

72 

84 

0 

9 

3 

32 

72 

84 

not  used 

33 

72 

84 

9 

21 

15 

34 

72 

84 

not  used 

35 

72 

84 

21 

33 

27 

36 

72 

84 

not  used 

37 

72 

84 

33 

42 

39 

7.4.2  Reference  source.  The  standard  references  for  the  UTM 
projection  are  DMA  TM  8358.1  (1990)  and  DMA  TM  8358.2  (1989). 

7.4.3  UTM  equations. 

7. 4. 3.1  Finding  UTM  coordinates  (xlw  vVTM)  .  Having  found 
the  central  meridian,  next  compute  Transverse  Mercator  coordinates 
(x,  y) ,  using  equation  sets  7.8  and  7.9.  Then,  the  UTM  coordinates 
can  be  found  from: 

(7.24) 

(Northern  hemisphere) 

(Southern  hemisphere) 


7. 4. 3. 2  Computing  geodetic  coordinates.  To  compute  geodetic 
coordinates  from  UTM  coordinates,  a  zone,  and  a  hemisphere,  first 
solve  equation  set  7.25  for  (x^,  yTM) ,  equation  set  7.23  for  XQi 
and  then  compute  (0,  X)  from  equation  set  7.12. 

■  (7.25) 

(Northern  hemisphere) 
(Southern  hemisphere) 


XTM 

Ytm 

Ytm 

kTM 

Ytm 


(  ^utm  —  500,000) 


0.9996 

0.9996  (yUTM> 

Q  gggg  (Yutm  -  10/000/000) 
- 1  k 
0.9996  u™ 

Yutm 


■'MJTM 

Yutm 

Yutm 

k0TM 

Yutm 


0.9996  XrM  +  500,000 
0.9996  yM 

0.9996  yTM  +  10,000,000 
0.9996 

Ytm 


TM 
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7. 4. 3. 3  Discontinuity .  Notice  that  UTM  coordinates  are 
discontinuous  in  y  at  the  equator. 

7.4.4  Accuracy .  The  accuracy  of  the  UTM  equations  is  based 
on  the  accuracy  of  the  Transverse  Mercator  equations.  Within  the 
area  of  coverage,  the  equations  for  x,  y,  <J),  and  X  have  an  error 
of  less  than  one  centimeter.  (See  Geological  Survey  Professional 
Paper  1395,  1987.) 


7.4.5  Area  of  coverage.  The  area  of  coverage  for  UTM 
coordinates  is  defined  by  zone  limits,  latitude  limits,  and 
overlap . 


Zone  limits:  6°  zones,  extending  3°  to  each  side 

of  the  central  meridian.  See  Table  7.1 
for  exceptions  to  this  general  rule. 

Latitude  limits:  North:  84° 

South:  80° 


Zone  overlap:  40  km  on  either  side  of  the  zone 

limits. 

Polar  overlap:  30’  toward  the  poles 

North:  84°  30' 

South:  80°  30' 


7.5  Lambert  Conformal  Conic  projection.  The  Lambert 
conformal  conic  projection  is  a  conformal  projection  in  which  the 
projected  parallels  are  arcs  of  concentric  circles  centered  at  the 
pole.  The  projected  meridians  are  radii  of  concentric  circles 
that  meet  at  the  pole.  Moreover,  there  are  one  or  two  parallels, 
called  standard  parallels,  along  which  the  point  scale  factor  is 
one.  (See  Figure  7.4.) 
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90°W  Pole  90°E 


FIGURE  7.4.  Meridians  and  parallels  in  the  Lambert  Conformal 
Conic  Projection.  (Not  drawn  to  scale.) 

7.5.1  Lambert.  Conformal  Conic  with  two  standard  parallels. 
This  presentation  of  the  mapping  equations  follows  Krakiwsky 
(1973).  The  first  case  to  be  considered  is  the  case  of  two 
standard  parallels,  <]>!  and  <j>2.  The  point  scale  factor  is  one 
along  both  standard  parallels.  Choose  a  fixed  central  meridian, 

X0 . 


7. 5. 1.1  Finding  (x.v)  .  Begin  by  fixing  a  latitude  (<j)0) 

that  is  below  the  area  of  interest.  For  a  given  (<(),  X)  the 
mapping  equations  are 


x  =  r  sin  LA 
y  =  rc  -  r  cos  LA 


where 


A  =  X  -  Xn 


r  =  K  e 


-Lq 


rG  =  K  e 


K  = 


L<Jo 
N2  COS 


N2  COS  (|); 


Le 


-Lq, 


Le 


-Lq. 


e  =  base  of  natural  logarithms  (  In  e  =  1  ) 
(()1  =  first  standard  parallel 
(j)2  =  second  standard  parallel 
a _ 

-  £2  sin2  0 


(7.26) 


(7.27) 
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q  =  In 


£2  = 


(TZ  0Yl  -  £  sin  0 
tanU+  2{!  +  e  sin  <|) 

a2  -  b2 
a2 


IndSfj  cos  (j^)  -  ln(N2  cos  <j)2) 

-  Qi 

qj  =  q  evaluated  at  ,  j  =  0,1,2 
Nj  =  N  evaluated  at  <j)j  ,  j  =  1,2 

7. 5. 1.2  Finding  {jjhAl  with  two  standard  parallels.  For  the 
reverse  case,  when  (x,  y)  is  given  along  with  <j)0,  tylt  <J)2 ,  and  X0i 
first  compute  L,  r0,  K,  Nx,  N2,  q^  qx,  and  q2  from  equation  set 
7.27.  Then  compute  X  and  q  from 


in  which 


r 

6 


y 

cos  0 


arctan 


x 

y 


(7.28) 


(7.29) 


Note:  When  evaluating  the  arc  tangent,  the  user  must  take 
care  to  insure  that  the  resulting  angle  lies  in  the  desired 
quadrant . 

Once  q  is  found,  then  <|>  can  be  determined  using  the  Newton  - 
Raphson  method,  which  is  explained  next.  Define  functions  f  and 
f '  ,  and  an  initial  approximation  <|) '  by 
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f  «)»  =  -q  +  j  In 


(1  +  sin  (JA 

(1  -  £  sin  (j)' 

Y  1 

_^1  -  sin  (|) J 

^1  +  £  sin  (j) 

JJ 

(7.30) 


£2 


^  (1-  £2  sin2  <j))cos  § 

=  2  arctan  (e^)  -  K/2 

Then  successive  iterations  for  (|>  are  found  from  the  equation 

f  ((l)1  ) 


(j)  •  •  =  (j)'  - 


f  1  «t>-  ) 


(7.31) 


In  particular,  substitute  (j)'  into  this  equation  to  obtain  (J)'1,  the 
first  approximation  to  0.  Then  set  (j>'  =  <)>'',  substitute  (J)'  into 
the  equation  to  obtain  (j)'  the  second  approximation  to  <|). 

Continue  this  procedure  until  the  differences  between  successive 
candidates  for  (|>  are  sufficiently  small. 

7.5.2  Lambert  Conformal  Conic  with  one  standard  parallel. 

In  the  second  case,  the  point  scale  factor  is  one  along  one 
standard  parallel,  (|)0.  Fix  a  central  meridian,  XQ. 

7. 5. 2.1  Finding  (x.v)  .  The  mapping  equations  are 


where 


x  =  r  sin  LA  (7.32) 

y  =  r0  -  r  cos  LA 


A  =  X  -  Xa 


L  =  sin  <|>0 


K  =  N0  cot  <|)0 

-Lq 

r  =  K  e 


<Io  sin(t>o 
e 


(7.33) 


rD  =  K  e 


-Lqo 


e  =  base  of  natural  logarithms  (  In  e  =  1  ) , 


N  = 


Vl  -  £2  sin2  4> 


q  =  In 


tan  ff  + 

f  1  -  £  sin  $ 

e 

2 

U  2)  \ 

1  +  £  sin  (() 

qo  =  q  evaluated  at  (j)0 
N0  =  N  evaluated  at  <])0 
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7. 5. 2. 2  Finding  (()).A)  with  one  standard  parallel.  To 
transform  in  the  other  direction,  suppose  that  x,  y,  <J)0 ,  and  A0 
are  given.  First,  compute  L,  K,  N0,  go,  r,  and  rQ  from  equation 
set  7.33.  Then,  as  in  the  previous  case,  (<|>,  A,)  can  be  found  from 


A  = 


K 


(7.34) 


in  which 


r 


-  y 

cos  0 


0  =  arctan 

¥ 


(7.35) 


Again,  <])  can  be  obtained  from  q  by  using  the  Newton-Raphson  method 
(see  equation  sets  7.30  and  7.31). 


7.5.3  Scale  factor  and  convergence.  The  scale  factor  and 
convergence  of  the  meridian  are  given  by 


k  = 


rL 

N  cos  (J) 


Y  =  lA 


(7.36) 


In  the  case  of  two  standard  parallels,  the  values  of  r,  L,  and  N 
are  given  in  equation  set  7.27.  For  one  standard  parallel,  the 
values  of  r,  L,  and  N  are  found  in  equation  set  7.33. 

7.5.4  Accuracy .  For  the  Lambert  conformal  conic  projection 
with  either  one  or  two  standard  parallels,  the  equations  for  x,  y, 
A,  q,  k,  and  Y  are  exact.  The  equation  for  (j)  can  be  iterated 
until  the  desired  accuracy  is  obtained. 


7.5.5  Area  of  coverage.  The  family  of  Lambert  conformal 
conic  projections  and  their  limiting  cases  is  used  worldwide.  For 
a  given  member  of  this  family,  the  point  scale  factor  increases  as 
a  point  moves  away  from  the  standard  parallel (s). 
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7.6  £olar_SJter.eographic  oroi ect ion .  The  Polar 
Stereographic  projection  is  a  limiting  case  of  the  Lambert 
conformal  conic  projection  when  the  one  standard  parallel 
approaches  a  pole.  In  this  conformal  projection  meridians  are 
straight  lines,  parallels  are  concentric  circles,  and  the  point 
scale  factor  is  one  at  the  pole  (see  Figure  7.5) . 


Figure  7.5.  Meridians  and  parallels  in  the  Universal  Polar 
Stereographic  Projection.  Left:  South  Zone, 
Right :  North  Zone  (Not  drawn  to  scale . ) 


7.6.1  Polar  Stereoaraohic  mapping  equations.  There  is 
considerable  diversity  in  the  treatment  that  various  authors  give 
to  the  polar  Stereographic  projection.  Some  authors  interchange 
the  x  and  y  axis;  others  change  the  sign  convention.  The 
discussion  here  is  largely  taken  from  NSWC/DL  TR-3624  (1977)  since 
it  most  closely  conforms  to  DMA  TM  8358.2  (1989). 

7. 6. 1.1  Finding  (x.v).  Given  (<}),  A,)  ,  the  Polar 
Stereographic  coordinates  are 

x  =  r  sin  X  {1 .31) 

y  =  -r  cos  X  (Northern  hemisphere) 
y  =  r  cos  X  (Southern  hemisphere) 


in  which 


r  =  K  tan 


(n  \ty\\ 

(  1  +  £  sin  |<1>|  \ 

U  "  2  J 

^1-6  sin  |<])| 

e 

2 


(7.38) 
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|<])l  denotes  the  absolute  value  of  (j>. 

(<}>,  A,)  is  expressed  in  radians. 

At  the  poles.  A,  is  not  defined.  Nevertheless,  these  equations  can 
be  used  at  the  poles  since  r  is  zero  there.  At  <j)  =  ±  n/2,  pick  an 
arbitrary  value  for  A,.  Then,  equation  sets  7.37  and  7.38  yield 
Polar  Stereographic  coordinates  of  (0,  0) . 

7. 6.1.2  Finding  ((|).A.).  The  discussion  of  the  conversion, 
from  (x,  y)  to  (<j),  X)  follows  Krakiwsky  (1973).  If 
(x,  y)  =  (0,  0),  then  <J)  =  ±  90°  and  X  is  undefined.  Otherwise, 
the  longitude  is  given  by 

X 
X 

If  (x,  y)  &  (0,  0),  then  x,  y  or  both  are  non-zero,  use  one  of 
these  non-zero  values  to  find  q  from  the  equations 


X 

=  -  arctan  —  (Northern  hemisphere) 


(7.39) 


=  arctan  —  (Southern  hemisphere) 


where 


_ LzJ _ 

K  I  cos  A,  I 

X 

K  sin  A, 


when  |xl  >  lyl 
when  lyl  >  lx  I 


(7.40) 


K 


2  a2 

b 


(i  -  e  A 


E 

2 


v1  +  E, 


(7.41) 


As  in  the  Lambert  conformal  conic  projection,  <j)  can  be  found  from 
q  using  the  Newton-Raphson  method.  Define  functions  f  and  f‘,  and 
an  initial  approximation  4>o  by 


f  (4>) 


-q 


In 


(1  +  sin  <]A 

fl  -  £  sin  (j)^ 

E 

_^1  -  sin  <j)^ 

^1  +  6  sin 

f  • (4» 


_ l  -  e2 _ 

(1-  £2  sin2  <]))  cos  <|) 


q 

2  arctan  e 


n 

2 


(7.42) 
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Then,  successive  iterations  for  0  are  found  from  the  equation 


<l>n  =  4>„-i  - 


f  '  <4>n-l> 


(7.43) 


In  particular,  substitute  (j)0  into  this  equation  to  obtain  the 
first  approximation  to  0.  Then  substitute  <}>!  into  the  equation  to 
obtain  <j)2,  the  second  approximation  to  (|).  Continue  this  procedure 
until  the  differences  between  successive  candidates  for  <J)  are 
sufficiently  small.  Even  though  q  approaches  infinity  at  the 
poles,  this  method  is  satisfactory  for  many  computers.  This 
procedure  yields  a  positive  value  of  <j).  In  the  southern 
hemisphere,  replace  this  positive  value  of  <j)  by  -(J). 

7. 6. 1.3  Alternate  method  for  finding  (ft.  Alternatively,  an 
approximate  equation  for  <})  is  given  in  Geological  Survey 
Professional  Paper  1395  (1987) . 


(|>  =  3C  +  A  sin2%  +  B  sin4%  +  C  sin6%  +  D  sin8% 


(7.44) 


where 


n 

X  -  2  -  z 


(in  radians) 


(7.45) 


K  sin 


f-l 


x  *  0 


K  cos 


I  y  *  o 


2a2  (1  -  £ 
b  i  +  8 


A  =  2*2  *  2^'  *  S*6  *  360*° 


B  .  Xg.  + 

a  ~  Ur  2401 


11520 


7  6  81 


C  "  1208  +  1120* 


_  _4279_^8 
U  161280° 


£2  = 


a2  -  b2 


53 


Since  this  procedure  yields  a  positive  value  of  <J),  replace  (J)  by  -<]) 
in  the  southern  hemisphere. 

7. 6. 1.4  Finding  the  point  scale  factor.  The  point  scale 
factor  is  given  by 


k  =  1  at  the  poles 


(7.46) 


where 


k  = 


r  =  K  tan 


N  cos  <j) 

7L  M 


elsewhere 


(?- 


K  = 


N  = 


2a2  (i  -  e" 


^  +  e, 
a 


Vl  -  £2  sin2  (J) 


1  +  £  sin  !(})  |  \ 
1  -  £  sin  id)  I 


J 


7. 6. 1.5  Finding  the  convergence  of  the  meridian, 
convergence  of  the  meridian  is 


(7.47) 


The 


Y=A,  (Northern  hemisphere)  (7.48) 

Y=  -X  (Southern  hemisphere) 

7.6.2  Accuracy .  For  the  Polar  Stereographic  projection, 
the  equations  for  x,  y,  X,  and  q  are  exact.  The  equation  for  <j) 
can  be  iterated  until  the  desired  accuracy  is  obtained. 

7.6.3  Area  of  coverage.  There  is  no  general  agreement  on 
the  area  of  coverage  for  the  Polar  Stereographic  projection. 
Section  7.7.5  gives  the  area  of  coverage  for  UPS  coordinates . 


7.7  Universal  Polar  Stereographic  (UPS)  projection.  The 
Universal  Polar  Stereographic  (UPS)  projection  is  the  standard 
military  grid  used  in  polar  regions  (see  Figure  7.6).  The  main 
reference  for  the  UPS  grid  is  DMA  TM  8358.2  (1989).  The  UPS  grid 
is  a  family  of  two  projections  that  differ  from  the  Polar 
Stereographic  projection  in  these  ways: 


a.  Both  the  x  and  y  values  (see  Figure  7.6),  called  easting 

and  northing,  respectively,  have  an  origin  of  2,000,000 
meters . 

b.  The  scale  factor  at  the  origin  is  0.994. 
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Figure  7.6. 


Meridians  and  parallels  imposed  on  a  Universal 
Polar  Stereographic  grid  (north  zone) . 


The  limits  of  the  system  are  north  of  84°  N  and  south  of  80°  S.  In 
order  to  provide  a  30 -minute  overlap  with  the  UTM  grid,  the  UPS 
grid  contains  an  overlap  that  extends  to  83°  30'N  and  79°  30' S. 
Although  the  UPS  grid  originally  was  defined  to  use  the 
International  Ellipsoid,  in  practice  several  reference  ellipsoids 
have  been  used.  The  ellipsoid  appropriate  to  the  selected  datum 
should  be  used. 


7.7.1  Universal  Polar  Stereograohic  (UPS)  mapping 
equations .  For  known  geodetic  coordinates,  the  UPS  easting  is 
given  by 

xups  =  2,000,000  +  0 . 994r  sin  X  (7.49) 

and  the  UPS  northing  is  given  by 

Yups  =  2,000,000  -  0.994r  cos  X  (Northern  hemisphere)  (7.50) 

yups  =  2,000,000  +  0.994r  cos  X  (Southern  hemisphere) 

The  parameter  r  is  computed  using  equation  set  7.38. 

7.7.2  Finding  (<j).X,).  If  a  hemisphere  and  UPS  coordinates 
(x,  y)  are  known,  convert  the  UPS  values  to  Polar  Stereographic 
coordinates  using  equation  set  7.51,  and  then  transform  to 
geodetic  coordinates  as  described  in  Sections  7. 6. 1.2  and  7. 6. 1.3. 
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^■ps  —  q  994  ^ups  2,000,000) 


(7.51) 


_  1 

Yps  “  0.994 


(y  ups 


2,000,000) 


In  particular,  at  the  poles  where  x  =  2,000,000  and  y  =  2,000,000 
then,  (j)  =  ±  90°  and  A  is  undefined.  Otherwise  longitude  and 
latitude  can  be  found  from 


A  = 
A  = 


-  arctan 

+  arctan 


(XUPS 

(Yups 

^xups 

(Yups 


2,000,000) 

2,000,000) 

2,000,000) 

2,000,000) 


( Northern  hemi sphere )  (7.52) 
(Southern  hemisphere) 


(|)  =  %  +  A  sin2%  +  B  sin4%  +  C  sin6%  +  D  sin8% 


in  which 


X  = 


n 

2 


z 


z 

tan  2 


(xups  -  2,000,000) 
0.994K  sin  A 


xups  *  2,000,000 


(7.53) 


Z  (y ups _ 2 , 000 , 000 ) 

2  0.994K  cos  A 

Yups  96  2,000,000  (Northern  hemisphere) 

z  _  (ynpfi  -  2,000,000) 

2  “  +  0.994K  cos  A 

Yups  ^  2,000,000  (Southern  hemisphere) 


K 


2a2  (1  -  E)  f 
b  1^1  +  ej 


£2  = 


a2  -  b2 
a2 


The  values  of  A,  B,  C,  and  D  are  found  using  equation  set  7.45. 


In  the  southern  hemisphere,  replace  the  computed  value  of  <|)  with 

-0. 
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7.7.3  Finding  point  scale  factor  and  convergence  of  the 
meridian.  The  point  scale  factor  and  convergence  of  the  meridian 
are 


k  =  0.994  at  the  poles  (7.54) 


k  =  0.994  - - t  elsewhere 

N  cos  (p 

y  =  X  (Northern  hemisphere) 
y  =  -X  (Southern  hemisphere) 

where  r  and  N  are  given  in  equation  set  7.47. 

7.7.4  Accuracy .  The  accuracy  of  the  UPS  equations  is  based 
upon  the  accuracy  of  their  Polar  Stereographic  counterparts.  For 
the  Polar  Stereographic  projection,  the  equations  for  x,  y,  X,  and 
q  are  exact.  The  equation  for  (J)  can  be  iterated  until  the  desired 
accuracy  is  obtained. 

7.7.5  Area  of  coverage.  The  area  of  coverage  for  UPS 
coordinates  is  defined  by  zone  limits  and  overlap. 

North  zone:  The  north  polar  area  84°  -  90° 

South  zone:  The  south  polar  area  80°  -  90° 

UTM  overlap:  30'  overlap, 

North:  83°  30' 

South:  79°  30‘ 

7.8  The  U.S.  Military  Grid  Reference  System  (MGRS) .  The 
MGRS  is  an  alpha-numeric  system,  based  upon  the  UTM  and  UPS  map 
projections,  for  identifying  positions.  An  MGRS  coordinate 
consists  of  a  zone  designation,  alphabetic  100,000-meter  grid 
square  designator,  and  numeric  coordinates  within  the  100,000 
meter  grid  square.  MGRS  coordinates  are  defined  for  the  primary 
UTM  and  UPS  areas,  but  not  for  the  overlap  areas. 

7.8.1  MGRS  coordinates  in  the  UTM  area.  In  the  UTM  area,  MGRS 

coordinates  are  based  on  the  ellipsoid,  geodetic  latitude,  UTM 
zone,  easting  and  northing  as  follows: 

a.  The  first  two  characters  of  the  MGRS  coordinate  are  the 
two  digits  of  the  numeric  UTM  zone  (see  sections  7. 4. 1.1  and 

7. 4. 1.2.)  Leading  zeros  must  be  included. 

b.  The  third  character  of  the  MGRS  coordinate  is  a  letter 
representing  a  band  of  geodetic  latitude.  Beginning  at  80°  South 
and  proceeding  northward,  20  bands  are  lettered  C  through  X, 
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omitting  I  and  0.  The  bands  are  all  8°  high  except  band  X  (72°N  to 
84°N) ,  which  is  12°  high. 

c.  The  fourth  and  fifth  characters  of  the  MGRS  coordinate 
are  a  pair  of  letters  representing  the  100,000-meter  grid  square. 
The  letter  pair  can  be  selected  from  figure  7.7  and  7.8  (from  DMA 
TM  8358.1.)  First,  reduce  the  UTM  northing  by  multiples  of 
2,000,000  until  the  resulting  value  is  in  the  range  of  0  to 
1,999,999  meters.  Second,  reduce  the  UTM  northing  and  easting 
values  to  the  nearest  100,000  meters.  Third,  locate  the 
corresponding  100,000  meter  northing  and  easting  grid  lines  for 
the  UTM  zone  number  on  the  figure.  The  grid  square  identifier 
will  be  immediately  above  and  to  the  right  of  the  intersection  of 
the  easting  and  northing  grid  lines  on  the  figure. 

Note:  The  alphabetic  method  of  designating  100,000- 
meter  grid  squares  has  changed  over ' time .  Older  products 
may  have  different  grid  square  designations  than  those 
described  above.  Also,  some  software  was  programmed  to 
compute  MGRS  coordinates  based  on  geographic  position  and 
a  map  of  preferred  ellipsoids.  The  boundaries  between 
preferred  ellipsoids  have  changed  with  time.  Therefore, 

MGRS  coordinates  computed  with  such  software  may  not  agree 
with  current  map  products.  MGRS  computation  software 
should  always  require  input  of  the  ellipsoid.  If  MGRS 
coordinates  are  to  be  used  in  conjuction  with  a  map 
product  they  should  be  checked  against  the  map  product 
to  verify  compatibility. 

d.  The  remainder  of  the  MGRS  coordinate  consists  of  the 
numeric  easting  and  northing  values  within  the  100,000-meter  grid 
square.  The  left  half  of  the  digits  is  the  easting  grid  value, 
which  is  read  to  the  right  of  the  100,000-meter  easting  grid  line 
established  in  step  c.  above.  The  right  half  of  the  digits  is  the 
northing  grid  value,  which  is  read  northward  of  the  100,000-meter 
northing  grid  line  established  in  step  c.  above.  Both  the  easting 
and  northing  grid  values  are  within  the  range  of  0  to  100,000 
meters  (i.e.,  zone  and  grid  square  designators,  only)  to  1  meter 
(i.e.,  zone  and  grid  designators  followed  by  10  digits,  five  for 
easting  grid  value  and  five  for  northing  grid  value) .  Both 
easting  and  northing  value  must  have  the  same  resolution  and  must 
include  leading  zeros. 

7.8.2  MGRS  coordinates  in  the  UPS  area.  In  the  UPS  area, 
MGRS  coordinates  are  based  on  the  ellipsoid,  geodetic  latitude  and 
longitude,  easting,  and  northing  as  follows: 

a.  The  first  character  of  the  MGRS  coordinate  is 
A  in  the  Southern  and  Western  Hemispheres, 

B  in  the  Southern  and  Eastern  Hemispheres, 

Y  in  the  Northern  and  Western  Hemispheres,  and 

Z  in  the  Northern  and  Eastern  Hemispheres  as  shown  in  figures  7.9 
and  7.10. 
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b.  The  second  and  third  characters  of  the  MGRS  coordinate  . 
are  letters  representing  a  100,000  meter  grid  square.  The  letter 
pair  can  be  selected  from  figure  7.9  or  7.10  (from  DMA  TM  8358.1.) 
Reduce  the  UPS  easting  and  northing  values  to  the  nearest  100,000 
meters.  Locate  the  corresponding  100,000  meter  northing  and 
easting  grid  lines  for  the  UTM  zone  number  on  the  figure.  The 
grid  square  identifier  will  be  immediatly  above  and  to  the  right 
of  the  intersection  of  the  easting  and  northing  grid  lines  on  the 
figure. 


c.  The  method  for  calculating  the  remainder  of  the  MGRS 
coordinate  is  the  same  as  in  paragraph  7.8.1  d. 

7.9  The  World  Geographic  Reference  System  (GEOREF) .  The 
GEOREF  is  an  alpha-numeric  system  for  reporting  positions  for  air 
defense  and  strategic  air  operations.  It  is  based  upon  geodetic 
coordinates.  GEOREF  is  described  in  DMA  TM  8358.1  (1990) . 

7.10  Non-standard  grids.  The  specifications  and  diagrams  of 
non-standard  grids,  used  as  the  primary  or  secondary  grid  on  DMA 
produced  maps,  are  contained  in  DMA  TM  8358.1  (1990). 

7.11  Modifications  to. .map  projections.  Providing  equations 
for  transforming  between  geodetic  coordinates  and  every  local  map 
projection  is  beyond  the  scope  of  this  Technical  Report.  However, 
many  local  map  projections  are  minor  modifications  of  the  map 
projections  given  in  this  Technical  Report.  Typically,  these 
modifications  include  use  of  false  Eastings  (E)  and  Northings  (N) 
at  the  projection  origin,  a  non-unity  scale  factor  (kQ)  along  the 
principal  line  of  the  projection,  and  units  other  than  meters. 
Sometimes  the  local  projection  origin  is  offset  from  the  origin 
used  in  this  Technical  Report.  Chapter  four  of  DMA  TM  8352.1 
(1991)  gives  parameters  for  several  local  projections. 

If  the  local  projection  is  based  on  one  of  the  projections 
given  in  this  Technical  Report  convert  geodetic  latitude  (<j))  and 
longitude  (A,)  to  coordinates  in  the  local  projection  by  using  the 
following  projection; 

1.  Make  sure  (<j).  A)  are  on  the  ellipsoid  for  the  local 
projection  and  determine  the  ellipsoid's  parameters,  a 
and  either  f  or  b. 

2.  Convert  (<|),  A.)  to  map  projection  coordinates  (Xp,  Yp) 
using  the  appropriate  map  projection  equations  given  in 
this  Technical  Report.  Note  that  the  units  for  Xp  and  Yp 
are  in  meters. 

3.  If  the  local  projection  has  false  coordinates,  E,  N,  at 
the  origin,  convert  them  to  units  of  meters  if  necessary. 
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4.  The  map  projection  coordinates  (Xp,  Yp)  and  scale  factor 
(kp)  can  be  converted  to  local  projection  coordinates 
using  equations  7.55  if  the  origin  of  the  local 
projection  is  the  same  as  the  origin  of  that  projection 
as  used  in  this  technical  report.  Those  origins  are 
central  meridian  (XQ)  and  equator  for  Transverse 
Mercator  and  Mercator  projections,  and  standard 
parallel (s)  and  equator  for  the  Lambert  Conformal 
projection. 

X  =  kcXp  +  E 

Y  =  kcYp  +  N  *  (7.55) 

k  =  k0kp 

Y  =  Yp 


where 

E  =  false  X-value  (in  meters)  at  the  origin 
N  =  false  Y-value  (in  meters)  at  the  origin 

ko  =  scale  factor  along  the  appropriate  line. 

(  Along  the  equator  for  the  Mercator  projection, 
along  the  standard  parallel (s)  for  the  Lambert 
Conformal  Conic  projection,  and  along  the  central 
meridian  for  the  Transverse  Mercator  projection.) 

5.  Convert  (X,  Y)  to  the  units  that  are  locally  used.  It 
is  the  user 1 s  responsibility  to  determine  if  this 
procedure  is  valid.  This  can  be  done  by  comparing 
computed  mapping  coordinates  to  test  points  where  the 
relationship  between  ((|),  X)  and  (X,  Y)  is  known, 
independently. 

The  above  procedure  is  intended  to  be  used  only  when  the 
origin  of  the  local  map  projection  coincides  with  the  origin,  used 
in  this  Technical  Report,  some  Transverse  Mercator  projections 
have  origins  of  (<|)o,  Xo)  instead  of  (0,  Xo)  on  the  equator.  For 
some  of  these  projections,  agreement  with  test  points  can  be 
obtained  by  changing  the  computation  of  Y  in  equation  set  7.55  to 
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Y  =  koYTM  +  N  -  koS(|>0  (7.56) 

where 

ko  =  scale  factor  on  the  central  meridian,  Xo 

Ytm  =  Transverse  Mercator  Y  value  computed  from  equation 
set  7.8 

N  =  false  Northing  in  meters  of  the  origin 

S«J>0  =  value  computed  using  equation  set  7.9 

(arc  length  in  meters  on  the  ellipsoid  from  the 
equator  to  the  latitude  <|)o.) 
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8.  SELECTION  OF  TECHNIQUES.  SOFTWARE  DEVELOPMENT.  AND  TESTING 


8.1  The  general  conversion  process.  The  general  coordinate 
conversion  process  can  be  divided  into  three  basic  parts: 

1.  Conversion  between  map  projection  coordinates  and 
geodetic  coordinates  (see  Figure  8.1). 

2.  Coordinate  conversion  between  geodetic  coordinates  and 
Cartesian  coordinates  within  a  datum  (see  Figure  8.2). 

3.  Coordinate  transformation  between  datums  (datum  shifts) 
(see  Figure  8.3) . 

The  general  conversion  process  which  includes  all  methods 
presented  in  this  Technical  Report,  is  described 
graphically  in  Figure  8.4. 

a.  Figure  8.1  represents  the  map  projection  coordinate  (x,y) 
to  the  geodetic  coordinate  (<j),A,)  conversion  described  in 
Section  7.  As  noted  in  Section  7.1,  map  projection 
coordinates  do  not  include  the  height  information 
necessary  for  finding  the  geodetic  height  (H)  or 
elevation  (h) .  For  conversions  from  map  projection  to 
geodetic  coordinates,  the  user  may  have  to  provide  this 
information.  Elevations  can  often  be  inferred  from 
contour  lines  on  the  map  or  they  can  be  obtained  from 
DTED  or  direct  measurement.  Geodetic  height  may  be 
measured  by  using  GPS  or  by  calculating  from  an  elevation 
and  geoid  separation,  or  for  some  applications,  may  be 
assumed  to  be  0  (zero) . 


A 

x,y 


Figure  8.1.  Map  projection  to  geodetic  coordinate 

conversion.  Height  information  for  H  and  h  is 
not  included  in  the  conversion. 


where 

x,y  =  map  projection  coordinates 
<))  =  geodetic  latitude 
X  =  geodetic  longitude 

A  =  map  projection  coordinates  to/from  latitude, 
longitude  conversion 
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b.  Figure  8.2  represents  coordinate  conversions  within  a 
datum  as  described  in  Section  4.3.4  and  6.5. 

x,  Y,  Z 

c  | 

o,  A,  h 

B  ^ 

O,  X,  H 

FIGURE  8.2.  Coordinate  conversions  within  a  datum, 
where 

Cartesian  coordinates 

geodetic  coordinates  (latitude,  longitude, 
geodetic  height) 
elevation 

geodetic  height  to/ from  elevation  conversion 
geodetic  coordinate  to/ from  Cartesian  coordinate 
conversion 

Note  that  a  conversion  between  Cartesian  coordinates  (X,Y,Z) 
and  local  datum  coordinates  ((|),  A,, H)  is  not  a  direct  conversion 
and  requires  going  through  the  geodetic  coordinates  (<j),A,,h)  . 

c.  Figure  8.3  represents  the  datum  shift  portion  of  the 
general  conversion  process  as  described  in  Section  6. 

Two  attributes  of  figure  8.3  are  important  to  note. 

1.  Direct  local  datum  to  local  datum  shifts  are  not 
performed. 

2.  Elevations  do  not  change  when  changing  datums. 
Exceptions  to  the  second  statement  are  explained 
in  Section  6.7. 

d.  In  general,  a  direct  local  datum  to  local  datum  shift 
should  not  be  done.  Instead,  the  local  datum  to  local 
datum  shift  should  be  done  in  two  steps. 

1.  An  (old)  local  datum  to  WGS  84  shift. 

2.  A  WGS  84  to  (new)  local  datum  shift. 


X,  Y,  Z  = 

<|>,  X,  h  = 

H  = 
B  = 
C  = 
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Figure  8.3.  Datum  shifts. 

Note  elevation  values  (H)  do  not  generally  change  when 
changing  datums, 

where 

X,  Y,  Z  =  Cartesian  coordinates 

<j),  X,  h  =  geodetic  coordinates  (latitude,  longitude, 
geodetic  height) 

H  =  elevation 

C  =  geodetic  coordinate  to/ from  Cartesian  coordinate 
conversion 

AA  =  Cartesian  coordinate-based  datum  shifts 
BB  =  geodetic  coordinate-based  datum  shifts 
CC  =  elevations  are  unchanged 

e.  The  general  conversion  process,  which  includes  all 

methods  presented  in  this  Technical  Report,  is  described 
graphically  in  figure  8.4.  Figure  8.4  is  the  combination 
of  figures  8.1  through  8.3, 
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h 


Figure  8.4.  The  general  conversion  process. 


where 


x,  y 
X,  Y,  Z 

<j>,  X,  h 

H 

A 

B 

C 

AA 

BB 

CC 


map  projection  coordinates 
Cartesian  coordinates 

geodetic  coordinates  (latitude,  longitude, 

geodetic  height) 

elevation 

map  projection  coordinates  to/ from  latitude, 
longitude  conversion 

geodetic  height  to/from  elevation  Conversion 
geodetic  coordinate  to/ from  Cartesian  coordinate 
conversion 

Cartesian  coordinate-based  datum  shifts 
geodetic  coordinate-based  datum  shifts 
elevations  are  unchanged 


f.  In  figures  8.1  through  8.4,  coordinate  conversions  are 
denoted  as  a  single  letter  (A,  B,  C) .  They  are  described 
in  Section  7,  Section  6.5  and  Section  4.3.4, 
respectively.  Datum  shifts  are  denoted  as  double  letters 
(AA,  BB,  CC)  and  are  described  in  Section  6. 


69 


g.  The  Cartesian  coordinate-based  datum  shifts  (AA)  include: 

Seven-parameter  geometric  transformation  (6.2) 

Ax, Ay,Az  direct  application  (6.3.1  &  6.3.2) 

Both  of  these  datum  shift  methods  are  geometric 
transformations . 

h.  The  geodetic  coordinate  based  datum  shifts  (BB)  include: 

Molodensky  shifts  (6.3.3) 

Ax,  Ay,  Az  three-step  method  (6. 3. 1.1  &  6. 3. 2.1) 

Direct  WGS  72  to  WGS  84  transformation  (6.4.2) 

i.  For  the  purposes  of  this  Technical  Report,  the  elevation 
to  elevation  shift  (CC)  is  the  equality  function  (=); 
that  is,  the  elevation  for  the  same  location  on  different 
datums  is  considered  to  be  the  same  (H  =  H)  .  An 
exception  to  this  rule  is  the  WGS  72 -derived  elevations 
to  WGS  84  elevations  shift.  Section  6.7  describes  the 
elevation  to  elevation  shift  process. 

8.1.1  Application  of  the  general  conversion  process.  The 
general  conversion  process  described  above  can  be  used  to  solve 
coordinate  conversion/datum  shift  problems  in  five  steps: 

a.  Locate  beginning  point  and  ending  point  of  the  problem  on 
figure  8.4. 

b.  Determine  valid  paths  between  beginning  and  end  points. 

c.  Evaluate  alternatives  and  select  a  path. 

d.  Locate  the  needed  equations  and  parameters  in  this 
Technical  Report. 

e.  Implement  the  selected  path. 

8.1.2  Procedural  examples. 

8. 1.2.1  Example  8.1.  Horizontal  geodetic  coordinate-based 
datum  shifts  ( (<|),A,,h)LGS  to  (<]),  X,)WGS  84)  (Local  datum  geodetic 
heights  (hLGS)  typically  are  not  available.  It  is  not  recommended 
they  be  transformed  to  WGS  84) . 

This  Technical  Report  has  discussed  two  methods  of  doing 
datum  shifts.  These  methods  are  the  three-step  method,  and 

Molodensky  shifts.  Either  the  three-parameter  Ax,  AY,  Az  shift  or 
the  seven -parameter  shift  may  be  used  in  the  three-step  method. 
This  example  corresponds  to  a  “BB"  shift  of  figures  8.3  and  8.4. 
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a.  Approach  1_-  Three-step  method  using  a  seven-parameter 

gfrif.t  t, 

1.  Compute  (X,Y,Z)LGS  from  (()>,  A,  h)  LGS  using  equation  sets 
4.4  and  4.5,  and  the  reference  ellipsoid  parameters 
from  Table  A.l. 

2.  Compute  (X,Y,Z)WGS84  from  (X,  Y,  Z)  Lgs  using  the  seven- 
parameter  shift  equation  set  6.1  and  shift  parameters 
from  Table  E.l. 

3.  Compute  { (J),  A)WGS  84  from  (X,Y,Z)WGS  84  using  equation 
sets  4.6  through  4.10. 

b.  Approach  2  -  Three- step,  method  using  Ax,  AY,  Az  shifts 

1.  Compute  (XjY/ZJlqs  from  (<J>,  A,  h)  using  equation  sets 
4 . 4  and  4.5.  Reference  ellipsoid  parameters  are 
found  in  Table  A.l. 

2.  Compute  (X,Y,Z)WGS84  from  (X,  Y,  Z)  Lgs  using  equation 
set  6.3,  and  shift  parameters  from  Table  B.l 

3.  Compute  (<}>,  A)Wgs  84  from  (X,Y,Z)WGS  84  and  equation 
sets  4 . 6  through  4.10. 

c.  Approach  3  -  Using  Molodenskv  shift 

1.  Compute  A(j)  and  AA  from  (<J),  A,  hLGS)  using  the 
Molodensky  equation  set  6. 5, shift  parameters  from 
Table  B.l,  and  the  local  datum's  reference  ellipsoid 
parameters  given  in  Table  A.l. 

2.  Compute  (<]),  A)  WGS  84  from  ((j),  A,)  LGS  using  equation 

-  set  6.6. 

8. 1.2. 2  Example  8.2.  WGS  84  Cartesian  coordinate  to  local 
latitude,  longitude,  elevation. 

Several  paths  for  converting  (X,Y,Z)Wgs84  to  (<)>lgs  ,  Algs,H)  can 
be  traced  in  figure  8.4.  The  steps  of  three  approaches  are 
provided  below. 

a.  Approach  1  -  Using  a  seven-parameter  datum  shift 

1.  Compute  (XjY/ZJlqs  from  (X,  Y,  Z)WGS  84  using  the  seven- 
parameter  shift  equation  set  6.2  and  parameters  from 
Table  E.l. 
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2.  Compute  ( <|),  X,)LGS  from  (X,Y,Z)LGS  using  equation  sets 
4.6  through  4.10  with  the  parameters  of  the  local 
ellipsoid. 

3.  Compute  hWGS  84  from  (X,Y,Z)WGS  84  using  equations  4.7 
through  4.11  (Use  equation  4.12  in  polar  regions). 

4.  Compute  H  from  hWGS  84  and  NWGS  84  using  equation  5.1, 
resulting  in  <|)LGS,  Xgs'  H. 

b.  Approach  2  -  Using  a  direct  AX,  Ay,  AZ  datum  shift 

1.  Compute  (X,Y,Z)lgs  from  (X,Y/Z)WGS84  using  the  AX, 

AY,  AZ  shift  equation  set  6.4  and  shift  parameters 
from  Table  B.l. 

2.  Compute  ( <j>,  X, ) LGS  from  (X/Y/Z)^  using  equation  sets 
4.6  through  4.11  with  the  parameters  of  the  local 
ellipsoid. 

3.  Compute  h^g  84  from  (X,Y ,Z)WGS  84  using  equations  4.7 
through  4.11.  (Use  equation  4.12  in  polar  regions). 

4.  Compute  H  from  hWGS  84  and  NWGS  84  using  equation  5.1, 
resulting  in  <j)LGS/  Xgs/  H. 

c.  Approach  3  -  Using  a  Molodenskv  datum  shift 

1.  Compute  (<j),  X,  h)WGs  84  from  (X,  Y,  Z)WGS  84  using  equation 
sets  4.6  through  4.11.  (In  polar  regions,  equation 
set  4.11  should  be  replaced  by  equation  set  4.12). 

2.  Compute  (([),  X,)LGS  from  (<|),  X,)Wgs  84  using  the 
Molodensky  equation  set  6.5,  equation  set  6.7  and 
shift  parameters  from  Table  B.l. 

3.  Compute  H  from  hWGS  84  and  NWGS  84  using  equation  5.1, 
resulting  in  <|)LGs,  X,Lgs#  H. 

8. 1.2. 3  Example  8.3.  Local  latitude,  longitude,  elevation 
( <t>LGs '  ^lgs'  H)  to  WGS  84  latitude,  longitude,  elevation 

(0WGS  84'  K  GS  84'  ' 

Example  8.3  is  a  extension  of  example  8.1.  Example  8.3  shows 
the  recommended  methods  for  determining  geodetic  height  used  in 
each  of  the  geometric  datum  shift  methods.  The  steps  of  three 
alternate  approaches  are  provided  as  follows. 
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1.  Compute  (X,Y,Z)LGS  from  (<j),  X,  H)  LGS  using  equation  sets 
4.4  and  4.5. 

2.  Compute  (X,Y,Z)WGS84  from  (X,  Y,  Z)  LGS  using  the  seven 
parameter  shift  equation  set  6.1  and  shift  parameters 
from  Table  E.l. 

3.  Compute  «|>,  X)WGs  84  from  (X,  Y,  Z)WGS  84  using  equation 
sets  4.6  through  4.10. 

4.  Take  H  from  the  input  coordinates  resulting  in  <j)wGs  84> 


1.  Compute  (XjYjZJlqs  from  (<J),  X,  H)  lqs  using  equation  sets 
4 . 4  and  4.5. 

2.  Compute  (X/Y,Z)Wgs84  from  (X,  Y,  Z )  LGS  using  AX,  AY, 

AZ  shift  equation  set  6.3,  and  shift  parameters  from 
Table  B.l. 

3.  Compute  «)>,  A,)WGS  84  from  (X,Y,Z)WGS  84  using  equation 
sets  4 . 6  through  4.10. 

4.  Take  H  from  the  input  coordinate  resulting  in  §uGs  84> 


1.  Compute  ((|),  X,)wgs  84  from  (<j>,  A,,  H)LGS  using  Mo lodensky 
equation  set  6.5,  equation  set  6.6,  and  shift 
parameters  from  Table  B.l. 

2.  Take  H  from  the  input  coordinate  resulting  in  (^g  84, 

^WGS  84/  H. 

8.2  Choosing  a  datum  shift  method.  Several  factors  should 
be  considered  in  selecting  a  datum  shift  method.  These  factors 
include:  shift  parameter  availability,  timing  requirements, 

accuracy  requirements,  software  availability,  hardware 
limitations,  standardization  considerations  and  potential  future 
requirements.  If  future  requirements  are  considered  as  part  of 
the  system  development,  future  product  improvements  may  be 
minimized.  These  requirements  may  include  the  need  to  provide 
data  to  other  systems  and  to  meet  higher  accuracy  requirements. 

The  discussion  below  applies  to  a  wide  range  of  Army  applications. 
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Special,  high  accuracy  applications  are  outside  the  scope  of  this 
Technical  Report  (see  Section  1.2). 

The  three-parameter  shift  is  recommended  for  mapping  and  charting 
and  most  other  military  application.  For  Western  Europe  (EUR-M) , 
and  Ordnance  Survey  of  Great  Britain  1936  (OGB-M) ,  STANAG  2211 
specifies  the  seven-parameter  shifts.  The  use  of  MRE's  is  not 
recommended,  although  it  should  be  noted  that  STANAG  2211 
specifies  the  use  of  MRE’s  for  geodetic  and  geophysical 
applications  where  “data  is  available  in  geodetic  coordinates." 

Parameter  sets  for  Ax,  Ay,  Az  shifts  are  found  in  Appendix  B,  and 
seven-parameter  geometric  shifts  are  contained  in  Appendix  E.  Be 
sure  to  use  the  same  method  throughout  a  project,  identifying  the 
method  chosen  in  the  project  documentation  so  that  future  users  of 
the  data  can  replicate  or  backtrace  your  results. 


8.2.1  Parameter  availability.  For  many  datums,  only  mean- 
value  three-parameter  AX,  AY,  and  AZ  shifts  are  given  in  this 
Technical  Report.  Based  on  system  implementation  tradeoffs,  the 
system  developer  can:  apply  the  shifts  directly  to  Cartesian 
coordinates  (see  Sections  6.3.1  and  6.3.2);  use  the  three-step 
method  (see  Section  6. 3. 1.1  and  6. 3. 2.1);  or  apply  shifts  to 
geodetic  coordinates  using  the  Molodensky  approximations  (see 
Section  6.3.3).  For  other  datums,  some  combination  of  seven- 
parameter  geometric  shifts,  and  mean  and  regional  AX,  AY,  and  AZ 
shifts  are  available.  Select  a  parameter  set  from  those  that  meet 
or  come  closest  to  meeting  the  accuracy  requirements  of  the  output 
product.  Base  the  final  selection  on  standardization  requirements 
and  system  implementation  tradeoffs. 

8.2.2  Accuracy .  Errors  in  the  shift  parameters  are  the 
dominant  error  source  in  datum  shifting.  Errors  introduced  by  any 
of  the  datum  shift  algorithms  given  in  this  Technical  Report  are 
small  compared  to  the  errors  in  the  shift  parameters  and  can 
generally  be  neglected.  However,  take  care  that  accuracy  is  not 
degraded  by  the  software  and  hardware  implementation  of  these 
algorithms . 

a.  For  the  large  datums, the  regional  AX,  AY,  AZ  shifts  and 
seven -parameter  shifts  are  generally  considered  the  most 
accurate,  and  mean  value  AX,  AY,  AZ  shifts  have  the 
poorest  accuracy.  The  estimated  accuracies  of  the  shift 
parameters  are  given  in  Appendices  B  and  E  for  each 
listed  datum.  The  errors  in  geodetic  coordinates  may  be 

approximated  from  the  AX,  AY,  AZ  errors  listed  in 
Appendix  B,  using  the  equation  set  6.8. 

b.  The  error  estimates  and  general  guidance  above  must  be 
used  judiciously.  For  a  given  datum,  all  types  of  shift 
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parameters  were  derived  from  essentially  the  same  set  of 
control  data,  and  errors  were  estimated  by  using  a 
comparison  to  the  same  data.  For  small  datums,  with 
poorly  distributed  control,  the  better  fitting  shift  type 
may  not  be  more  accurate.  The  seven-parameter  model 
cannot  accurately  represent  datum  rotations  and  scale 
changes  unless  the  control  is  geometrically  well 
distributed  across  the  datum. 

c.  Regional  AX,  AY,  AZ  shifts  were  developed  when  there  were 
sufficient  data  to  identify  regional  variations  over  the 
datum  as  a  whole.  Regional  AX,  AY,  AZ  shifts  are  more 
'  accurate  than  mean-value  AX,  AY,  AZ  shifts. 

8.2.3  Implementation  issues.  System  implementation  issues 
may  influence  the  selection  of  transformation  method (s) . 

a.  Computational  speed  and  efficiency  are  important  for  some 
applications,  such  as  real-time  navigation.  The  shift 
type  and  method  might  be  selected  from  those  that  meet 
accuracy  requirements  to  minimize  computational  time. 
Alternately,  shifts  might  be  computed  in  a  background 
process  and  applied  as  fixed  values  for  a  period  of  time 
or  over  a  given  area.  The  system  designer  must  analyze 
the  impact  of  such  design  strategies  on  system  accuracy 
and  performance. 

b.  For  some  applications,  system  computational  and  data 
storage  capability  are  considerations  in  selecting  and 
implementing  transformations.  The  system  designer  must 
determine  that  word  precision,  range,  and  accuracy  of 
mathematical  functions  preserve  the  required  accuracy. 

c.  Code  complexity,  modularity,  and  software  maintenance  are 
important  system  design  concerns.  While  the  datum  shift 

,  algorithms  themselves  are  relatively  simple,  selection 
logic  often  is  complex. 

d.  Both  future  requirements  and  present  needs  should  be 
considered.  Factors  for  consideration  include:  timing 
requirements,  current  accuracy  requirements,  potential 
future  accuracy  requirements  (such  as  requirements  to 
provide  data  for  future  systems),  effects  of  improved 
databases,  and  the  effects  of  improvements  in  other- 
inputs.  In  other  words,  the  savings  in  future  planned 
product  improvements  should  be  evaluated  against  today's 
costs. 

e.  It  is  possible  that  a  software  developer  may  choose  to 
use  a  transformation  technique  that  is  more  accurate  than 
is  required.  In  such  a  case,  the  results  are  not 
necessarily  superior  to  those  arising  from  a  less 
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accurate,  but  acceptable,  technique.  The  data  are  simply- 
degraded  less  by  the  more  accurate  method;  the 
transformed  data  can  be  no  better  than  the  input  data. 

For  example,  techniques  that  achieve  geodetic  accuracy 
must  be  used  when  geodetic  accuracy  is  needed.  Such 
techniques  can  be  used  with  a  1 : 250, 000-scale  product, 
but  the  transformed  data  still  have  the  accuracy  of  a 
1:250,000  product. 

8. 2. 3.1  Preferred  shift  methods.  Preferred  shift  methods 
and  shift  parameters  may  change  as  additional  information  becomes 
available.  Shift  parameters  and  identification  information  should 
be  treated  as  data  apart  from  the  algorithm  code  so  that  they  may 
be  updated  without  modifying  and  recompiling  the  code.  Some 
examples  of  implementation  tradeoffs  are: 

a.  If  Cartesian  coordinate  are  already  available,  AX,  AY,  AZ 
shifts  can  be  applied  directly  (see  Sections  6.3.1  and 

6.3.2)  rather  than  through  the  Molodensky  approximations 
(see  Section  6.3.3). 

b.  If  the  program  contains  algorithms  for  converting  between 
Cartesian  and  geographic  coordinates  (see  Section  4.3.4), 
the  three-step  method  ( see  Sections  6 . 3 . 1 . 1  and  6 . 3 . 2 . 1 ) 
can  be  used  instead  of  Molodensky  approximations  (see 

Section  6.3.3)  to  apply  AX,  AY,  AZ  shifts. 

c.  Both  the  Molodensky  approximations  (Section  6.3.3)  and 
Cartesian-to-geographic  conversion  equations  (Section 

4. 3. 4. 2)  become  singular  at  the  poles.  Selection  of 
one  technique  or  the  other  will  reduce  the  logic  needed 
to  guard  against  singularity.  An  alternate  equation  to 
compute  geodetic  height  from  Cartesian  coordinates  in 
polar  regions  is  given  as  equation  4.12. 


8. 2. 3. 2  Software  testing.  When  developing  software  based  on 
methods  discussed  in  this  Technical  Report,  it  is  advisable  to 
check  the  software  against  test  points.  Test  points  can  be  found 
in  Section  8.5. 

a.  Issues  relating  to  software  development  and  system 
performance  might  force  a  software  developer  to  use  a 
special-purpose  transformation  technique  that  is  not 
mentioned  in  this  Technical  Report.  The  results  of  the' 
special-purpose  software  should  be  compared  to  those  of 
the  techniques  in  this  Technical  Report  to  insure  that 
the  differences  are  acceptable. 

b.  Another  useful  debugging  aid  is  circular  testing.  This 
is  done  by  using  the  output  from  a  transformation  as  the 
input  for  the  inverse  function;  in  other  words  trying  to 
verify  the  original  numbers  are  returned  within 
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acceptable  accuracy.  Since  the  code  for  a  transformation 
and  its  inverse  might  contain  mistakes  that  cancel  each 
other,  circular  testing  can  not  prove  that  code  is 
correct.  At  best,  it  can  demonstrate  that  code  is 
faulty. 

c.  The  following  is  a  list  of  points  that  might  prove 

troublesome  to  coordinate  transformation  software.  Test 
the  code  at  these  points,  if  it  is  appropriate: 

Poles 

Polar  regions 

Discontinuity  in  X  (transition  between  0,360  degrees  or 
-180,  180  degrees) 

Origin  expressed  in  ((|),  A,) 

<j)  =  0 

X  =  o 

Both  hemispheres  in  <|)  (North  and  South) 

Both  hemispheres  in  X  (East  and  West) 
x  =  0 
y  =  0 

All  four  x-y  quadrants 

Discontinuities  in  x  and  y  at  zone  boundaries,  equator, 
etc. 

8.3  Error  analysis.  Error  analysis  is  an  important  part  of 
the  system  design  process.  Error  analysis  will  estimate  the 
accuracy  of  the  system’s  output  products.  This  is  a  prediction 
of  how  well  the  system  will  meet  its  accuracy  requirements.  The 
system  analyst  can  use  error  analysis  to  perform  system 
implementation  studies.  The  major  steps  in  error  analysis  are: 

a.  Develop  the  system  model. 

h.  Determine  the  estimated  errors  for  each  source  in  the 
model . 

c.  Reformulate  the  estimated  errors  to  a  common  description. 

d.  Propagate  the  errors  through  the  model. 

Each  of  these  topics  is  discussed  below.  An  example  is  given  in 
8.3.5. 


8.3.1  System  model.  The  system  model  simply  is  a 
description  of  the  data  used  in  the  process  and  all  the  steps 
(processes)  used  to  process  the  data.  The  data  include  the  raw 
input  data  and  parameters,  like  datum  shift  parameters,  used  in 
the  processing  of  the  data.  Both  hardware-driven  and  software- 
driven  processes  should  be  included  in  the  model.  For  example, 
errors  in  digitizing  and  plotting  are  primarily  hardware  related; 
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coordinate  conversion  and  datum  shifting  are  primarily  software 
processes.  For  the  datum  shifting  and  coordinate  conversion,  the 
system  model  will  probably  be  a  simple,  serial  process.  That  is, 
the  output  of  one  process  will  be  the  input  for  the  next,  with  no 
feedback  to  prior  processes.  For  trade-off  studies,  alternate 
system  models  can  be  developed  to  evaluate  accuracy  versus 
implementation . 

8.3.2  Error  estimates.  Once  the  system  error  model  has  been 
described,  error  estimates  must  be  assigned  to  each  process  and 
data  source  in  the  model.  When  this  is  done,  the  analyst  will 
quickly  determine  which  error  sources  can  be  neglected  and  which 
must  be  propagated  through  the  model. 

a.  Estimates  of  the  expected  errors  in  the  input  data  may  be 
provided  with  the  data  or  in  product  specifications . 

Also,  the  original  source  or  intended  use  of  the  data  may 
provide  an  indication  of  its  accuracy  (e.g.  digital 
data  derived  from  a  cartographic  source  will  have  an 
accuracy  no  better  than  the  cartographic  source) . 

b.  Without  proper  guidance,  it  is  often  difficult  to 
estimate  the  error  in  the  data.  This  is  especially  true 
for  cartographically  derived  elevations,  because  map 
accuracy  standards  for  elevation  allow  for  error  in  both 
the  horizontal  location  and  the  vertical  direction. 
Therefore,  the  allowable  elevation  error  is  dependent  on 
the  topography  and  may  vary  within  a  map  sheet. 

c.  In  addition  to  providing  datum  shift  parameters. 
Appendices  B  and  E  give  the  estimated  errors  in  the 
parameters  for  AX,  AY,  AZ;  and  seven-parameter 
shifts,  respectively.  (Note:  Correlations  between 
parameters,  which  may  be  significant,  were  not  available. 
Therefore,  they  are  not  given  in  the  appendices.) 

d.  The  accuracy  and  behavior  of  the  equations  presented  in 
this  Technical  Report  are  well  defined.  Table  8.1 
identifies  the  location  of  accuracy  guidance  for  the 
coordinate  conversion,  datum  shift,  and  map  projection 
equations.  Additional  references  describing  equation 
accuracy  and  behavior  are  noted  in  the  text  and  listed  in 
section  2. 

e.  The  errors  introduced  by  system  implementation  factors  • 
such  as  computer  word  size,  accuracy  and  precision  of 
mathematical  functions,  truncation  of  terms,  etc.  must  be 
quantified. 

f.  The  errors  introduced  by  hardware  processes  must  be 
determined  from  equipment  specifications,  analysis, 
and/or  testing. 
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Table  8.1. _ Accuracy  guidance  for  equations 


Sytblect. 


Accuracy  guidance  section (s) 


Coordinate  Conversion  Equations 

Cartesian/geodetic  conversion  4.3.4 

Height  conversion  4.4.3  - 


Datum  shift  Equations 

Seven-parameter  geometric  transformation 
AX , AY , AZ  shifts 

Direct  WGS  72  to  WGS  84  transformation 
Vertical  datum  shifts 


6.2.4 

6.3.1,  6.3.2 
6. 4. 2. 2 

6.7,  6.7.1,  6.7.2 


7.2.3 

7.3.3 

7.4.4 

7.5.4 
7.6.2' 

7.7.4 


Appendix  B 
Appendix  E 
Section  6. 4. 2. 2 

8. 3. 2.1  Error  types.  Depending  on  the  source,  error 
estimates  may  be  described  in  several  ways.  Error  descriptions 
may  have  a  bias  component  and  a  random  (noise)  component.  A  set 
of  error  estimates  describing  the  accuracy  in  multiple  dimensions, 
for  example,  <TAy,  C^,  may  have  correlations  between  the 

parameters.  Generally,  with  MC&G  products,  the  error  descriptions 
do  not  differentiate  between  biased,  correlated,  and  independent 
random  errors;  the  assumption  is  made  that  the  error  estimates  are 
unbiased,  independent  and  random.  The  resultant  error  analysis 
will  be  less  accurate  than  if  the  true  statistical  natures  of  the 
error  sources  were  known. 


Wap  projection  Equations 
Mercator  projection 
Transverse  Mercator  projection 
Universal  Transverse  Mercator  (UTM)  grid 
Lambert  Conformal  Conic  projection 
Polar  Stereographic 
Universal  Polar  Stereographic 

Shift  Parameter  Accuracy  values 
AX , AY , AZ  shifts 

Seven-parameter  geometric  transformation 
Direct  WGS  72  to  WGS  84  transformation 


8.3.3  Reformulation.  Once  the  error  sources  have  been 
quantified,  they  must  be  converted  to  a  common  frame  of  reference. 
The  form  used  for  the  system  output  error  estimate  is  often  a 
convenient  frame  of  reference. 
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a.  The  reference  for  the  error  estimates  must  be  established 
(i.e.  the  ellipsoid  or  the  map  sheet) .  It  may  be 
necessary  to  scale  from  one  reference  to  another. 

b.  The  error  estimates  should  be  converted  to  the  same 
units . 


c.  The  datum  shift  parameter  error  estimates,  G^,  GAy,  GAz, 

given  in  Appendix  B,  can  be  reformulated  to  G^, 

using  equation  set  6.8. 

d.  Error  estimates  typically  are  linear,  circular,  or 
spherical  statistical  values  with  associated 
probability  levels.  Error  estimates  should  be  converted 
to  the  same  type  of  statistic  at  the  same  probability 
level.  Two  independent  linear  standard  deviations 
(68.27%  probability),  Gx  and  Gy,  can  be  converted  to  a 
circular  standard  deviation  (39.35%  probability),  using 


(Gx  +  Gv) 

°c  ~  2 

if  the  values  of  Gx  and  0>y  are  reasonably  close. 
Factors  to  convert  between  different  levels  of 
probability  are  given  in  Table  8.2. 


(8.1) 


8.3.4  Error  propagation.  Assuming  the  error  estimates  are 
statistically  independent  and  have  been  converted  to  the  same 
units  and  statistical  type,  they  may  be  propagated  through  a 
simple,  serial  system  model  using  the  Root -Sum- Square  (RSS) 
technique . 

l 

RSS  =  [  (Nx) 2  +  (N2)2  + -  +  (Nn) 2 ]  2  (8.2) 


If  error  estimates  have  known  biases  or  correlations,  then  a  more 
sophisticated  error  propagation  technique  should  be  used. 
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F 


Error 

Form 

Standard  Error 
(One  Sigma) 

Probable 

Error 

Map  Accuracy 
Standard 

Near  Certainty 
Error 

Linear 

G  (68.27%) 

PE  (50%) 

MAS  (90%) 

30  (99.73%) 

Circular 

Cc  (39.35%) 

CEP  (50%) 

CMAS  (90%) 

3 . 5O0  (99 .78%) 

Spherical 

as  (19.9%) 

SEP  (50%) 

**SAS  (90%) 

40s  (99.89%) 

**SAS  =  Spherical  Accuracy  Standard 

Linear  Error  (One  dimensional) 


68.27% 


99.73% 


50% 


68.27 


90 


99.73 


1.0000 


0.6745 


0.4101 


0.2248 


1.4826 


1.0000 


0.6080 


0.3333 


2.4387 


1.6449 


1.0000 


0.5483 


4.4475 


3.0000 


1.8239 


1.0000 


Circular  Error  (Two  dimensional) 


To 

From 

39.35% 

39.35% 

1.0000 

50 

0.8493 

90 

0.4660 

99.78 

0.2857 

1.0000 


0.5486 


0.3364 


1.8227 


1.0000 


0.6131 


99.78% 


2.9726 


1.6309 


1.0000 


Spherical 

Error  (Three 

Dimensional) 

From 

To 

19.9% 

50% 

90% 

99.89% 

19.9% 


50 


90 


99.89 


1.000 


0.650 


0.400 


0.250 


1.538 


1.000 


0.615 


0.385 


2.500 


1.625 


1.000 


0.625 


2.600 


1.600 


1.000 


(Note:  Errors  are  assumed  to  be  random  and  uncorrelated.) 

8.3.5  Example.  The  commander  wants  a  printed  product 
showing  the  transportation  routes  near  Seoul,  Korea  with  WGS  84, 
UTM  coordinates,  at  a  scale  of  1:50,000.  The  data  source  is  a 
Class  A,  1:50, 000-scale  topographic  map  on  the  Tokyo  datum.  The 
product  will  be  created  by  digitizing  the  transportation  routes  on 
the  map  sheet,  shifting  the  coordinates  from  the  Tokyo  datum  to 
WGS  84  using  the  three-step  method,  and  plotting  the  shifted  data. 


81 


a.  The  system  model  is  shown  in  Figure  8.5.  In  this  model, 
it  is  assumed  that  computer-related  issues  will 

not  degrade  coordinate  conversion  and  datum  shift 
calculations.  Also,  it  is  assumed  that  the 
representation  of  the  digitized  data  will  not  change 
through  the  rest  of  the  process.  Often  this  assumption 
is  not  valid  because  digital  data  representation  must  be 
changed  to  a  predetermined  format.  Such  data 
modifications  can  introduce  error. 

b.  The  error  estimates  for  each  of  the  elements  in  the  model 
are 


Element 


Error 


Map  accuracy  (horizontal) 
Digitizer 

UTM  to  Geodetic  conversion 
Geodetic  to  XYZ  conversion 
Datum  shift  equations 
Datum  shift  parameters: 

°Ay 

®Az 

XYZ  to  Geodetic  conversion 
Geodetic  to  UTM  conversion 
Plotter 


0.02"  CMAS (90%) 
+0.005" 

<1  cm 
0.0 
0.0 
5  m 
3  m 
3  m 

0.1  m 
<1  cm 
+0.01“ 


Only  the  map,  digitizer,  datum  shift  parameter,  and  plotter 
errors  are  significant  and  will  be  considered  in  the  rest  of 
the  analysis. 

c.  The  error  estimates  must  be  reformulated  to  the  same 

units  and  statistical  basis.  Since  the  commander  wants 
to  know  the  estimated  error  on  the  ground,  we  will 
convert  to  circular  error,  in  meters,  on  the  ellipsoid, 
at  50  percent  probability  (CEP) . 

1.  Map  error  must  be  scaled  from  the  map  to  the 
ellipsoid  and  changed  from  CMAS  to  CEP. 

Map  Error  =  0.02"  CMAS  x  0.0254  m/in  x  50,000 
=  25.4  m  CMAS  x  0.5486 
=  13.9  m  CEP 

2 .  It  is  assumed  that  the  digitizer  error  is  a  near 
certainty,  circular  error  specification  that  must  be 
scaled  from  the  map  to  the  ellipsoid  and  converted  to 
CEP. 

Digitizer  error  =  0.005"  max  x  0.0254  m/in  x  50,000 

=  6.4  m  max  x  0.3364 
=  2.1  m  CEP 
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3.  The  datum  shift  parameter  error  estimates,  0^,  0Ay, 
0Az,  must  be  converted  to  0^  and  0^,  then  to  CEP. 

Using  the  above  values  for  0^,  0^Y/  and  0Az,  (j)  =  37.5° 
and  A  =  127°,  equation  set  6.8  yields  0^  =  3.34  m 
and  0^  =  4.38  m. 

CEP  =  1.1774  x  0C  »  1.1774  x  (0^  +  0^) /2*  =  4.5  m 

4.  It  is  assumed  that  the  plotter  error  is  a  near 
certainty,  circular  error  (99.78%  probability) 
specification  that  must  be  scaled  from  the  map  to  the 
ellipsoid  and  converted  to  CEP. 

Plotter  error  =  0.01"  max  x  0.0254  m/in  x  50,000 
=  12.7  m  max  x  0.3364 
=  4.3  m  CEP 

5.  The  RSS  error  for  the  complete  system  model  is 

l 

RSS  =  [(13.9)2  +  (2.1)2  +  (4.5)2  +  (4 . 3 ) 2 ] 2 
=  15.4  m  CEP 

d.  From  the  above  analysis,  it  is  evident  that  the  error  in 
the  source  data  is  the  prime  contributor  to  error  in  the 
output  product  for  this  application.  For  other  types  of 
applications,  such  as  surveying,  this  may  not  be  the 
case. 
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1:50,000 
Class  A  Map 


Transportation 
Route  Map 


8.4  Numerical  Examples. 

8.4.1  Example  8.4.  Convert  WGS  84  geodetic  coordinates 
(<j),  X,,h)WGS  84  to  Universal  Transverse  Mercator  coordinates 

( xutm / Y utm )  the  NAD  27  reference  system  using  a  three-step  AX, 

AY,  AZ  datum  shift. 

This  example  converts  WGS  84  geodetic  coordinates  to 
Universal  Transverse  Mercator  projection  coordinates  in  two  steps. 
First,  the  WGS  84  geodetic  coordinates  are  converted  to  NAD  27 

geodetic  coordinates  using  the  three-step  AX,  AY,  AZ  datum  shift. 
Geodetic  height  is  needed  for  this  conversion.  Next,  the  NAD  27 
geodetic  latitude  and  longitude  are  converted  to  UTM  coordinates. 
Geodetic  height  is  not  needed  for  this  step. 

Beginning  WGS  84  coordinates: 

$  =  42°  56'  52 . 163 "N  =  0.749580918883376  radians 
X  =  288°  22'  24 . 326 "E  =  1.891473118  radians 
h  =  203.380  m 

Step  1:  Convert  WGS  84  geodetic  coordinates  to  NAD  27  geodetic 
coordinates  using  the  three-step  AX,AY,AZ  datum  shift. 

1.1  Convert  (<]), X, h)WGS  84  to  (X,Y,Z)WGS84  using  equation  sets 
4.4  and  4.5. 

a  =  6378137 

f  =  0.00335281066474 

£2  =  2f  -  f2  =  0.00669437999013 

b  =  a(l  -  f)  =  6356752.3142 

Rn  =  ,  . ===  6388070.57383009 

N  i/l  -  £2  sin2  <j) 

XWGS  84  =  (Rn  +  h)  cos<t>  cosA,  =  1473933.5413 
ywgs  84  =  (Rn  +  h)  cos(|>sinA,=  -4437679.0666 
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+  hj  sin(j)  =  4323399.2717 

1.2  Convert  ( X ,  Y ,  Z )  34  to  ( X ,  Y ,  Z )  27  ts  ing  e  qua  t  ion  sot 

6.4  and  AX,AY,AZ  datum  shift  values  from  Table  B.l  NAS-C. 

XNAD  27  =  XWGS  84  -  AX  =  1473933.5413  -  (-8.0)  =  1473941.5413  m 

ynad  27  =  ywgs  84  -  AY  =  -4437679.0666  -  160.0  =  -4437839.0666  m 


ZNAD  27  =  ZWGS  84  ~  Az  =  4323399.2717  -  176  =  4323223.2717  m 

1.3  Convert  (X,Y,Z)NAD27  to  ((|),  27  following  equations 

given  in  Section  4 . 3 . 4 . 2 . 


A-nad  27  =  arctan  2?>)  =  5.033056560575  radians 

V-^NAD  27  ) 

Calculating  <|)  is  an  iterative  process. 

1.3.1  Compute  an  initial  approximation  to  a  variable  P, using 
equation  4.7.  Parameters  a  and  b  for  the  Clarke  1866  ellipsoid 
are  found  in  Table  A. 

Po  =  arctan  [ - ,a  Z  ]  =  0.74788820090716142  radians 

b  VX2  +  Y2 

1.3.2  Substitute  PQ  for  p  in  equation  4.8  to  get  an 
approximation  for  <|). 


x  r  Z  +  8'2  b  sin3  p  n 

O  =  arctan  I  ,  - - I 

>/x2  +  Y2  -  a  £2  cos3  p 


0.74958142853674489  radians 


£2  =  2f  -  f2  =  0.0067686579973 


and 

02 

£'2  =  - - -  =  0.0068147849459 

1  -  e2 


For  most  applications,  this  first  iteration  is  sufficient. 

Step  2:  Convert  the  NAD  27  geodetic  latitude  and  longitude  to 
Universal  Transverse  Mercator  projection  coordinates. 

2.1  Find  the  UTM  zone  z  and  the  central  meridian  XQ  from 
equations  7.22  and  7.23. 
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z  =  greatest  integer  < 
z  =  19 


180  X 


-  29 


X0  =  (6z  +  177)  —  =  5.078908123303  radians 


2.2  Calculate  Transverse  Mercator  projection  coordinates 
(X^Y^)  using  equations  7.8  and  7.9. 

2.2.1  Calculate  terms  in  equation  set  7.9  using  £2  and  8'2  as 
calculated  in  Step  1.3.2. 

*•  ■  1  -  ¥  -  h'  -  Ji86  -  lift?8  -  ”8305681856 


A  -  1 
A2  -  8 


A  -  -15- 

4  "  256 


6  (e‘  + 


ig,  +  .  -i55^8 

4T  128  4096fc 


')’  0 


.2542555420215  x  10'2 


64  +  t£6  -  =  0.2698010542695  x  10'5 


A  - 

"  3072 


/  41  l  \ 

£6  -  J^E8  =  0.3502448582027  x  10“8 


A8  =  -  J31Q72  E8  =  -0.5044416191579  x  lO'*1 


S(|)  =  a  [ A0<|>  -  A2  sin  2<|)  +  A4  sin  4<j)  -  A6  sin  6(|)  +  A8  sin  8<})] 
S(j)=  4756711.680282 


N  =  RN  =  .  =  6388250.562973 

yi  -  £2  sin2  <|) 

Note  that  N  =  Rn  differs  from  RN  calculated  in  Step  1.1,  as  Step  2 
is  in  the  NAD  27  reference  frame  and  Step  1.1  is  in  the  WGS  84 
reference  frame. 


T[  =  8'  cos  <])  =  0.6042571576496  x  10'1 
t  =  tan(J)  =  tan(0. 74958094822852639)  =  0.9308149267480 
A  =  X-X  =  -0.04585156272818  radians 

0 
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2.2.2  Calculate  the  Transverse  Mercator  projection 
coordinates  (x  ,y  )  . 

TM  J  TM 


xTM  =  N  A  cos<J) 

N  A5cos5(|) 
+  120 

N  A2 

Ytm  -  +  2 


N  A3  cos3  d)  ... 

+  - 6 -  (1  "  t2 

(5  -  18t2  +  t4  +  14T|2 
sin  (j)  cos  (|) 


+  T|2) 

-  58t2T[2) 


N  A4 

+  ~~24~  S3-n  0  cos3  4*  ^  “  t2  +  9T|2  +  4T|4) 

+  Nrj^  sin  (|)cos5  (|)  (61  -  58t2  +  t4  +  270T12  -  330t2Tl2) 

XrM  =  -214408.9715411 
yTM  =  4760061.989246 

2.3  Convert  the  Transverse  Mercator  projection  coordinates 
(Xt^Ytm)  to  Universal  Transverse  Mercator  projection  coordinates 
(XuTMrYuTM)  US3Ln9  equation  set  7.24  (Northern  hemisphere). 


xUTM  =  0.9996  X^  +  500,000  =  0.9996  (-214408.9715411)  +  500,000 
yUTM  =  0.9996  y^  =  0.9996  (4760061.989246) 


The  resulting  Universal.  Transverse  Mercator  coordinates  (xUTM,yUTM) 
in  the  NAD  27  reference  system  are: 

Xutm  —  285676.792 


Yutm  =  4758157.964 

8.4.2  Example  8.5.  Convert  WGS  84  geodetic  coordinates 
(<j), A,,h)  to  geodetic  coordinates  ((j),A„h)  in  the  NAD  27  reference 
system  using  the  standard  Molodensky  method. 

This  example  converts  WGS  84  geodetic  coordinates  to  NAD  27 
coordinates  using  the  standard  Molodensky  method. 

Beginning  WGS  84  coordinates: 

<))  =  42°  56  1  52 . 163 "  N  =  0.749580918883376  radians 
X  =  108°  22' 24. 326"  W  =  1.891473118  radians 
h  =  203.380 
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Step  1:  Compute  Aa,  Af,  W2,  M,  N,  £2,  and  £'2  using  equation  set 
6.5.  Obtain  aWGS  94*  Anad  27'  ^wgs  84'  and  fnAD27  from  Table  A.  1. 

Aa  =  aWGS  84  -  aCLARKE  1866  =  "69.4 

Af  =  f^Qs  84  -  f CLARKE  1866  =  3.7  2  6  4  6  3  8  6  3  X  10'3 

£2  =  2f  -  f2  =  0.006694379990 
£2 

£’ 2  =  — 7  =  0.006739496742 
1-E2 

W  =  Vl  '  e2  sin2<|)  =  0.9984449805 
N  =  RN  =  “  =  6388070.574 


M  = 


a(l  -  £2) 
W3 


6365086.681 


Step  2:  Compute  A(|)  and  AX,  and  Ah  using  equation  set  6.5,  and 
AX,  AY,  and  AZ  datum  shifts  from  Appendix  B  for  datum  code  NAS-A. 

In  this  example,  the  signs  of  AX,  AY,  and  AZ  are  reversed  from 
those  of  Appendix  B  due  to  the  direction  of  the  datum  shift  (WGS 
84  to  local  coordinates) . 

A<|)  =  [  (-sin(])cosXAx  -  sin<j)sinAAy  +  cos(j)Az  +  ^  si^cos<^a  + 

sin<|>cos<|) (2N  +  E,2Msin2<|))  (1  -  f)Af]/[M  +  h] 


A(j)  =  0.131" 


aX 

Ah 


[-sinXAx  +  cosXAy] 

— rr^ — r; — rr~  =  2.615" 

{ (N  +  h) cos(j)] 

cos(])cosAAx  +  cos(|)sinXAy  +  sin(|)Az 


a{l  -  f )  .  „ 

WAa  +  - — -  sm2f  Af 

w 


Ah  =  28.65  m 


Step  3 :  Compute 
equation  set  6.7 


the  NAD  27  geodetic  coordinates  ((|),A,,h) 
that  results  in  the  following  valuesl; 


from 
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^NAD  27  _  ^WGS  84  +  “  ^2°  56' 52.294"N 

^NAD  27  =  ^WGS  84  +  ^  =  108°  22 '21. 711  "W 
h-NAD  27  =  ^WGS  84  +  =  232.03  Itl 


8.4.3  Example  8.6.  Convert  WGS  84  geodetic  coordinates  (([>, 
A,)WGS  84  to  Mercator  projection  coordinates  (x,y)  . 

This  example  converts  WGS  84  geodetic  coordinates  to  Mercator 
projection  coordinates  within  the  WGS  84  reference  system. 

Geodetic  height  (h)  is  not  used  in  this  conversion.  The  Mercator 
projection  central  meridian  XQ  is  288°  =  5.026548246  radians 

Beginning  WGS  84  coordinates  are 

<1>wgs  84  =  42°  56  •  52. 163  "N  =  0.7495809189  radians 

Kgs  84  =  288°  22'  24.326"E  =  5.033057346  radians 

Step  1:  Use  equation  7.1  and  WGS  84  datum  parameters  a  and  f 
from  Table  A.l  to  calculate  x. 

X  =  a{X  -  K)  =  6378137(5.033057346  -  5.026548246)  =  41569.357 


Step  2:  Calculate  £  from  f  using  equation  4.9,  then  use  equation 
7.2  to  calculate  y. 

8  =  *s[2f  -  f2  =  0.0818191908426 


y  =  a  In 


p 

i  A  jh 

(  1  -  £  sin  0  ^ 

e" 

2 

tan 

i  +  2)  ' 

(^1  +  £  sin  (j)^ 

=  5274911.868 


The  resulting  Mercator  projection  coordinates  (x,y)  are 


x  =  41569.357 


y  =  5274911.868 

8.4.4  Example  8.7.  Converst  Mercator  projection  coordinates 
(x,y)  in  the  WGS  84  reference  system  to  WGS  84  geodetic  latitude 
and  longitude  (Cj), X,)  . 
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This  example  converts  WGS  84  Mercator  coordinates  to  WGS  84 
latitude  and  longitude. 


Beginning  Mercator  coordinates  and  central  meridian  (XQ) : 

X  =  41569.3572 
y  =  5274911.8684 

A,0  =  288°E  =  5.0265482457  radians 

Step  1:  Compute  the  longitude  {X)  using  equation  7.4  where  a  is 
the  semi -major  axis  of  the  WGS  84  ellipsoid  listed  in  Table  A.l. 

X  =  ~  +  A,0  =  5.0330657221  radians  =  288°  22'24.326,,E  or 

Cl 

71°  37'  35 . 674 "W 

Step  2:  Compute  an  intermediate  value  t  from  equation  7.5. 

„  5274911.8684 

t  =  e<-  P  =  e("  6378137  ]  =  .43734619376064 


Step  3:  Compute  an  initial  value  of  the  geodetic  latitude  ((j)0) 
from  equation  7.6. 

•  7C 

<Po  =  2  -  2  arctan(t)  =  0.746233651  radians 


Step  4:  Iterate  equation  7.7  until  the  updated  value  of  (j)  is 
sufficiently  close  to  the  previous  value  of  (J). 

£2  =  2  f  -  f 2  =  0.00669437999013 


8  =  0.0818191908426 
^n+l  =  f  -  2  arctan 


1  -  £  sin  <J)n 

e  \ 
2 

t 

V 

1  +  £  sin 

J 

(j)!  =  0.7495688569418636 


The  geodetic  latitude  will  be  given  by  the  last  value  of  <])  • 

calculated. 

Final  §  =  0.7495809188833704  radians  =  42°  56'  52.163 "N 

The  resulting  WGS  84  geodetic  latitude  and  longitude  ((jbA,)  are: 

<j)=  42°  56  ■  52. 163  "N 
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X  =  288°  22'  24 . 326  "E  or  71°  37'  35.674"W 


8.4.5  Example  8.8.  Convert  Universal  Transverse  Mercator 
coordinates  (Xutm/Yutm)  in  the  NAD  27  reference  system  to  NAD  27 
geodetic  latitude  and  longitude (<|),  A,). 

This  example  converts  NAD  27  UTM  coordinates  to  NAD  27 
geodetic  latitude  and  longitude.  The  UTM  coordinates  are 
converted  to  Transverse  Mercator  projection  coordinates.  The 
Transverse  Mercator  projection  coordinates  are  then  converted  to 
NAD  27  latitude  and  longitude. 

Beginning  coordinates  are 
Xutm  =  285,677.332 
Y  utm  =  4,758/154.856 
Zone  =  19 


Step  1:  Convert  Universal  Transverse  Mercator  coordinates 
(Xutm/Yutm)  to  Transverse  Mercator  coordinates  (xTM,yTM)  using 
equation  set  7.25. 

Xtm  =  q  999^  (^utm  -  500,000)  —  -214,408.4314 


Ytm 


0  9996  ^ mm  =  4,760,058.8796 


Step  2:  Convert  Transverse  Mercator  coordinates  (xTM,yTM)  to 
geodetic  latitude  and  longitude  (<|),  X)  . 


2.1  Calculate  the  UTM  zone  central  meridian  (X0)  using 
equation  7.23. 


XQ 


(6z  +  177) 


5.078908123303  radians 


2.2  Compute  the  eccentricity  (£)  and  second  eccentricity  (£') 
of  the  Clarke  1866  ellipsoid,  the  reference  ellipsoid  for  NAD  27, 
using  equation  set  7.9. 

£2  =  2f  -  f2  =  0.6768657997291  X  10'2 


£2 

£12  =  - - -  =  0.6814784945915  X  10’2 

1  -  £2 

2.3  Calculate  the  subscripted-A  terms  of  equation  set  7.9, 
using  £2  calculated  in  Step  2.1.  The  subscripted-A  terms  are  used 
in  the  calculation  of  the  footpoint  latitude. 
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13  5  175 

A°  =  1  '  4  62  -  64  8‘  ‘  256  ®6  ’  16384  8®  =  0-998305681856 


A,  =  T 


3  /  1  15  455  \ 

g  fe2  +  j  e4  +  Y28  £6  “  4^  E8J  =  0.2542555420215  X  10'2 
^ r  ^E4  +  |  E6  -  68^  =  0.2698010542695  X  10-5 


Aa  = 


As  = 


3072 


f86-#88 


3502448582027  X  10'8 


Ar  = 


315 

■131072  e8  =  -  0.5044416191579  X  lO’11 


2.4  Compute  the  footprint  latitude,  <J>i,  using  equation  7.10 
and  a  from  Table  A.l. 

Make  an  initial  estimate  for  0i: 

0i  =  =  0.7463005398433  radians 

a 


2.5  Calculate  using  equation  set  7.9  and  the  estimate  of 


0,  0i. 

=  a[Ao0  -  A2  sin  20  +  A4  sin  40  -  A6  sin  60  +  A8  sin  80] 

S4>  =  0.4735829138159  X  107 

2.6  Calculate  S  ’ 4,  using  equation  set  7.11  and  the  estimate  of 

0,  0i. 

S  1  (j,  =  a[A0  -  2A2cos  20  +  4A4  cos  40  -  6Ag  cos  60  +  8A8  cos  80] 

S  ’ «!>  =  0.6364798134276  X  107 

2.7  Find  an  updated  value  of  the  footpoint  latitude,  0r  using 
equation  7.10. 

0i  (updated)  =  0i  (current)  ~  g  , ~ —  =  0.7501073756426  radians 

2.8  Iterate  Steps  2.5  through  2.7  until  successive  values  of 
0i  are  as  close  as  is  desired. 
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For  this  example,  the  final  value  of  (|>i  is: 
Final  (]>!  =  0.7501073020563  radians 


2.9  Calculate  values  used  in  equation  set  7.12  using 
equation  set  7.13. 

tj.  =  tan  <j>!  =  tan  (0.7501073020563) 

Ni  =  5  6,388,261.949762 

Vl  -  £2  sin2  <j)i 


Ri  =  3  =  6,365,044.225886 

\yl  -  £2  sin2  <()i  ) 


Tlx2  =  £’ 2  cos2  <)>!  =  0.3647692467262  X  10'2 

B3  =  1  +  2ti2  +  Tlx2  =  2.740138642336 

B4  =  5  +  3tx2  +  T|i2  -  4T|i4  -  9tx2T|i2  =  7.579827062321 

B5  =  5  +  28tx2  +  24tx4  +  6TU2  +  8ti2T|i2  =  47.45050110665 

B6  =  61  +  90tx2  +  46T|x2  +  45tx4  -  252ti2TU2  -  3Tlx4  -  - 

90tx4Tli2  +  225tx4T|i4  =  172.1890097418 

2.10  Calculate  the  geodetic  latitude  and  longitude  (<{>,A,) 
using  equation  set  7.12. 


x  x  tiNi  Tl  /x\2  B4  rxy 

*  ■  -  ~ft  u  y  -  si  y 

<j>  =  0.74958094314009481  radians 


+ 


42°  56.  52.168" 


^0  +  sec 

X  =  5.03305669671  radians =  288°  22'  22.464“ 

8.4.6  Example  8.9.  Convert  WGS  84  latitude  and  longitude 
(<{>,A,)  to  Lambert  Conformal  coordinates  (x,y)  using  two  standard 
parallels . 

This  example  converts  WGS  84  latitude  and  longitude  to  Lambert 
conformal  projection  coordinates  in  the  WGS  84  reference  system. 


94 


Beginning  WGS  84  coordinates  are 

<1>WGS  84  =  42°  56-  52.163 "N  =  0.749580918883376  radians 
A«gs  84  =  288°  22-  24.326 "E  =  5.033065722110382  radians 


Standard  parallels: 

<J>!  =  42°  30-  00.000-  =  0.7417649320975901  radians 

<j)2  =  43°  00-  00.000-  =  0.7504915783575617  radians 
Origin: 

<J>0=  42°  30-  00.000-  =  0.7417649320975901  radians 

Xq=  288°  00-  00.000-  =  5.026548245743669  radians 

Step  1:  Calculate  eccentricity  £  from  flattening  f. 

£  =  V2f  -  f2  =  0.0818191908426 

Step  2:  Calculate  q  at  <]),  <|)i,<|)2and  <j>o  using  equation  set  7.27. 


q(0o)  =  0.8164293462094778 

Step  3:  Calculate  N  at  <}>i  and  fusing  equation  set  7.27  and  a 
from  Table  A-l. 
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N«t>j)  =  6387903.46785173 
N(02)  =  6388090.05758619 


Step  4 :  Calculate  L  from  equation  set  7 .27 . 
ln(Nicos0i)  -  In  (N9COS02) 

L  =  - - - - —  =  0.6788029306559900 

02  -  qi 

Step  5:  Calculate  K  from  equation  set  7.27. 

N,  cos  0.  Nn  cos  0n 

K  =  t  ~ .  T1  =  I  ~  .  =  12076169.25203636 

l  e-1^  l  e*^ 

Step  6:  Calculate  r  and  rG  from  equation  set  7.27. 
r  =  KB'1^)  =  6888432.18075871 


r0  =  K£-L£i<V  =  6938179.32109766 

Step  7:  Calculate  the  Lambert  coordinates  (x,y)  using  equation 
set  7.26. 

x  =  r  sin  L(X-  X0)  =  30474.8898082 
y  =  r0  -  r  cos  L (A.  -  Xq)  =  49814.5521555 

8.4.7  Example  8.10.  Convery  Lambert  Conformal  projection 
(with  two  standard  parallels)  coordinates  (x,y)  in  the  WGS  84 
reference  system  to  WGS  84  geodetic  latitude  and  longitude  (0,A,)  . 

This  example  converts  WGS  84  Lambert  Conformal  coordinates  to 
WGS  84  latitude  and  longitude. 

Beginning  Lambert  Conformal  coordinates,  standard  parallels 
(0i,02)  and  origin  (0O'^o)  are 

x  =  30474.890 

y  =  49814.552 

0X  =  42°  30 1 00 "  =  0.7417649320976  radians 

02  =  43°  OO-OO"  =  0.7504915783576  radians 

0O  =  42°  30'00“  =  0.7417649320976  radians 

X0  =  288°  001 00“  =  5.0265482457437  radians 
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Step  1:  Compute  £,  L,r0,k,N1,N2,q0,q1,  and  q2  using  equation  set 
7.27.  Obtain  a  from  Table  A.l. 


£  =  \2f  -  f 2  =  0.0818191908426 


In 

tan  ff  ♦  *f)  • 

fl  -  £  sin  (J)^ 

e 

2 

\4  2  J 

3-  +  £  Sin  <l>i; 

<3o  =  0.8164293462094778 
<2i  =  0.8164293462094778 

<h  =  0.8282703592274435 
a 


N<  = 


yi  -  £2 


smJ 


Nj.  =  6387903.467851727 
N2  =  6388090.057586191 


L  = 


K  = 


In  (Nxcos^!)  -In  (N2cos({)2) 


Qk  -Qi 

NlCOSCl)!  N2COS(})2 


=  0.6788029306559713 


=  12076169.25203651 


Le-^  Le-Lti2 

Step  2 :  Compute  rG  using  equation  set  7.27. 

rD  =  Ke-^<V  =  6938179.321097849 
Step  3:  Compute  r  and  0  from  equation  set  7.29. 


0  =  arctan 


x 


ro  -  y 


=  4.4240820859501107  x  10'3  radians 


r  =  — — jf-  =  6888432.18091526 
COS0 

Step  4:  Compute  geodetic  longitude  (A,)  using  equation  set  7.28. 

Q 

X  =  -  +  XQ  =  5.0330657221511255  radians  =  288°  22'  24.326" 


97 


Step  5:  Compute  q  using  equation  set  7.28. 

/k  \ 

In  (“) 

q  =  - 7 -  =  0.8270301920802524 


Step  6:  Find  an  initial  approximation  ((j)')  for  the  geodetic 
latitude  from  equation  set  7.30. 

<t>'  =  2  arctan  (e^)  -  j  =  0.746233650926621 


Step  7:  Compute  f  (<|)  ’ )  and  f'«|>')  from  equation  set  7.30. 


f(<j)')  = 
£(<))■)  = 

f  •  (<|)  • )  = 


fl+sin<])’  ^ 

( 1  -  £  sin  <1> '  ^  e 

_^l-sin(j) '  J 

JL  +  £  sin  <J)'J 

-4.5493507336696604  x  10'3 
1  -  £2 


(1  -  £2  sin2  <J)')cos<f)’ 


=  1.357001593456580 


Step  8:  Find  an  updated  approximation  ((j)")  for  the  geodetic 
latitude  using  equation  7.31. 


0" 


f«D-) 

f  •  ((j); ) 


0.7495861530575 


Step  9:  Iterate  steps  7  and  8  until  successive  values  of  (J)"  are 
sufficiently  close  after  many  iterations. 


Final  «|>)  =  0.74958091885881052  radians  =  42°  56*  52.163“. 
The  resulting  WGS  84  coordinates  are: 


<|>  =  42°  56'  52.163" 
X  =  288°  22-  24.326" 


8.4.8  Example  8.11.  Convert  geodetic  latitude  and  longitude 
(<j),A,)  on  the  International  Ellipsoid  to  Universal  Polar 
Stereographic  coordinates  (x,y) . 

This  example  converts  geodetic  latitude  and  longitude  to  UPS 
projection  coordinates.  Coordinates  are  computed  on  the 
International  Ellipsoid. 
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Beginning  coordinates  are 

<|>  =  87°  17*  14.400*'  S  =  -1.523451361952464  radians 
X  =  132°  14*  52.303“  E  =  2.308160619653466  radians 


Step  1:  Compute  the  eccentricity  (8)  of  the  International 
Ellipsoid  using  equation  set  7.38  and  a  and  b  from  Table  A.l. 


82  =  2f  -  f2  =  0.0067226700623 
8  =  0.08199189022 

Step  2:  Compute  k  using  equation  set  7.38. 

e 

0^2  c  \  — 

K  =  -r—  [  "  ] 2  =  12713920.17558161 

b  u  (1  +  8) 


Step  3:  Compute  r  using  equation  set  7.38. 


r 


Ktan 


(n  i  <b  i  > 

(l  +  8  sin  101^ 

U  '  2  J 

^1  -  8  sin  1  <j>  1 ; 

303059.087879879 


Step  4:  Compute  UPS  coordinates  using  equations  7.49  and  7.50. 

x  =  2,000,000  +  0.994  r  sin  X  =  2222991.410 

y  =  2,000,000  -  0.994  r  cos  X  =  1797464.051 

8.4.9  Example  8.12.  Convert  Universal  Polar  Stereographic 
coordinates  (X,Y)  to  geodetic  latitude  and  longitude  (([>,A,). 


This  example  converts  UPS  coordinates  calculated  on  the 
International  Ellipsoid  to  geodetic  latitude  and  longitude  on  the 
International  Ellipsoid. 

Beginning  UPS  coordinates  in  the  Southern  hemisphere  are 


x  =  2222991.410 
y  =  1797464.051 

Step  1:  Convert  Universal  Polar  Stereographic  coordinates 
(XUPS,YUPS)  to  Polar  Stereographic  coordinates  (XPS,YPS)  . 
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Xups  -  2000000 

Xps  =  ~  0  994 -  =  224337. 434205231 

Yups  -  2000000 

yps  =  ~  — 5—  9~9^ -  =  -203758.500100604 

Step  2:  Compute  the  longitude  (A)  from  equation  7.39. 

A,  =  K  +  arctan 

Yps 

A  =  2.308160619653466  radians  =  132°  14' 52. 303"  E 

Step  3:  Compute  k  using  equation  7.41,  and  a  and  b  from  Table 
A.  1 . 

„  e 

2a2  r  (1  -  £)  1  o 

K  =  ~b~  I  (1  Ve)  •*  =  12713920. 17558161 

Step  4:  Compute  q  using  equation  set  7.40. 

e.q  _  — DLLes - 

“  K  I  cos  A| 

q  =  3.73652495299514 

Step  5:  Compute  an  initial  estimate  (<j)0)  of  the  geodetic  latitude 
using  equation  set  7.42. 

x  It 

<j)0  =  2  arctan  =  1.52313176628106  radians 

Step  6:  Compute  f((|>0)  and  f '  (<|)0)  using  equation  set  7.42. 


f((J)0)  =  -0.006730109061328 


1  -  £2 


f '  (<*>)  ~  (1  -  E2  sin2(])) 


COS(J) 


=  20.9875717033244 


Step  7:  Compute  an  updated  geodetic  latitude  (j>lf  using  equation 
set  7.43. 

d)n  =  0n-i  "  =  1.523452437445710 

Yn  Yn  1 
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Step  8:  Iterate  steps  6  and  7  until  successive  values  of  the 
geodetic  latitude  are  sufficiently  close. 

Final  (<j))  =  1.523451361952464  radians  =  87°  17'  14.4"S. 

The  resulting  geodetic  latitude  and  longitude  ((]),A,)  coordinates 
are: 

0  =  87°  17'  14. 4  "S 
X  =  132°  14 1 52 . 303  "E 

8.5  List  of  test  points.  The  test  points  listed  here  are 
taken  from  DMA  TM  8358.2  (1989)  and  Snyder  (1987) . 

Mercator 

Central  Meridian  =180° 

Ellipsoid  =  Clarke  1866 

<|)  =  35°  N 

X  =  75°  W  (285°  E) 

x  =  11,688,673.7  m 
y  =  4,139,145.6  m 

k  =  1.2194146 

Universal  Transverse  Mercator 

Central  Meridian  =  285°  E  (75°  W) 

Ellipsoid  =  Clarke  1866 

<)),  =  40°  30'  N 

X  =  73°  30'  W  (286°  30'  E) 

x  =  627,106.5  m 

y  =  4,484,124.4  m 

k  =  0.9997989 

Lambert  Conformal  Conic 

Ellipsoid  =  Clarke  1866 

Standard  parallels  =  33°  N 

=  45°  N 

Origin:  ((>0  =  23°  N 

?l0  =  96°  W  (264°  E) 
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>1  X 


0  =  35°  N 
X  =  75°  W  (285°  E) 

x  =  1,894,410.9  m 

y  =  1,564,649.5  m 

k  =  0.9970171 

Universal  Polar  Stereographic 

Ellipsoid  =  International 

<|)'  =  87°  17'  14 . 400  "S 

X  =  132°  14'  52. 303 "E 

=  1,797,464.051  m 
=  2,222,991.410  m 
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9 .  NOTES 


9.1  Intended  use.  This  Technical  Report  provides  guidance 
on  coordinate  conversion  and  datum  shifting  for  a  wide  range  of 
Army  mapping,  charting,  and  positioning  applications.  Guidance  is 
provided  to  aid  the  Army  and  others  in  developing  software  for 
applications  of  datum  shifts  and  other  transformations  and  in 
implementing  this  software  in  support  of  operational  units. 

9 .2  Subject  (kev  word)  listing 

Cartography 

Charts 

Coordinate  Systems 

Datums 

Elevation 

Geodetic  Satellites 
Geodesy 

Geodetic  Coordinates 
Geographic  Coordinates 
Grids 

Hydrographic  Surveying 
Latitude 
.  Longitude 
Map  Projection 
Mapping 
Maps 
MC&G 

Mean  Sea  Level 

Military  Standards 

Photogrammetry 

Projections 

Surveying 

Topographic  Maps 

Topography 

UTM 
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APPENDIX  A 


REFERENCE  ELLIPSOID  PARAMETERS 


10.  GENERAL 

10.1  Scope.  This  appendix  identifies  reference  ellipsoids  with 
associated  geometric  parameters. 

20.  APPLICABLE  DOCUMENTS. 

The  ellipsoids  and  associated  parameters  are  taken  from  DMA 
TR  8350.2  (Edition  2,  Insert  1.)  Note:  Use  semi -major  axis  a  and 
flattening  f  to  calculate  E  and  E' .  Do  not  use  semi-minor  axis  b, 
which  is  a  derived  geometric  constant. 

30.  GENERAL  REQUIREMENTS. 

30.1  Application.  The  ellipsoid  furnishes  a  simple,  consistent, 
and  uniform  reference  system  for  geodetic  coordinates. 


Preceding  Page  Blank 
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TABLE  A.  1 

REFERENCE  ELLIPSOIDS 
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APPENDIX  B 


AX,  AY,  and  AZ  DATUM  SHIFT  PARAMETERS 

10 .  GENERAL 

10.1  Scope.  This  appendix  identifies  the  ellipsoid  center  shift 
transformation  parameters  used  to  transform  a  local  geodetic  system 
to  WGS  84.  An  identifier,  the  reference  ellipsoid,  AX,  AY,  AZ 

shifts,  and  the  estimated  errors  (G)  in  the  shifts'  are  listed  for 
each  local  geodetic  system. 

20.  APPLICABLE  DOCUMENTS. 

The  datum  shift  parameters  and  associated  information  for  each 
local  geodetic  system  are  taken  from  DMA  TR  8350.2  (Second  Edition, 
Insert  1. ) 

30.  GENERAL  REQUIREMENTS 

30.1  Application.  The  AX,  AY,  and  AZ  shifts  are  used  for  datum 
transformations  using  the  three-step  method  discussed  in  section 
6.3,  and  shown  by  equation  sets  6.3  and  6.4. 


The  following  symbology  is  used  in  Table  B.l. 

*  WGS  84  minus  local  geodetic  system. 

**  Derived  from  non-Dobbler  sources,  accuracy  unknown. 

***  Derived  from  a  single  Doppler  station,  accuracy -not  verified. 

+  See  Appendix  F  to  transform  Old  Hawaiian  Datum  coordinates  on 
the  International  Ellipsoid. 

++  Also  known  as  Hito  XVIII,  1963. 

+++  Use  A  =  6,377,483.865  meters  for  the  Bessel  1941  Wllipsoid  in 
Namibia . 

&  Derived  from  non-satellite  sources,  accuracy  unknown. 
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APPENDIX  C 


DATUM  LIST 

10.  GENERAL 

10.1  Scope .  This  appendix  consists  of  a  list  of  countries  and 
geographic  areas  and  their  associated  datums. 

20.  APPLICABLE  DOCUMENTS. 

This  section  is  not  applicable  to  this  appendix. 

30.  GENERAL  REQUIREMENTS 

30.1  Application.  The  intended  use  of  this  appendix  is  to  establish 
table  of  countries  or  geographic  areas  for  quick  reference  to  their 
respective  datum. 


TABLE  C-l 

Countries  and  Their  Associated  Datums 


~ NTRY / GEOGRAP H I C  AREA 


Afganistan 


Alaska 


Alberta 


Antigua 


Argentina 

Ascension  Island 

Australia 

Austria 

Bahama  Islands 

Bahrain  Island 

Baltra,  Galapagos 
Islands 

Bangladesh 

Barbados 

Barbuda 


DLEZM 

Herat  North 

North  American  1927,  1983 

North  American  1927,  1983 

Antigua  Island  Astro  1943 
North  American  1927 

Campo  Inchauspe, 

South  American  1969 

Ascension  Island  1958 

Australian  Geodetic  1966, 

European  1950,  1979 

Cape  Canaveral, 

North  American  1927 
Ain  el  ABD  1970 

South  American  1969 

Indian 

North  American  1927 

North  American  1927 


1984 


TABLE  C-l  (cont'd) 

Countries  and  Their  Associated  Datums 


7NTRY  L GEQGPAP H I C . .  AREA 


datum 


Belgium 

Belize 

Bolivia 

Botswana 


European  1950 

North  American  1927 

Provisional  South  American  1956 
South  American  1969 
ARC  1950 


Brazil 

British  Columbia 


Corrego  Alegre, 

South  American  1969 
North  American  1927,  1983 


Brunei 


Timbalai  1948 


Burkina  Faso 


Adindan 
Point  58 


Burundi 

Caicos  Islands 

Cameroon 

Canada 

Canal  Zone 


ARC  1950 

North  American  1927 

Adindan 

Minna 

North  American  1927,  1983 

North  American  1927 


Canary  Islands 


Pico  de  Las  Nieves 


TABLE  C-l  (cont'd) 

Countries  and  Their  Associated  Datums 


DATUM 

Caribbean 

North  American  1927 

Caroline  Islands 

Kusaie  Astro  1951 

Cayman  Brae  Island 

L.C.  5  Astro  1961 

Central  America 

North  American  1927,  1983 

Chatham  Island 

Chatham  Island  Astro  1971 

Chile 

Provisional  South  American  1956, 

South  American  1969,  PS  Chile  1963 

Cocos  Island 

Anna  1  Astro  1965 

Colombia 

Bogota  Observatory, 

South  American  1969 

Congo 

Point  Noire  1948 

CONUS 

North  American  1927,  1983 

Corvo  Island  (Azores) 

Observatorio  Meteorologico  1939 

Costa  Rica 

North  American  1927 

Cuba 

North  American  1927 

Cyprus 

European  1950 

TABLE  C-l  (cont’d) 

Countries  and  Their  Associated  Datums 


COUNTRY /GEOGRAPHIC  AREA 


DATUM 


Denmark 
Diego  Garcia 
Djibouti 

Dominican  Republic 
East  Canada 
East  Falkland  Island 
Easter  Island 
Eastern  United  States 
Ecuador 

Eft ate  Island 
Egypt 

El  Salvador 
England 


European  1950 

ISTS  073  Astro  1969 

Aybella  Lighthouse 

North  American  1927,  1983 

North  American  1927,  1983 

Sapper  Hill  1943 

Easter  Island  1967 

North  American  1927,  1983 

Provisional  South  American  1956 
South  American  1969 

Bellevue  (IGN) 

European  1950 
Old  Egyptian 

North  American  1927,  1983 

European  1950,  Ordnance  Survey 
of  Great  Britain  1936 


Erromango  Island 


Bellevue  (IGN) 


TABLE  C— 1 


(cont '  d) 

Countries  and  Their  Associated  Datums 


COUNTRY /GEOGRAPHIC  AREA 

Eritrea 

Espirito  Santo  Island 
Ethiopia 

Faial  Island  (Azores) 

Federal  Republic  of 
Germany  (before  1990) 

Federal  States  of 
Micronesia 

Finland 

Flores  Island  (Azores) 
Florida 

France 

Gabon 

Ghana 


DATUM 

Massawa 

Santo  (DOS)  1965 
Adindan 

Graciosa  Base  SW  1948 
European  1950 

Kusaie  1951 

European  1950,  1979 

Observatorio  Meteorologico 

Cape  Canaveral 

North  American  1927,  1983 

European  1950 

M'poraloko 

Leigon 


1939 


Gibraltar 


European  1950 


TABLE  C-l  (cont’d) 

Countries  and  Their  Associated  Datums 


Gizo  Island  (New 
Georia  Islands) 


DOS  1968 


Graciosa  Island 


Graciosa  Base  SW 


Grand  Canyon 


North  American  1950 


Greece 


European  1950 


Greenland  (Hayes 
Peninsula) 


North  American  1927 


Guadacanal  Island 


Gux  1  Astro 


Guam 


Guam  1963 


Guatemala 


North  American  1927,  1983 


Guinea-Bissau 


Bissau 


Guinea 


Dabola 


Guyana 


Provisional  South  American  1956 
South  American  1969 


Hawaii 


Old  Hawaiian 
North  American  1983 


Honduras 


North  American  1927,  1983 


TABLE  C-l 


(cont 'd) 

Countries  and  Their  Associated  Datums 


COUNTRY / GEOGRAP H IC  AREA 

E&XIZM 

Hong  Kong 

Hong  Kong  1963 

Iceland 

Hjorsey  1955 

India 

Indian  . 

Iran 

European  1950 

Iraq 

European  1950 

Ireland 

European  1950,  Ireland  1965 

Isle  of  Man 

Ordnance  Survey  of  Great  Britain 
1936 

Israel 

European  1950 

Italy 

European  1950 

Iwo  Jima 

Astro  Beacon  "E" 

Jamaica 

North  American  1927 

Japan 

Tokyo 

Johnston  Island 

Johnston  Island  1961 

Jordan 

European  1950 

Kalimantan  Island 
(Indonesia) 

Gunung  Segara 

TABLE  c-l  (cont'd) 

Countries  and  Their  Associated  Datums 


PNTRY / GEOGRAPHIC  AREA 


DATUM 


Kauai 


Old  Hawaiian 
North  American  1983 


Kenya 

Kerguelen  Island 
Korea 

Kuwait 

Lebanon 

Leeward  Islands 

Liberia 
Luxembourg 
Madagascar 
Madeira  Islands 
Mahe  Island 


Arc  1960 

Kerguelen  Island  1949 
Tokyo 

European  1950 

European  1950 

Fort  Thomas  1955 

Antigua  Island  Astro  1943 

Montserrat  Island  Astro  1958 

Liberia  1964 

European  1950 

Tananarive  Observatory 

Porto  Santo  1936 

Mahe  1971 


Malawi 


Arc  1950 


TABLE  C-l 


(cont 'd) 

Countries  and  Their  Associated  Datums 


>?NTRY/ GEOGRAPHIC  AREA 


DATUM 


Malaysia 

Maldives,  Republic  of 
Mali 

Malta 

Manitoba 

Marcus  Islands 

Marshall  Islands 

Mascarene  Island 

Masirah  Island 
(Oman) 

Maui 

Mexico 

Micronesia 
Midway  Island 


Timbalai  1948 

Gan 

Adindan 
European  1950 
North  American  1927,  1983 
Astronomic  Station  1952 

Wake  Eniwetok  1960 

Reunion 

Nahrwan 

Old  Hawaiian 
North  American  1983 

North  American  1927,  1983 

Kusaie  1951 
Midway  Astro  1961 


Mindanao  Island 


Luzon 


TABLE  C-l  (cont'd) 

Countries  and  Their  Associated  Datums 


Montserrat 

Morocco 

Namibia 

Nepal 

Netherlands 

Nevis 

New  Brunswick 

New  Foundland 

New  Zealand 
Nicaragua 
Niger 
Nigeria 

Northern  Ireland 

Northwest  Territories 
Norway 


cam* 

Montserrat  Island  Astro  1958 

Merchich 

Schwarzeck 

Indian 

European  1950,  1979 
Fort  Thomas  1955 
North  American  1927,  1983 
North  American  1927,  1983 
Geodetic  Datum  1949 

North  American  1927 

Point  58 

Minna 

Ireland  1965 

North  American  1927,  1983 
European  1950,  1979 


TABLE  C-l  (cont'd) 

Countries  and  Their  Associated  Datums 


»?NTRY/gSQSRAF.flIC.JmEA 

datum 

Nova  Scotia 

North  American  1927,  1983 

Oahu 

Old  Hawaiian 

North  American  1983 

Okinawa 

1 

Tokyo 

Oman 

Oman 

Ontario 

North  American  1927,  1983 

Paraguay 

Chua  Astro,  South  American  1969 

Peru 

Provisional  South  American  1956 

South  American  1969 

Philippines 

Luzon 

Phoenix  Islands 

Canton  Astro  1966 

Pico  Island  (Azores) 

Graciosa  Base  SW 

Pitcairn  Island 

Pitcairn  Astro  1967 

Porto  Santo  Island 

Porto  Santo  1936 

Portugal 

European  1950 

/ 


Puerto  Rico 


Puerto  Rico 


TABLE  C-l  (cont'd) 

Countries  and  Their  Associated  Datums 


COUNTRY/GEOGRAPHIC  AREA 

DAKZM 

Qatar 

Qatar  National 

Quebec 

North  American  1927,  1983 

Salvage  Islands 

Selvagem  Grade  1938 

San  Salvador  Island 

North  American  1927 

Santa  Maria  Islands 

(Azores) 

Sao  Braz 

Sao  Jorge  (Azores) 

Graciosa  Base  SW 

Sao  Miguel  Island 

Sao  Braz 

Sarawak  and  Sabah 

Timbalai  1948 

Sardinia 

European  1950,  Rome  1940 

Saskatchewan 

North  American  1927,  1983 

Saudi  Arabia 

Ain  El  Abd  1970,  European  1950, 

Nahrwan 

Scotland 

European  1950,  Ordnance  Survey  of 

Great  Britain  1936 

Senegal 

Adindan 

Shetland  Islands 


European  1950,  Ordnance  Survey  of 
Great  Britain  1936 


TABLE  C-l 


(cont 'd) 

Countries  and  Their  Associated  Datums 


Sicily 


European  1950 


Singapore 


Kertau  1948,  South  Asia 


Somalia 


Afgooye 


South  Africa 


Cape 


South  Chile 


Provisional  South  Chilean  1963 


South  Geogia  Islands 


ISTS  Astro  1968 


South  Greenland 


Qornoq 


Spain 


European  1950,  1979 


Sri  Lanka 


Kandawala 


St.  Helena  Island 


Astro  Dos  71/4 


St.  Kitts 


Fort  Thomas  1955 


Sudan 


Adindan 


Sumatra  Island 
(Indonesia) 


Djakarta  (Botavia) 


Surinam 


Zander i j 


Swaziland 


Arc  1950 


TABLE  C-l 


(cont 1 d) 

Countries  and  Their  Associated  Datums 


Sweden 


European  1950,  1979 


Switzerland 


European  1950,  1979 


Syria 


European  1950,  1979 


Taiwan 


Hu-Tzu-Shan 


Tanzania 


Arc  1960 


Tasmania  Island 


Australian  Geodetic  1966,  1984 


Terceira  Islands 
(Azores) 


Graciosa  Base  SW 


Tern  Island 


Astro/Tern  Island  (Frig)  1961 


Thailand 


Indian 


Trinidad  and  Tobago 


Naparima,  BWI;  South  American  1969 


Tristan  da  Cunha 


Tristan  Astro 


Tunisia 


Carthage 


Turk  and  Caicos 
Islands 

United  Arab  Emirates 


North  American  1927 


Nahrwan 


Uruguay 


Yacare 


TABLE  C-l  (cont'd) 

Countries  and  Their  Associated  Datums 


Venezuela 


Vietnam 

Virgin  Islands 

Viti  Levu  Island 
(Fiji  Island) 

Wales 


Wake  Atoll 


Western  Europe 
West  Malaysia 
Western  United  States 


Yukon 


Zaire 


Zambia 


Provisional  South  American  1956 
South  American  1969 

Indian 

Puerto  Rico 

Indian 


Ordnance  Survey  of 
Great  Britain  1936 

Wake  Island  Astro  1952 

European  1950 

Kertau  1948 

North  American  1927,  1983 
North  American  1927,  1983 
Arc  1950 


Arc  1950 


Zimbabwe 


Arc  1950 


10.0  GENERAL. 


APPENDIX  D 

WFS  84  Geoid  Separation  Computation 


10.1  Scope.  This  appendix  provides  the  bilinear  method  for 
interpolating  geoid  separation  and  WGS  84  Geoid  Separations  on  a 
10°  x  10°  grid.  These  separations  are  expressed  in  meters. 

20.  APPLICABLE  DOCUMENTS. 

This  section  is  not  applicable  to  this  appendix. 

30.  GENERAL  REQUIREMENTS. 

30.1  The  Bilinear  Interpolation  Method  of  Calculating  WGS  84 
Geoid  Separation  (Nvrr  ^  .  The  bilinear  interpolation  method  uses 

interpolation  in  two  directions  ((j),X)  to  approximate  the  geoid 
separation  Nwaq  S4,  from  four  known  WGS  84  geoid  separations  given 
in  a  look-up  table.  The  bi-linear  interpolation  method  has  the 
form  given  in  equation  set  D.l.  Figure  D.l  shows  the  relationship 
between  the  positions  of  the  known  NWGS  84  values  and  the 
interpolation  point. 

Np  ((j),  X,)  =  a0  +  axX  +  a2Y  +  a3XY  (D.l) 


where 

Np(<j>,X)  =  Geoid  Separation  (N)  to  be  interpolated  at  Point 
P((J)A)  . 

a0  =  N3 
ax  =  N2  -  Nx 
a2  =  N4  -  N3 
ci^  =  Nx  +  N3  —  N2 

_  ~ 

(X2  -  Xj 

Y  =  (fr  - 

($2  -  ^1) 


145 


<|)  =  Geodetic  latitude  of  Point  P 
A,  =  Geodetic  longitude  of  Point  P 
Nx,  N2,  N3,  N4  =  Known  geoid  heights  at  grid  points  used  in  the 

interpolation  process. 

♦ 


M>.,M 


Figure  F-l  Coordinate  System  Associated  with 
Geoid  Seperation  Bi-Linear  Interpolation  Method 

30.2  NWGS  84  Tables.  Nwgs  84  values  are  tabulated  on  a  10°  x  10° 
grid  in  Table  D.l.  The  values  in  Table  D.l  are  suitable  for  use 
in  moderate  accuracy  mapping  and  charting  applications.  Qualified 
users  may  obtain  digitally-stored  tables  of  Nwgs  84  values  on 
30'  x  30'  and  1°  x  1°  grids  for  survey  and  high  accuracy  mapping 
and  charting  applications. 
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WGS  84  Geoid  Separations 

(N  =  M  =  180,  10°  x  10°  Grid,  Units  =  meters 
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WGS  84  Geoid  Separations 

(N  =  M  =  180,  10°  x  10°  Grid,  Units  =  meters) 
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APPENDIX  E 


Seven-Parameter  Geometric  Datum  Shifts 


10 .  GENERAL 

10.1  Scope.  This  appendix  provides  the  data  needed  to  perform 
seven-parameter  geometric  datum  shifts  between  two  local  geodetic 
systems  and  WGS  84. 

20.  APPLICABLE  DOCUMENTS. 

This  section  is  not  applicable  to  this  appendix. 

30.  GENERAL  REQUIREMENTS 

30.1  Application.  Table  E-l  lists  the  coverage  area,  ellipsoid 
applicable  to  the  local  geodetic  system,  ellipsoid  center  shifts 
(AX,  AY,  AZ),  axes  rotations  (£,  \|C,  CO),  and  scale  factor  change 
(AS)  for  use  in  seven-parameter  geometric  datum  shifts  between 
two  local  geodetic  systems  and  WGS  84.  The  root-mean 
square  differences  between  Doppler-derived  WGS  84  geodetic 
coordinates  and  WGS  84  coordinates,  computed  from  local  system 
coordinates  using  the  seven-parameter  model,  are  provided  as  an 
indication  of  the  datum  shift  consistency. 
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TABLE  E.l 


Seven-Parameter  Local  Geodetic  System  to  WGS  84 
Datum  Transformations 


DATUM:  EUROPEAN  1950  (EUR-M) 

COVERAGE  AREA:  Western  Europe 
ELLIPSOID:  International 

AX  =  -102  m  AY  =  -102  m  AZ  =  -129  m 

8  =  0.413"  \|/  =  -0.184"  CO  =  0.385"  AS  =  2.4664  x  10'6 

RMS  DIFFERENCES  :  WGS  84  surveyed,  minus  computed  coordinates 

A(|)  =  2  m  A  A,  =  3  m  AH  =  2  m 


DATUM:  ORDNANCE  SURVEY  OF  GREAT  BRITAIN  1936 

COVERAGE  AREA:  England,  Isle  of  Man,  Scotland,  Shetland  Islands, 

and  Wales 

ELLIPSOID:  Airy 

AX  =  446  m  AY  =  -99  m  AZ  =  544  m 

8  =  -0.945"  \|/  =  -0.261"  (0  =  -0.435"  AS  =  -20.8927  x  10'6 

RMS  DIFFERENCES  :  WGS  84  surveyed  minus  computed  coordinates 
A(j)  =  2  m  AA  =  2  m  AH  =  1  m 
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APPENDIX  F 


Old  Hawaiian  Datum  with  International  Ellipsoid 


10.  GENERAL 

10.1  Scope .  This  appendix  provides  the  information  needed  to 
shift  coordinates  between  the  Clarke  1866  and  International  in  the 
Old  Hawaiian  datum. 

20.  APPLICABLE  DOCUMENTS 

DMA  SGG  MEMORANDOM  FOR  RECORD 
SUBJECT:  Old  Hawaiian  Datum  with 
International  Ellipsoid 
. 13  June  1990 

Army  Map  Service  GEODETIC  MEMO 
No.  687,  SUBJECT:  Hawaiian  Islands, 

UTM  Coordinates  and  International  Spheroid 
Geographic  Positions.  P.O.#  23832-002, 

29  November  1950. 

30.  GENERAL  REQUIREMENTS 

30.1  Background .  Positions  on  the  Old  Hawaiian  datum  may  be  on 
either  the  Clarke  1866  ellipsoid  or  the  International  ellipsoid. 
Users  of  Old  Hawaiian  datum  positions  must  determine  which 
ellipsoid  applies.  Datum  shift  parameters  from  Old  Hawaiian  datum 
to  WGS  84  have  been  determined  only  for  the  Clarke  1866  ellipsoid. 
Old  Hawaiian  datum  coordinates  on  the  international  ellipsoid  must 
be  converted  to  the  Clarke  1866  ellipsoid  before  being  shifted  to 
WGS  84.  The  conversion  between  International  and  Clarke  1866 
ellipsoids  has  been  defined  in  terms  of  UTM  coordinates.  Geodetic 
latitude  and  longitude  must  be  converted  to  UTM  before  changing 
ellipsoids.  Height  conversions  should  not  be  attempted  as  no 
conversion  method  is  available. 

30.2  Application.  To  convert  Old  Hawaiian  datum  UTM  coordinates 
on  the  International  ellipsoid  (EINT,  Nint)  to  UTM  coordinates  on 
the  Clarke  1866  ellipsoid  (EC66/  NC66)  use: 

In  UTM  Zone  4 

EC66  =  Eint  -  3  (1.1) 

Nc66  =  Nint  -  169 

In  UTM  Zone  5 

Ec66  =  Ejnt  +13  (1.2) 

Nc66  =  Nint  ~  169 
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To  convert  Old  Hawaiian  datum  UTM  coordinates  on  the  Clarke 
1866  ellipsoid  to  the  International  ellipsoid  use: 


In  UTM  Zone  4 

EINt  =  Ec66  +  3  (1.3) 

Nint  =  Nc66  +  169 

In  UTM  Zone  5 

Eint  =  Ec66  ~  13  (1.4) 

Nint  =  Nc66  +  169 

30.3  Example .  Convert  geodetic  survey  coordinates  on  the  Old 
Hawaiian  datum.  International  ellipsoid,  to  WGS  84  geodetic 
coordinates . 

<J)INT  =  21°  19'  37 . 431  "N 
A,int  =  157°  58'  25. 610 "W 

This  position  is  on  the  island  of  Oahu. 

a.  Convert  <|)int/  ^-int  to  UTM: 

Zone  =  4 

Eint  =  606,428 
Nint  =  2,358,722 

b.  Convert  Eint/  Nint  to  the  Clarke  1866  ellipsoid  using 
equation  set  1.1,  for  zone  4. 

Ec66  =  606,425 
Nc66  =  2,358,553 

c.  Convert  EC66»  NC66  to  geodetic  coordinates  on  the 
Clarke  1866  ellipsoid. 

(|>C66  =  21°  19'  37 . 425  "N 
A<c66  =  157°  58  '  25 . 631  "W 

d.  AX,  AY,  AZ  shift  parameters  for  datum  identifier  code 
OHA-D  should  be  used  for  shifting  geodetic  survey 
coordinates  on  Oahu. 

e.  Shift  (]>c66 >  ^c66  to  WGS  84  using,  for  example,  the 
three-step  method. 
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<1>wgs  84  =  21°  19*  26 . 069  "N 
Kgs  84  =  157°  58 1  15.766"W 
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